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Method and Apparatus for Routing Sets of Nets 

CLAIM OF BENEFIT TO PRIOR APPLICATION 

This patent application is a continuation of the U.S. Patent Application 10/215,896, filed 
8/9/2002, which is incorporated herein by reference. This application also claims the benefit of the 
U.S. Provisional Patent Application 60/396,571, filed 7/15/2002; U.S. Provisional Patent 
Application 60/388,518, filed 6/12/2002; U.S. Provisional Patent Application 60/385,975, filed 
6/4/2002; U.S. Provisional Patent Application 60/351,459, filed 1/22/2002. 

FIELD OF THE INVENTION 

The invention is directed towards method and apparatus for routing sets of nets. 

BACKGROUND OF THE INVENTION 

An integrated circuit ("IC") is a device that includes many electronic components (e.g., 
transistors, resistors, diodes, etc.). These components are often interconnected to form multiple 
circuit components (e.g., gates, cells, memory units, arithmetic units, controllers, decoders, etc.) on 
the IC. The electronic and circuit components of IC's are jointly referred to below as "components." 

An IC also includes multiple layers of wiring ("wiring layers") that interconnect its electronic 
and circuit components. For instance, many IC's are currently fabricated with metal or polysilicon 
wiring layers (collectively referred to below as "metal layers," "interconnect layers" or "wiring 
layers") that interconnect its electronic and circuit components. One common fabrication model uses 
five metal layers. In theory, the wiring on the metal layers can be all-angle wiring (i.e., the wiring can 
be in any arbitrary direction). Such all-angle wiring is commonly referred to as Euclidean wiring. In 
practice, however, each metal layer typically has a preferred wiring direction, and the preferred 
direction alternates between successive metal layers. IC designs often penalize non-preferred- 
direction wiring on a layer. 

Many IC's use the Manhattan wiring model, which specifies alternating layers of preferred- 
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direction horizontal and vertical wiring. In this wiring model, the majority of the wires can only 
make 90° turns. However, occasional diagonal jogs are sometimes allowed on the preferred 
horizontal and vertical layers. 

Design engineers design IC's by transforming circuit description of the IC's into geometric 
descriptions, called layouts. To create layouts, design engineers typically use electronic design 
automation ("EDA") applications. These applications provide sets of computer-based tools for 
creating, editing, and analyzing IC design layouts. 

EDA applications create layouts by using geometric shapes that represent different materials 
and devices on IC's. For instance, EDA tools commonly use rectangular lines to represent the wire 
segments that interconnect the IC components. These tools also represent electronic and circuit IC 
components as geometric objects with varying shapes and sizes. 

Also, in this document, the phrase "circuit module" refers to the geometric representation of 
an electronic or circuit IC component by an EDA application. EDA applications typically illustrate 
circuit modules with pins on their sides. These pins connect to the interconnect lines. 

A net is typically defined as a collection of pins that need to be electrically connected. A list 
of all or some of the nets in a layout is referred to as a netlist. In other words, a netlist specifies a 
group of nets, which, in turn, specify the required interconnections between a set of pins. 

The IC design process entails various operations. Some of the physical-design operations that 
EDA applications commonly perform to obtain the IC layouts are: (1) circuit partitioning, which 
partitions a circuit if the circuit is too large for a single chip; (2) floor planning, which finds the 
alignment and relative orientation of the circuit modules; (3) placement, which determines more 
precisely the positions of the circuit modules; (4) routing, which completes the interconnects 
between the circuit modules; and (5) verification, which checks the layout to ensure that it meets 
design and functional requirements. 

Routing is a key operation in the physical design cycle. It is generally divided into two 



phases: global routing and detailed routing. For each net, global routing generates a "loose" route for 
the interconnect lines that are to connect the pins of the net. After global routes have been created, 
the detailed routing creates specific individual routing paths for each net. 

While some commercial routers today might allow an occasional diagonal jog, these routers 
5 do not typically explore diagonal routing paths consistently when they are specifying the routing 
geometries of the interconnect lines. This, in turn, increases the total wirelength {i.e., total length of 
interconnect lines) needed to connect the nets in the layout. 

In addition, routers today are mostly gridded. The manufacturing processes for designing ICs 
specify a manufacturing grid that specifies manufacturable resolution. The boundary of all circuit 
10 elements is defined by the straight-line connections between adjacent manufacturing points. Gridded 
routers typically define arbitrary grids of intersecting lines to specify the available locations for 
routing interconnects. These arbitrary grids are often much coarser than the manufacturing grids 
(e.g., they are typically line-to-via spacing). Consequently, they arbitrarily limit the locations of 
interconnect lines and impose arbitrary spacing between the items in the layout. These arbitrary 
15 limits increase the size and efficiency of a design. The routing grids also discourage using arbitrary 
widths or spacing for interconnect lines. 

Furthermore, existing routers primarily utilize preferred-direction wiring to route their 
designs. Many IC layouts are designed by penalizing the use of interconnect lines in each particular 
layer when the interconnect lines are not in the preferred wiring direction of the particular layer. 
20 Such preferred-direction wiring leads to IC layouts and ICs that have most of their interconnect lines 
and wiring on each of their metal layers traverse in the same direction. Such IC layouts and ICs do 
not efficiently use the available spacing on the interconnect layers, and this adversely affects the size 
and efficiency of the layouts and the ICs. 
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SUMMARY OF THE INVENTION 

Some embodiments of the invention provide a method of routing a set of nets. The method 
specifies a first order for the set of nets. It then routes the nets according to the specified first order. 
The method then specifies a second order for the set of nets, where the second order has the fewest 
5 possible number of differences with the first order. The method then routes the nets according to the 
specified second order. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The novel features of the invention are set forth in the appended claims. However, for 
purpose of explanation, several embodiments of the invention are set forth in the following figures. 

Figure 1 illustrates a wiring model of some embodiments of the invention. 
5 Figure 2 illustrates a process that is used by some embodiments to produce ICs that utilize 

NPD-wiring models. 

Figure 3 presents a conceptual illustration of a routing process used by some embodiments of 
the invention. 

Figures 4 and 5 illustrate two layers of Gcells that have been combined to produce a sub- 

10 region. 

Figure 6 illustrates a portion of a layer that has been triangulated into four triangular faces, 
each of which has a space, three nodes, and three edges, while Figure 7 provides examples of 
topological routes that are formed by via nodes, Steiners, walls, joints, and nodes. 

Figure 8 illustrates a process that provides the overall flow of the Q* topological engine in 
15 some embodiments. 

Figure 9 illustrates a triangulation process that is used in some embodiment of the invention 
by the Q* topological routing engine. 

Figure 10 illustrates the layout of Figure 4 after nodes have been defined at each sub-region 
corner and at each port or obstacle geometry point. 
20 Figure 11 illustrates a triangulation technique. 

Figures 12 and 13 illustrate why maximizing the minimal angle of decomposed triangles 
improves the likelihood that generated topological routes can be geometrized. 

Figures 14 and 15 illustrate one manner for performing an edge-flipping operation. 

Figure 16 provides a pictorial example of a constraining operation. 
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Figure 17 illustrates a hole-identification process that can be used to try to specify holes 
within triangles. 

Figure 18 illustrates an example of how the layout of Figure 4 might look after triangulation. 
Figure 19 and 20 illustrate how some embodiments calculate edge capacities. 
5 Figure 21 illustrates the overall flow of a solving engine of some embodiments of the 

invention. 

Figure 22 illustrates a route generation process of some embodiments. 

Figures 23-30 provide an example of an A* path search that uses such an F cost. 

Figure 31 illustrates an example of a line PLF that has four line segments. 
10 Figure 32 illustrates an example of a surface PLF. 

Figure 33 presents an example of filtering a first filtered PLF by a second filter PLF, where 
both PLFs are defined across a line. 

Figure 34 illustrates the minimum PLF for the two PLFs of Figure 33. 

Figure 35 illustrates a Q* path search process of some embodiments of the invention. 
15 Figures 36-39 provide examples for describing the process of Figure 35. 

Figure 40 illustrates a process that propagates a PLF that is defined over a point to a line or a 
surface. 

Figures 41 and 42 illustrate examples of propagating a PLF from a point P to a line L and to 
a surface S. 

20 Figure 43 illustrates a process for propagating a PLF from a line to another line or from a 

surface to a line, and Figure 44 and 45 provides examples for describing this process. 

Figures 46-51 illustrate how some embodiments identify the propagation vectors that 
emanate from the knot locations of a line PLF or surface PLF. 

Figure 52 illustrates a process for propagating a G PLF from a line to a surface. 

25 Figure 53 illustrates an example for propagating a G PLF from a line to a surface. 
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Figure 54 illustrates a process for propagating a PLF from a line to a point or from a surface 
to a point. Figures 55 and 56 describe the process of Figure 54. 

Figure 57 presents an example that illustrates an expansion from a start surface to a 
destination surface. 

5 Figures 58 A and 58B illustrate how to compute the flow across an edge after a potential 

expansion. 

Figures 59 and 61 illustrate processes that generate PLFs that express costs of expansions to 
an edge or a hole, while Figures 60 and 62 present examples that describe these processes. 
Figures 63-65 illustrate how to add two surface PLFs. 
10 Figure 66 illustrates a process that performs filtering and minimum operations for an 

expansion to a line. 

Figures 67 and 68 illustrate how to filter two surface PLFs. 

Figures 69 and 70 illustrate how some embodiments compute an H function. 
Figure 71 illustrates a process that specifies a topological path. 
15 Figures 72-76 illustrate data structure that defines a face, an edge, a node, an edge item, and 

a face item in some embodiments of the invention. 

Figures 77 and 78 illustrate an example of topological routes. 

Figure 79 illustrates a process that provides the overall flow of an IDA* topological engine 
in some embodiments of the invention. 
20 Figure 80 pictorially illustrates an example of a solving engine's IDA*-searching operation 

for a set of three nets. 

Figure 81 illustrates a more detailed process used by the solving engine in some 

embodiments of the invention. 

Figure 82 illustrates a process that the solving engine uses to generate topological routes. 

25 Figure 83A and 83B illustrate a process for inserting Steiner-tree face items in face. 
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Figure 84 illustrates a process for generating paths between one or more sources and one or 
more targets for a selected pin-pair. 

Figures 85-89 illustrate possible expansions from edge items, nodes, and face items. 
Figure 90 illustrates three types of legality checking. 
5 Figures 91 and 92 illustrate two via-checking processes. 

Figures 93-105 present several examples that illustrate the process of Figure 91. 
Figure 106 illustrates a process for computing the cheapest-route cost for a net. 
Figure 107 illustrates a computer system used in some embodiments. 



DETAILED DESCRIPTION OF THE INVENTION 

In the following description, numerous details are set forth for purpose of explanation. 
However, one of ordinary skill in the art will realize that the invention may be practiced without the 
use of these specific details. In other instances, well-known structures and devices are shown in 
5 block diagram form in order not to obscure the description of the invention with unnecessary detail. 
I. NON-PREFERRED DIRECTION ARCHITECTURE 

Some embodiments of the invention utilize non-preferred-direction ("NPD") wiring models 
for designing IC layouts. An NPD wiring model has at least one NPD interconnect layer that has 
more than one preferred routing direction. Specifically, in such a wiring model, each NPD 

10 interconnect layer has at least two wiring directions that are equally preferable with respect to one 
another, and that are as preferable as or more preferable than the other wiring directions on that layer. 

An NPD router that uses an NPD wiring model generates IC layouts with one or more NPD 
interconnect layers, with each having two or more preferred wiring directions. In some embodiments, 
such a router does not impose arbitrary penalty costs for routing (i.e., does not penalize wiring) in the 

15 two or more preferred directions of an NPD interconnect layer of its layout. Rather, such a router 
costs wiring {i.e., interconnect lines) in the preferred directions proportionately to the metric costs 
that they introduce in the design. The metric cost is a non-arbitrary cost that is based on putative chip 
performance, manufacturing considerations, or other real-world design consideration, and not an 
artifical characteristic of the router. The metric costs can be based on a number of properties of the 

20 wiring, such as resistance, delay, manufacturing yield, etc. A router can account for these metric 
costs by using costing functions that are based on the interconnect-line attributes, such as length, 
area, etc. 

For instance, horizontal and +45° wiring directions might be the preferred wiring directions 

of an interconnect layer. On such a layer, an NPD router would not impose an arbitrary penalty cost 

25 on routing in either the horizontal or +45° directions, in order to deter routes in one of these 
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directions. Instead, some embodiments might cost horizontal and +45° lines based solely on the 
length of the lines in these directions, while other embodiments might cost horizontal and 445° lines 
based solely on the area of the lines in these directions. Costing the lines based on their areas more 
accurately estimates the metric resistance and delay costs of the lines when the widths of the lines 
have to be different in different preferred directions. For instance, manufacturing design rules might 
require the preferred direction +45° wires to be wider than preferred direction horizontal wires on a 
layer. Such 45° lines would be more expensive than horizontal lines if area were used to cost the 
routes. This additional expense, however, is not an artificial penalty cost, but rather is an actual cost 
that relates directly to the metric cost that the 445° lines introduce into the design. It should be noted 
that such 45° lines would be less expensive than horizontal lines if delay were used as the primary 
metric for costing the routes. 

Figure 1 illustrates a wiring model 100 of some embodiments of the invention. This wiring 
model has five interconnect layers 105-125. In this five-layer model, none of the layers has a single 
preferred wiring direction. Instead, each layer allows 4 different preferred directions of wiring. 
Specifically, as illustrated by the top view 130 of the fifth interconnect layer 125, each layer of the 
wiring model 100 can have horizontal, vertical, and ±45° diagonal interconnect lines. Each layer is 
an octilinear layer, as it allows interconnect lines to traverse in eight separate vector directions from 
any given point. 

In some embodiments, an interconnect line is "diagonal" if it forms an angle other than zero 
or ninety degrees with respect to the layout's Cartesian coordinate axes, which typically align with 
the layout's boundary and/or the boundary of the layout's expected IC. On the other hand, an 
interconnect line is "horizontal" or "vertical" if it forms an angle of 0° or 90° with respect to one of 
the coordinate axes of the layout. In the wiring model of Figure 1, (1) the horizontal interconnect 
lines are parallel (i.e., are at 0°) to the x-axis, which is defined to be parallel to the base of the layout, 
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(2) the vertical interconnect lines are parallel to the y-axis, which is defined to be parallel to the 
height and perpendicular (Le. 9 is at 90°) to the base of the layout, (3) the +45° diagonal interconnect 
lines are at +45° with respect to the x-axis, and (4) the -45° diagonal interconnect lines are at -45° 
with respect to the x-axis. 

Other embodiments of the invention use different NPD wiring models. For instance, some 
embodiments of the invention's NPD wiring model include only diagonal interconnect lines, other 
embodiments use only horizontal and vertical interconnect lines, and yet other embodiments use 
diagonal interconnect lines with either horizontal or vertical interconnect lines but not both. Also, 
some embodiments use non-45° diagonal wiring. For example, some embodiments use horizontal, 
vertical and ±120° diagonal interconnect lines. 

In addition, some embodiments have more than five layers, while other embodiments have 
fewer. Some embodiments also assign a single preferred direction for some of the layers, while 
allowing other layers not to have a single preferred wiring direction. For instance, some 
embodiments have preferred-direction Manhattan wiring for the first three layers (e.g., horizontal 
preferred-direction wiring for the first layer, vertical preferred-direction wiring for the second layer, 
and horizontal preferred-direction wiring for the third layer), and NPD wiring for the fourth and fifth 
layers. 

By generating IC layouts with one or more NPD interconnect layers, some embodiments 
fabricate ICs that have NPD wiring for one or more of their metal layers (i.e., fabricate IC's with one 
or more NPD metal layers, with each having two or more preferred wiring directions). For instance, 
the wiring model 100 of Figure 1 can be used to generate an IC layout with five NPD-interconnect 
layers. Such a layout can then be used to generate a five metal layer IC, where each of the metal 
layers has four equally preferable wiring directions. 

Figure 2 illustrates a process 200 that is used by some embodiments to produce ICs that 

utilize NPD-wiring models. This process initially describes (at 205) the specification of an IC. The 
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IC specification is a high-level representation that considers various factors, such as performance, 
functionality, physical dimension, etc. The process then specifies (at 210) the circuit design of the IC. 
This operation typically entails one or more sub-operations, such as specifying the ICs architectural 
design, functional design, logic design, etc. Numerous EDA tools exist for assisting in the creation of 
the system's circuit design. 

After specifying the circuit design, the process specifies (at 215) a physical design layout for 
the IC. According to the invention, the process specifies the physical design layout based on an 
NPD-wiring model. The generation of physical design layouts that use the NPD-wiring model will be 
described below. At 215, the physical-design operation might result in the partitioning of the IC into 
multiple ICs (Le., might result in the generation of one or more layouts that represent several ICs). 
Also, at 215, the process verifies the generated physical design layouts. Verification typically entails 
various operations, such as design-rule check, extraction, etc. These operations are to ensure that the 
generated layouts specify designs that meet the fabrication rules and perform the desired 
functionalities. 

Once physical design and verification are completed, the process converts (at 220) the 
physical design layout into one or more sets of masks for one or more ICs. Each set of masks 
typically includes one mask for each wiring layer. Currently, each mask is typically a photo- 
lithographic mask that identifies spaces on a semiconductor wafer where certain materials need to be 
deposited, diffused and/or removed. Accordingly, at 225, the process uses the masks to fabricate one 
or more ICs. Specifically, the fabrication operation uses the masks to create devices on ICs and to 
define wiring that connect the devices. When the physical design layout uses an NPD-wiring model, 
the masks define NPD wiring on one or more layers of the ICs. After fabrication, the process 200 
tests (at 230) the fabricated ICs. The process then (at 235) packages any IC that passes its tests, and 
then ends. 

Some embodiments fabricate ICs that have at least one wiring layer that does not have a 
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wiring direction with more that 50% of the wiring on that layer. Other embodiments fabricate ICs 
that have at least one wiring layer that does not have a wiring direction with more than 70% of the 
wiring on that layer. For instance; one such IC might predominantly include horizontal and vertical 
direction wiring on the wiring layer that does not have any wiring direction with more than 70% of 
the wiring. 

n. GRIDLESS ARCHITECTURE 

A gridless routing process is described below. This routing process is gridless as it does not 
require the interconnect lines to be positioned with respect to any grid that is coarser than the 
manufacturing grid, if any. In other words, the only grid that the interconnect lines might have to be 
aligned with is the manufacturing grid. The generated gridless layouts can be used, in turn, to 
fabricate ICs that have their metal lines aligned with the finer manufacturing grids instead of coarser 
non-manufacturing grids. 

The gridless routing process described below generates gridless NPD octilinear layouts. 
However, one of ordinary skill will realize that this routing process can be used to generate other 
gridless layouts. In addition, even though a gridless NPD routing process is described below, some 
embodiments can use the invention's NPD routing in a gridded router that generates gridded NPD 
layouts. 

III. NPD AND GRIDLESS ROUTING 
A. Conceptual Flow, 

Some embodiments generate gridless NPD layouts by using a detailed routing technique that 
does not specify a preferred wiring direction for any of its interconnect layers. The detail-routing 
embodiments described below use the NPD wiring model 100 of Figure 1. However, one of ordinary 
skill will realize that other embodiments of the invention use different NPD wiring models. 

In the embodiments described below, the detailed routing is performed after a global routing 
stage, which (1) partitions the routing region into global routing cells ("Gcells") and (2) defines, for 
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each net, global routing paths that connect the Gcells containing the net's pins. One hierarchical 
global routing approach recursively divides the routing region into smaller sub-regions, and defines 
routing paths at each hierarchical level until reaching the lowest recursive level's sub-regions, which 
are the Gcells. Another global routing approach flatly divides the routing region into numerous 

5 Gcells and then defines the routing paths between the Gcells. Under either approach, the global 
router can use either an NPD wiring model or a preferred-direction wiring model. 

Figure 3 presents a conceptual illustration of a detail-routing process 300 used by some 
embodiments of the invention. This routing process defines detail routes for nets within a region of 
the IC layout. This region can be the entire IC layout, or a portion of this layout. As shown in this 

10 figure, this process initially selects (at 305) a sub-region of the IC layout region to detail route. The 
sub-region can be a portion of the IC layout, or the entire IC layout. Several manners for selecting 
such a region will be described below in Section in.B. 

Next, for each particular net in the selected sub-region, the process identifies (at 310) a 
topological route that connects the particular net's routable elements in the sub-region. In the 

15 embodiments described below, a net has two or more pins, a pin can have one or more ports, and 
each port can have one or more geometries. In these embodiments, a net's routable elements are the 
port geometries, and a net is typically routed along one port of each of its pins. One of ordinary skill 
will realize, however, that other embodiments may define the routable elements of the nets 
differently. 

20 A topological route is a route that is defined in terms of its relation to other layout items, such 

as pins, obstacles, boundaries, and/or other topological routes of other nets. As such, a topological 
route provides a general plan for how to route a net, without necessarily providing a specific 
geometric path to do so. One topological route represents a set of diffeomorphic geometric routes 
(/.«., a set of geometric routes that can be morphed into one another through a continuous sequence 

25 of perturbations without changing the route's path relative to any other pin, path or obstacle). A 
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geometric route is one explicit realization of a topological route. A geometric route is defined in 
terms of exact coordinates that define the route as it travels through the interconnect layers. Several 
manners for identifying topological routes for each net within the selected sub-region will be 
described below. 

5 After 310, the process determines (at 315) whether the identified topological routes identified 

at 310 are geometrically routable (Le. t whether there exists a design-rule-correct geometric route for 
each identified topological route). If so, the process transitions to 320, which will be described 
below. Otherwise, if the process determines (at 315) that the identified topological routes for some of 
the nets are not routable, it initially directs the topological router to generate additional topological 

10 routes that are more likely to have design rule-correct geometric routes. If the topological router 
repeatedly fails to generate geometrically routable topological routes, the detail-routing process flags 
one or more nets as unroutable, re-defines topological routes for some or all the nets, and then 
transitions to 320. Performing the routability checking is described in U.S. Patent Application 
. "Method and Apparatus for Routing Nets in an Integrated Circuit Layout", with the serial number 

15 10/215,563 and filed on 8/9/2002. This application is incorporated in the present application by 
reference. 

At 320, the process generates these geometric routes and stores these routes in a detail- 
routing storage structure (such as a database). Generating geometric routes from topological routes is 
described in the above-mentioned patent application. At 320, the process also converts the generated 

20 geometric detail routes into global routing paths, which it stores in a global routing storage structure 
(such as a database). This is done just in case the router has to detail route some Gcells again. At 
325, the process then determines whether it has generated detail routes for all the sub-regions of the 
IC region. If not, the process returns to 305 to select another sub-region and to repeat 310-320 to 
compute geometric routes in the newly selected sub-region. Otherwise, the process ends. After 325, 

25 some embodiments might repeat process 300 for certain congested sub-regions in order to alleviate 
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the congestion in these regions, improve wiring quality, or fix violations left by previous attempts. 
B. Repion Selection. 

As mentioned above, the detail-routing process 300 selects (at 305) a sub-region of the IC 
layout region to detail route. In some embodiments of the invention, this selection involves selecting 

5 several contiguous Gcells, and generating a sub-region by combining the selected Gcells. Several 
manners for selecting contiguous Gcells are disclosed in United States Patent Application entitled 
"Method and Apparatus for Generating Multi-Layer Routes" and having serial number 10/076,121. 
One of ordinary skill will realize that other embodiments might select and combine Gcells in order to 
identify a sub-region, but rather might simply specify a portion or the entire IC layout as the sub- 

10 region. 

In generating a sub-region by combining several contiguous Gcells, the process 300 adds any 
virtual pins ("vpins") of the Gcells in the periphery of the sub-region as geometries of the sub-region. 
Virtual pins are artificial pins that are set to account for the global route propagation into Gcells from 
other Gcells or higher-level slots. In the embodiments described below, the virtual pins are 

1 5 represented as single point geometries. 

For example, Figures 4 and 5 illustrate two layers of 16 Gcells that have been combined to 
produce a sub-region 400. For sake of simplicity, this example assumes that the generated sub-region 
only has pins, virtual pins, and obstacles on the two layers shown in Figures 4 and 5. As shown in 
Figure 4, the sub-region 400 includes two obstacles 410 and 415. This sub-region also includes port 

20 geometries of three nets A, B, and C. Again, for sake of simplifying the example, the port geometries 
420-445 and 505 of nets A, B, and C are all from the same ports of their respective nets in the 
example illustrated in Figure 4. As shown in Figures 4 and 5, the obstacle and port geometries can 
have a variety of shapes (e.g., convex or non-convex polygonal shapes, etc.). Also, these geometries 
have horizontal, vertical, and diagonal sides. As shown in Figure 4, three of the twelve Gcells on the 
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periphery include virtual pins 450, 455, 460, and 510 for nets A and C. During the sub-region 
generation, these four virtual pins 450, 455, 460, and 510 are added to the sub-region data structure 
as port geometries of nets A and C. 

The detailed routing process 300 defines a sub-region based on various attributes. These 
attributes include a list of all obstacle and port geometries in the region, and a list (/.e., a netlist) that 
specifies the nets in the sub-region. In the embodiments described below, each net specifies one or 
more pins, each pin refers to one or more ports, and each port refers to one or more geometries. 
Other embodiments might specify a net differently. A sub-region is also defined by information 
about net properties on different layers. For different layers, this information can include minimum 
wire size, minimum spacing, minimum via size, minimum cost per unit length, etc. In some 
embodiments, this information also includes (1) different net widths on different layers, (2) different 
spacing between nets on different layers, and. (3) different spacings between nets and unrelated 
geometries on different layers. 

C. Topological Route Generation. 

In some embodiments, a topological routing engine receives the sub-region defined at 205 
(Le. t receives the problem instance). For each net in the received sub-region, the topological routing 
engine generates (at 310) a topological representation of a route that connects the net's routable 
elements. Different embodiments use different topological routing engines. Below, two topological 
routers are described that can be used as the topological routing engine. One topological routing 
engine is referred to as the Q* topological routing engine, while the other is referred to as the IDA* 
topological routing engine. Both engines are multi-layer topological routers. Accordingly, for each 
net in the sub-region, each engine generates a topological route (i.e. f a topological representation of a 
route) that connects the net's routable elements on one or more layers. In other words, each router 
selects a net and, for the selected net, defines a topological route that connects the selected net's 
routable elements on one or more interconnect layers, before selecting another net for routing. To 
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facilitate its multi-layer approach, each topological router uses vias that are defined topological^. 
These vias are referred to below as topological vias. 

In addition, both topological routing engines can route sets of nets as an ensemble. 
Specifically, the Q* engine selects a set of nets in the sub-region. For the selected set, this engine 
5 might generate several different routing solution sets, where each solution includes a topological 
route for some or all of the nets in the set. From the identified set of solutions, this engine then 
selects the best solution set for the selected set of nets. Similarly, the IDA* routing engine selects a 
set of nets in the sub-region and deterministically traverses through the solution space to identify the 
best possible combination of topological routes for the selected set of nets. Both topological engines 
10 can consider the estimated routing cost of unrouted nets in the selected set while they are generating 
a solution set. These engines consider such estimated costs so that they can terminate examining a 
solution set when an estimated cost exceeds an acceptable threshold. 

The IDA* routing engine uses an IDA* path search engine that identifies paths between sets 
of sources and targets in a depth-first, iterative-deepening manner. The Q* routing engine, on the 
15 other hand, uses a Q* path search process that finds the shortest topological route for each net. The 
Q* path search process also has a powerful shove operation that allows the topological routes of 
later-routed nets to shove the topological routes of earlier-routed nets. In the embodiments described 
below, this operation allows the earlier-routed nets to bend around the topological routes of later- 
routed nets. Accordingly, this operation reduces the routing constraints placed on later-routed nets by 

20 the routes of earlier-routed nets. 

In the embodiments described below, both topological routers are NPD routers that use NPD 
wiring models (such as the wiring model of Figure 1). In some embodiments, both topological 
engines allow nets to have different widths on different layers. They also can impose different 
spacing constraints , between pairs of nets. The spacing constraint for a pair of nets can also be 
25 different on different layers. Both topological engines can also base their topological routes on 
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different wiring models, such as a wiring model that employs only Manhattan lines, one that uses 
Manhattan and ±45° lines, one that uses Manhattan and ±120° lines, etc. 
1. Q* Topological Routing Engine 
a. Overview 

5 The Q* topological engine decomposes each layer of the sub-region that it receives into 

numerous polygons and embeds topological routes into the decomposed region. Different 
embodiments decompose the sub-region into different types of polygons, such as triangles, 
rectangles, etc. The embodiments described below use a triangulation decomposition operation that 
decomposes each layer of the sub-region into numerous triangular faces, each with three nodes and 
10 an edge between each pair of nodes. The decomposed sub-region is at times referred to below as the 
tessellated region. 

The initial triangulation also specifies one space within each face. Figure 6 illustrates a 
portion of a layer that has been triangulated into four triangular faces, each of which has a space, 
three nodes, and three edges. For each layer of the sub-region, the initial triangulation operation adds 

15 a topological graph to the sub-region definition. Each layer's topological graph includes several 
nodes, edges, and faces. In the embodiments described below, nodes are defined at the obstacle- 
geometry vertices, port-geometry vertices, the four corners of the sub-region on each layer, and some 
or all of the vpins. The embodiments described below define edges between certain pairs of nodes in 
order to triangulate the sub-region. An edge can be part of up to two triangular faces. For each edge, 

20 the Q* topological engine stores the edge's capacity and length. It further stores the wire flow (route 
flow) across each edge. 

Different embodiments define edges differently. In the embodiments described below, the Q* 
topological engine defines each edge as an ordered list of topological particles in the tessellated 
region. In these embodiments, the topological particles that can be on an edge (i.e., that can define an 
25 edge) are nodes, joints, and exits, which are described below. 
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Nodes are vertices of the tessellated region, as described above. Two nodes start and end each 
edge's ordered list of topological particles. A node can be part of multiple edges. A joint is a 
topological particle that is defined at the intersection of an edge and a net's route. An exit is a 
topological particle that is specified on an edge between each two non-exit particles on the edge. As 
illustrated in Figure 6, an exit is the only topological particle between an edge's nodes, after the 
initial triangulation before any route has been embedded in the region. Whenever a joint is defined 
on an edge, it divides a previously undivided exit of the edge into two exits. Joints and exits are 

further described below. 

Each triangular face has 3 edges, 3 nodes, and a set of "spaces" representing topological 
regions inside the face. After the initial triangulation, each face has one space within it, as shown in 
Figure 6. When topological routes are embedded across a face, additional spaces are defined in the 
face. Every particle that is within a face or is on an edge of the face is associated with one or more 
spaces of the face. Some particles can be part of multiple spaces. Each space contains its face's nodes 
and exits just after the initial triangulation. 

After the initial triangulation, the Q* engine embeds topological routes in the decomposed 
topological region. The Q* engine defines a topological route as an ordered list of topological 
particles in the decomposed topological region. Each embedded topological route is a graph that has 
several vertices and one or more wire segments that connect some or all the vertices. The wire 
segments are referred to as walls in the embodiments described below. In other words, a wall is a 
wire segment {i.e., an interconnect segment) that connects two vertices in a topological route. In the 
embodiments described below, the topological particles that can serve as vertices of a route are 
nodes, joints, "holes," "via nodes," and "Steiners." 

A hole is a topological particle that represents a potential via between two spaces on two 
different layers. In the embodiments described below, a hole is defined between each pair of spaces 
that are on adjacent layers and that have sufficiently large overlap to accommodate a via. In other 
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embodiments, however, holes might be defined between space pairs that are on different but not 
adjacent layers. A hole is part of both spaces that it connects. A space can contain several holes, one 
for each sufficiently overlapping space in an adjacent layer. Figure 6 illustrates holes in several 
spaces. 

When a net's topological route is embedded, a hole is converted into two via nodes, one in 
each space connected by the hole. In the embodiments described below, the two via nodes have an 
association (e.g., have a reference to one another) but do not have a wall (Le. t a wire segment) 
between them. Other embodiments might specify walls between related via nodes. A Steiner is a 
Steiner point (i.e., a vertex that connects to more than two walls) that is generated at the juncture 
between a pre-existing path of a net and a newly embedded path of the net. 

Figure 7 provides examples of topological routes that are formed by via nodes, Steiners, 
walls, joints, and nodes. Specifically, this figure illustrates two embedded topological paths 700 and 
702 (where a path in this instance is portion of a route). Both of these paths traverse several 
triangular faces on layer 2. Path 700 also traverses two triangular faces on layer 3, as it is a multi- 
layer route. 

Path 700 is specified by the following ordered list of particles: node 704, wall 706, joint 708, 
wall 710, joint 712, wall 714, node 716, node 718, wall 720, joint 722, wall 724, and node 726. Path 
702 is specified by the following ordered list of particles: node 728, wall 730, joint 732, wall 734, 
joint 736, wall 738, Steiner 740, wall 742, and node 744. Figure 7 also illustrates a wall 746 of 
another portion of path 702's route. 

As shown in Figure 7, the embedded routes divide the faces that they completely cross into 
several spaces. For instance, these paths divide the face formed by vertices 704, 748, and 750 into 
three spaces. Also, a joint specifies the intersection of each route with each edge. Each joint connects 
to two exits, as each joint divides a previously undivided exit into two exits. 

As mentioned above, path 702 includes a Steiner 740. This Steiner connects walls 738, 742, 
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and 746. In addition, multi-layer route 700 includes two via nodes 716 and 718. These via nodes 
enable this path to traverse layers 2 and 3. These via nodes correspond to the hole 605 illustrated in 
Figure 6. Specifically, the Q* engine converts this hole into the two via nodes 716 and 718 when it 
embeds the path 700. 

5 Each net has a representation of its connectivity and topology. Each edge has an ordered list 

of the routes that cross it. In other words, a route's topological particles that are on edges appear on 
their edge's ordered list (e.g., linked list) in the topological order that their route crosses their edge. 
For the example illustrated in Figure 7, the Q* engine specifies the edge between nodes 748 and 750 
as an ordered list of the following particles: node 750, exit 752, joint 708, exit 754, joint 736, exit 
10 756, and node 748. In this edge's ordered list, joints 708 and 736 of paths 700 and 702 appear in the 
order that their routes cross this edge. 

As mentioned above, a topological route is a route that is defined in terms of its relation to 
other layout items. The Q* engine uses two different sets of associations to specify the topology of 
each topological route. The first set is the ordered list of topological particles on the route, and the 
15 second set is the ordered list of topological particles on the edges intersected by the route. 

In the embodiments described below, certain particles (such as joints, via nodes, walls, and 
Steiners) that define a route are moveable, while other route-defining particles (such as port- 
geometry nodes) have a fixed location. The topological engine does not need to specify and store the 
location of moveable route-defining particles, as the routes are topologically defined. However, in 
20 the embodiments described below, the Q* topological engine specifies and stores locations for the 
route-defining particles in order to be able to compute wirelength accurately, 
b. Overall Flow of the Q* Topological Router 
Figure 8 illustrates a process 800 that provides the overall flow of the Q* topological engine 
in some embodiments. As shown in this figure, the process initially triangulates (at 805) the sub- 
25 region defined at 205. This triangulation will be further described below by reference to Figure 9. 
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After the triangulation, the process defines (at 810) an initial grouping of the nets in the sub-region. 
Different embodiments define the initial grouping differently. Some embodiments use a clustering 
approach that groups the nets that maximally interfere with each other. Some of these embodiments 
cluster the nets by (1) defining, for each net, a three-dimensional bounding box that contains the net's 
pins, (2) pair wise intersecting the bounding box of each net with each of the other nets, 
(3) computing the volume of intersection between the intersections, and (4) grouping the nets that 
have the greatest overlapping volume (Le., that are likely to have the most amount of overlap). In 
some embodiments, the bounding boxes are rectangular bounding boxes that are defined with respect 
to a coordinate system that is rotated by 45° with respect to the coordinate system that is aligned with 
the sides of the layout. Some embodiments might use a multi-clustering approach that allows some 
or all nets to be in more than one clustered group. 

After grouping the nets, the process selects (at 815) one of the groups identified at 810 for 
routing. Different embodiments make this selection differently. Some embodiments select the groups 
based on an initial ranking that depends on descending entropy values (i.e. t descending estimated 
number of shortest routes) for all the nets in each group. Some embodiments specify the entropy of a 
net by (1) defining a grid of horizontal and vertical lines that are spaced apart by the minimum line to 
via spacing, (2) identifying the number of horizontal grid lines m and vertical grid lines n between 
the pins of each net, (3) computing the number of different n-element sets in an m-n element set (Le. 9 
f m-n^ 



computing 



n j 



), and (4) setting the net's entropy value equal to the log base 2 of this 



computation. The process 800 iterates through 815 several times. During some of the subsequent 
iterations, some embodiments select a group at 815 based on the number of unrouted or illegally 
routed nets in the group. 

After selecting a group at 815, the process 800 calls (at 820) a solving engine to find a 
topological route for each net in the selected group. In the embodiments described below, the solving 
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engine typically defines a legal or illegal route for each net in the group. In the circumstances that it 
cannot even find an illegal route for one or more nets, the solving engine notifies the process 800, 
which then takes this failure into account at 825 and 830 described below. The operation of the 
solving engine will be further described below by reference to Figure 21. 
5 Once the solving engine returns its solution (which may be incomplete) at 820, the process 

updates (at 825) the net groupings, if necessary. When necessary, the process updates the net groups 
based on the solution returned by the solving engine. The process 800 identified (at 810) the initial 
groupings based only on an estimated measure of interference between the nets. However, each time 
the process 800 receives a solution for a group of nets from the solving engine at 820, this process 

10 can receive or derive a more accurate measure of the interference between the nets in the groups. 
Based on this measure, the process might update the net groupings at 825. For instance, if the 
solution returned by the solving engine specifies that a net A uses the same edge as a net B, then the 
process updates the net groupings so that nets A and B are more likely to be clustered. 

At 830, the process determines whether it should terminate. In some embodiments, the 

15 process terminates only if it has found a legal route for all the nets in the region, or if it has 
unsuccessfully tried a particular number of times to find a legal solution for all the nets. If the 
process determines at 830 that it should not terminate, it returns to 815 to select another group of 
nets, and then repeats 820-830 for this newly selected group of nets. The process ends when it 
decides to terminate at 830. 

20 c. Triangulation 

Different embodiments use different decomposition techniques to create a topological 
structure for embedding topological routes. The embodiments described below use triangulated 
graphs for each layer of the routed sub-region. Specifically, these embodiments use a constrained 
Delaunay triangulation ("CDT") technique. Several such techniques are disclosed in C.L. Lawson, 

25 "Transforming triangulations", Discrete Math, 3:365-372, 1972; C.L. Lawson, "Software for C 



Surface Interpolation," In J.R. Rice, editor, Math Software m . pp 161-194, Academic Press, New 
York, 1977; L. J. Guibas, D. E. Knuth, and M. Sharir, "Randomized Incremental Construction of 
Delaunay and Voronoi Diagrams", Algorithmica, 7:381-413, 1992. 

Figure 9 illustrates a triangulation process 900 that is used in some embodiment of the 
5 invention by the Q* topological routing engine. This process decomposes each layer of the sub- 
region into several triangular faces. Accordingly, this process initially selects (at 905) one of the 
layers of the sub-region. Next, the process (1) defines (at 910) a graph node at the location of each 
corner of the sub-region on that layer, and (2) defines (at 915) a graph node at the location of each 
geometry point of a port or obstacle in the selected sub-region layer. 

10 Some embodiments also define one node for each vpin in the selected sub-region layer. Other 

embodiments define a graph node for only some of the vpins. For instance, some embodiments select 
as triangulation nodes vpins that are near interior geometry nodes that are close to the boundary. Of 
the remaining vpins, these embodiments designate as nodes every nth (e.g., 5 th ) vpin around the 
boundary. The vpins that are not specified as nodes, are specified as joints in some embodiments. 

15 Also, in some embodiments, the process 900 defines (at 915) nodes on the sides of certain 
geometries to ensure that there are sufficient edges to detect congestion in the IC region. Figure 10 
illustrates the layout of Figure 4 after nodes have been defined at each sub-region corner and at each 
port or obstacle geometry point. 

Next, as shown in Figure 9, the process creates (at 925) two triangles by dividing the region 

20 along a diagonal line connecting two of the corner node vertices. The process then successively 
inserts (at 930) individual port or obstacle nodes in the created triangles to divide the triangles into 
smaller sub-triangles. Specifically, when a new node is inserted, the triangle containing that node is 
identified, and that triangle is further triangulated by connecting the newly inserted node to the 
vertices of the identified triangle. Figure 11 illustrates this triangulation technique. In this example, 

25 two triangles 1120 and 1140 are created by connecting two diagonal nodes 1105 and 1110 of the 
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sub-region layer. Next, a node 1115 is inserted in the sub-region. The triangle 1 120 that contains this 
newly inserted node is then further triangulated into three smaller triangles 1 125, 1 130, and 1 135 by 
connecting the newly inserted node to the vertices of triangle 1120. Each time the Q* topological 
engine defines (e.g., at 930) an edge between two nodes, it defines an exit between the two nodes. 
Also, each time the Q* topological engine defines a triangle (e.g., at 930), it defines the space within 
the triangle. 

After defining a set of triangles at 930, the process performs (at 935) an edge-flipping 
operation to maximize the minimal angle of each triangle, i.e., to make the triangulation Delaunay. 
This operation is done in order to improve the likelihood that the topological routes produced by the 
Q* topological engine can be converted into specific geometric routes. To have an absolute 
guarantee that the generated topological routes can be geometrized, a visibility graph needs to be 
constructed to analyze every edge between any two graph nodes that have unobstructed views of 
each other in the relevant routing metric to ensure that each such edge is not overcongested. 
However, such an approach would be computationally expensive. Hence, instead of examining the 
edge between each such pair of graph nodes, some embodiments use an edge-flipping Delaunay 
triangulation as a good approximation of the visibility graph. 

Figures 12 and 13 illustrate why maximizing the minimal angle of each triangle improves the 
likelihood that the generated topological routes can be geometrized. In these figures, nodes 1205 and 
1210 do not have an edge between them. Hence, the Q* topological engine cannot measure the 
congestion of the straight-line path between these two nodes. It can, however, measure the 
congestion on edges 1215, 1220, 1225, 1230, and 1235. In Figure 12, the triangles are equilateral 
triangles, and therefore have the largest minimal angles. As illustrated in Figure 12, it is relatively 
unlikely that a set of topological routes exist that could overfill the capacity of the straight-line path 
between nodes 1205 and 1210 without overcongesting the capacity of edge 1235. However, as 
illustrated in Figure 13, when the adjoining triangles have small minimal angles, it is quite possible 
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that a set of paths overcongest the straight-line path between nodes 1205 and 1210 without 
overcongesting the capacity of the adjoining edge 1235. 

Some embodiments perform an edge-flipping operation by identifying, for each triangle, a 
circle that encompasses all of the triangle's vertices. If that circle encompasses the vertex of another 

5 triangle as well, and if the two triangles do not jointly form a non-convex polygon, then the common 
edge between the two triangles is flipped. Flipping an edge between two triangles means deleting the 
existing common edge between the two triangles and defining a new common edge between the two 
vertices of the triangles that were not previously connected. The edge flipping operation typically 
results in a new pair of triangles that has larger minimal angles than the original pair of triangles. 

10 This operation defines necessary exits and spaces for the new triangles, and discards the original pair 
of triangles and their associate spaces. When a pair of abutting triangles form a non-convex structure, 
the common edge between them is not flipped. . 

Figures 14 and 15 illustrate one manner for performing the edge-flipping operation. Figure 
14 illustrates two triangles 1405 and 1410. Circle 1415 is the identified circle that encompasses all of 

15 the vertices of triangle 1405. This circle also includes the vertex 1420 of triangle 1410. As the two 
triangles do not form a non-convex polygon, the common edge 1425 between these two triangles is 
flipped. Figure 15 illustrates the newly defined pair of triangles 1430 and 1435 and the edge 1440 
between them. The edgel425 between the triangles 1405 and 1410 would not have been flipped if 
the circle 1415 did not contain the vertex 1420 of triangle 1410. 

20 After performing the edge-flipping operation at 935, the process performs (at 940) a 

constraining operation to ensure that a triangle edge exists (1) between each pair of successive points 
of an obstacle or port, and (2) at the sub-region boundary edges. Figure 16 provides a pictorial 
example of this operation. This figure illustrates two successive pairs of obstacle or port points 1605 
and 1610 that do not share an edge after the edge-flipping operation of 935. Figure 16 illustrates the 

25 constraining operation flipping three different edges 1615, 1620, and 1625 until an edge 1630 exists 
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between the two node points 1605 and 1610. Once this edge exists, it is identified as a constrained 
edge that should not be removed. Accordingly, this operation is referred to as a constraining 
operation, as it defines edges at sub-region boundaries and between successive pairs of geometry 
points, and then marks the edges as constrained. 

Next, the process performs (at 945) a follow-up edge-flipping operation without flipping the 
edges constrained at 940. This follow-up operation is to try to maximize the minimal angle of the 
triangles as much as possible while respecting the constrained edges. The process then determines (at 
950) whether it has examined all the sub-region's layers. If not, the process transitions back to 905 to 
select another layer and to repeat 910 through 945 to triangulate the newly selected layer. 

Otherwise, the process identifies (at 955) holes in the spaces of the triangles that it defined. 
Figure 17 illustrates a hole-identification process 1700 that the Q* topological routing engine can 
use to try to specify holes for the spaces of a triangle (referred to as the examined triangle). After 
triangulation, the topological engine performs (at 955) the process 1700 for each triangle that it has 
defined. This engine also uses this process at other stages in the routing process. In general, this 
engine uses all or some of this process each time it creates a new space (e.g., when it creates new 
triangles, or when modifies previously defined triangles). 

The process 1700 initially identifies (at 1705) all triangular faces that overlap the examined 
triangular face in the layer above, and if applicable, in the layer below. The overlapping faces can be 
identified based on the x- and y- coordinates of the nodes of the faces. In some embodiments, the 
topological engine stores for each face a list of all faces that overlap it in the layers above and, if 
applicable, below. 

The process then identifies (at 1710) the overlapping faces that are potential via destinations 
from the examined face. Specifically, for each triangle, the Q* topological engine in some 
embodiments stores a polygon that approximates the legal area where the center of vias can be in the 
face in the absence of any route crossing the face. The engine identifies the polygon for each face by 
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identifying the closest point within the face that a via can be to each obstacle on the layer, accounting 
for the via size and/or shape, the required spacing to obstacles, etc. If the result is disjoint or non- 
convex, the engine approximates by using the convex hull of the legal locations. 

The process 1700 uses the face polygons to identify quickly the overlapping faces that might 
5 be potential via destinations. In particular, at 1710, the process 1700 (1) intersects the examined 
face's polygon with the polygon of each overlapping face that it identified at 1705, and (2) identifies 
each overlapping face as a potential via destination if the intersection polygon computed for it has 
sufficient space for a via. The process 1700 then determines (at 1715) whether it identified at 1710 
any overlapping face as a potential via destination. If not, the process ends. Otherwise, the process 

10 selects (at 1720) a potential via-destination face that it identified at 1710. 

For each space within each face, the Q* topological routing engine stores a polygon that 
specifies the legal area where the center of vias can be in the space. The engine identifies the polygon 
for each space by identifying the closest point within the space that a via can be to each obstacle on 
the layer, accounting for the via size and/or shape, the required spacing to obstacles, the required 

15 spacing to other routes crossing the face, etc. If the result is disjoint or non-convex, the engine 
approximates by using the convex hull of the legal locations. A space's polygon might be identical to 
its face's polygon (which, in some embodiments, is the case after the initial triangulation) or might be 
different than its face's polygon (which, in some embodiments, is the case when topological routes 
cross the face). 

20 The process 1700 specifies a hole for each pair of overlapping spaces when the intersection 

of the polygons of these spaces is at least a certain size (e.g., it is larger than a particular threshold 
value). Specifically, after selecting (at 1720) a potential via-destination face, the process 1700 
identifies (at 1725) an intersection polygon for each combination of space pairs in the examined and 
selected faces. When the intersection polygon for a space pair is at least a certain size, the process 

25 then identifies (at 1730) a hole for the space pair, if one such hole does not already exist. The process 
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then determines (at 1735) whether it has examined all the potential via-destination faces that it 
identified at 1710. If not, the process transitions back to 1720 to select another such face. Otherwise, 
the process ends. 

The end result of the triangulation operations 900 is a set of triangulated sub-region layers, 
5 which can be used to embed and remove topological routes in the sub-region. Figure 18 illustrates 
an example of how the layout of Figure 4 might look after triangulation. 

In some embodiments, the triangulation operation defines the capacity of each edge in the 
triangulated region. The above-described triangulation operation defines each edge's capacity 
whenever it creates the edge at 925 through 945. Figure 19A illustrates a process 1900 that the 
10 topological engine can call each time it wants to identify the capacity of an edge in the triangulated 
sub-region. This process is designed to support multiple wiring models. 

As shown in Figure 19A, the process 1900 determines (at 1905) whether either node of the 
edge belongs to a vpin. If so, the process identifies (at 1910) the edge as the capacity vector. If not, 
the process identifies (at 1915) the capacity vector as the vector that traverses the shortest distance 
15 between the sides of the geometry abutting one of the edge nodes and the sides of the geometry 
abutting the other edge node. 

After 1910 or 1915, the process 1900 defines the edge capacity as the length of the largest 
projection of the capacity vector onto one of the legal routing directions. The projection P of the 
capacity vector C onto a legal routing direction D is given by 
20 P = C * cos a, 

where a is the angle between the capacity vector C and the legal routing direction D. Accordingly, 
the edge capacity is the magnitude of the projection vector. The largest projection of the capacity 
vector can be identified (at 1920) in a variety of ways. Some embodiments compute the magnitude of 
the capacity vector's projection onto all the legal routing directions and then select the largest. Others 
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identify the routing direction that makes the smallest angle with the capacity vector, define this 
routing direction as the direction of projection, and then compute the projection of this capacity 
vector onto the identified routing direction. 

Other embodiments might compute the edge capacities slightly differently from process 
1900. For instance, some embodiments might define each edge (including an edge that does not 
connect to vpins) to be its own capacity vector. Some of these embodiments then specify each edge's 
capacity as the edge's largest projection onto one of the legal routing directions. 

Figure 19B illustrates the different ways discussed above for defining the capacity of an 
edge. Figure 19B illustrates an edge 1925 between two nodes 1930 and 1935 on an octilinear layer 
like the one shown in Figure 1. In this example, vector 1965 is the vector that traverses the shortest 
distance between sides 1945 and 1950 abutting node 1930 and sides 1955 and 1960 abutting node 
1935. Some embodiments would select the edge 1925 as the capacity vector. For these embodiments, 
the projection 1975 is the largest projection of this capacity vector onto one of the legal routing 
directions, which in this case is the +90° direction. On the other hand, some embodiments would 
select the vector 1965 as the capacity vector. For these embodiments, the projection 1970 is the 
largest projection of the capacity vector 1965 onto a legal routing direction, which in this case is the 
+45° direction. Defining an edge's capacity based on its capacity vector's largest projection onto a 
legal routing direction provided by the wiring model enables the Q* topological engine to produce 
topological routes that are customized for different wiring models. 

The embodiments described above identify an edge's capacity as the largest projection of the 

edge's capacity vector onto one of the legal routing directions. Other embodiments, however, might 

use other techniques to compute an edge's capacity from its capacity vector. For instance, some 

embodiments that use the octilinear wiring model of Figure li (l)for each layer, identify eight 

sectors about the eight routing directions, (2) identify the sector of the edge's layer that contains the 

edge's capacity vector, and (3) project the capacity vector onto the routing direction associated with 
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the identified sector. 

Figure 20A illustrates eight sectors that some embodiments identify for a particular layer 
about the eight routing directions of Figure 1. These eight sectors are defined between eight 
constraining angles: 22.5°, 67.5°, 112.5°, 157.5°, 202.5°, 247.5°, 292.5°, and 337.5°. Some 
embodiments assume that these angles are the angles of edges that constrain the embedding of 
geometric routes about unrelated geometries in the layout. Some of these embodiments assume that 
net routes have equal widths and spacings in the Manhattan and diagonal directions on the particular 
layer. Computing the edge capacities by using the eight constraining angles and eight sectors 
illustrated in Figure 20A is equivalent to defining the edge capacities as the largest projection of the 
capacity vectors onto the legal routing directions. 

Other embodiments, however, define the eight constraining angles of a particular layer 
differently. Therefore, they compute some edge capacities that are different than the largest 
projection of the edge capacity vectors onto the legal routing directions. These embodiments allow 
the widths and spacings of nets to be different in the Manhattan and diagonal directions on the same 
layer. These embodiments define the constraining angles on a layer based on the worst combination 
of net spacings and widths in the Manhattan and diagonal directions on the layer. Figure 20B 
illustrates eight such constraining angles for a particular layer. In this figure, the angle a can be 
computed based on the equation (1) below. 

S • +W - 

a = arc tan^ =£. (1) 

To understand this equation, the following variables need to be defined: 

• i is a number from 1 to n representing one of n nets; 

• k is a number from 1 to p representing one of p spacing constraints between nets 
and between nets and obstacles; 

• SJJ 1 is a k lh spacing between a net and another net or an unrelated geometry in a 
Manhattan direction of the particular layer; 

• is a width of the i" 1 net in a Manhattan direction of the particular layer; 

• Si is the k 1 * 1 spacing constraint between a net and another net or an unrelated 
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geometry in a diagonal direction of the particular layer, 

• Wi d is the width of the i* net in the diagonal direction of the particular layer. 

• The values S k and W; y are computed as follows: 

S ky =ceil(s^*V2-S^);and (2) 
5 W uy =cti\(w?*j2-Wr), (3) 

where the "ceil" operation in the above equations signifies rounding up to the next manufacturing 
grid, as the spacings and widths are expressed in terms of manufacturing grid units and the routes are 
defined with respect to the manufacturing grid. Given the above definitions, the variables S^ y and 

10 S™ n are the related k lh spacing constraints S ky and that produce the smallest ratio 

Similarly, the variables and W* n are the related i lh width constraints Wj y and W, m that 

W iy 

produce the smallest ratio — jj- . 

Wj 

Figure 20C provides an example that illustrates some of the variables mentioned above, and 
describes the derivation of the equation (1). This figure illustrates two nets 2005 and 2010 bending 
15 around an obstacle 2015 of a particular layer. In this figure, 

• is the 80 lh spacing constraint specifying the spacing between the net 2005 and 
the obstacle 2015 in the Manhattan directions of the particular layer, 

• S£> the 80 lh spacing constraint specifying the spacing between the net 2005 and 
the obstacle 2015 in the diagonal directions of the particular layer, 

20 • W 5 m is the width of the net 2005 in the Manhattan directions of the particular 

layer, 

• W 5 D is the width of the net 2005 in the diagonal directions of the particular layer, 

• S"^ is the 122 nd spacing constraint specifying the spacing between the nets 2005 
and 2010 in the Manhattan directions of the particular layer, and 

25 • Sf 22 is 122 nd spacing constraint specifying the spacing between the nets 2005 and 

2010 in the Manhattan directions of the particular layer, 

• W,o is the width of the net 2010 in the Manhattan directions of the particular 
layer, 

• W,q is the width of the net 2010 in the diagonal directions of the particular layer, 
30 • S"J 6 is the 156 lh spacing constraint specifying the spacing between the net 2010 

and an obstacle 2050 in the Manhattan directions of the particular layer, and 

• S, D 56 is 156 lh spacing constraint specifying the spacing between the net 2010 and 
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the obstacle 2050 in the diagonal directions of the particular layer. 
Figure 20C illustrates S B>y , W 5y , S 122 y , W 10y , S IM y . In this figure, the points that define the route 

bends in this case, the transitions from the -45° direction to the vertical direction) are specified 
5 on the manufacturing grid. As shown Figure 20C, 

• Sgoy is the difference in the y-coordinates of the vertex 2020 and the bend- 
defining point 2025 of the net 2005, 

• W s is the difference in the y-coordinates of the bend-defining points 2025 and 

2030 of the net 2005, 

10 • S 122 is the difference in the y-coordinates of the bend-defining point 2030 of the 

net 2005 and the bend-defining point 2035 of the net 2010, 

• W 10 y is the difference in the y-coordinates of the bend-defining points 2035 and 

2040 of the net 2010, and 

• S 156 y is the difference in the y-coordinates of the bend-defining point 2040 of the 

15 net 20 10 and the vertex 2045 of the obstacle 2050. 

The five equations below specify how S 122y , S my , W Sy , W 10y , and S 156 y are computed in order to 

define the bend-defining points on the manufacturing grid. 

S a , y =ceilfe*>^-SS); 
20 W J>y = ceil(w 5 d * V2 - W 5 m ) ; 

S.aj-calte^VS-SZ,); 
W l0>y =ceil(w, d 0 *V2-W™); 

S l56 , y =ceil(s? 56 *V2-SrJ. 
In the example illustrated in Figure 20C, the angle a 0 of a line connecting the obstacle 



25 vertex 2020 and the first bend point 2025 is arctan 



°80.y 



The angle a, of a line connecting the 



obstacle vertex 2020 and the bend point 2035 is arctan 
connecting the obstacle vertices 2020 and 2045 is: 
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^O.y _ S.y J l22.y 



The angle oc 2 of a line 



a, = arctan 



m . cm 



ss+w 3 n, +s™ 22 + w™+sr 56 

This angle a, will always be greater or equal to the angle oc p : 



a p = arctan 



{ ps: n +(p-i)wi j* 

where p-1 is the number of nets crossing a constraining edge. Equation (1) is obtained from the 
5 above equation, by realizing that for all p the angle ct p will be greater or equal to the angle provided 

by equation (1). The equation (1) reduces to: 

a = arctan V W t y 

s m + wi 

when the particular layer has only one spacing constraint S m in the Manhattan directions and one 
spacing constraint S d in the diagonal directions. If the particular layer has only one width W m in the 
10 Manhattan directions and one width W d in the diagonal directions for all nets, the equation (1) 
reduces further to: 

S y + W y 

a = arctan 

S m +W m 

After using equation (1) to compute the angle a for a particular layer, some embodiments 
specify the eight constraining angles for the particular layer. These eight constraining angles are 

15 illustrated in Figure 20B. These eight angles specify eight sectors, where each sector is associated 
with a particular direction. These embodiments identify the capacity of each edge that is within the 
particular layer by first identifying the sector that contains the edge's capacity vector, and then project 
each edge's capacity vector onto the routing direction associated with the identified sector. For 
instance, if the angle of the capacity vector is between a and 90 -a, the capacity vector falls within 

20 the sector between constraining angles a and 90- a. As shown in Figure 20B, this sector's 
associated routing direction is the 45° direction. Accordingly, the capacity vector will be projected 
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onto the 45° direction. 

d. Solving Engine 
(1) Overall Flow 

As mentioned above, the Q* topological routing engine calls (at 820) its solving engine to 
identify topological routes for a group of nets. Figure 21 illustrates a process 2100 that the solving 
engine performs in some embodiments each time that the process 800 calls it at 820 and specifies a 
selected group of nets to route within a given sub-region. For each group of nets that it receives, the 
process 2100 generates several solutions based on different net ordere, and then returns the best 
solution that it generates for the group. This exploration of the solution space is further enhanced by 
the modifications of the net groupings at 825. 

As shown in Figure 21, the process 2100 starts by setting (at 2105) a Best_Cost variable to a 
large number. This process uses this variable to identify the best solution for the nets that it 
examines. It also uses this variable to determine whether it should terminate generating routes based 
on a particular net order. 

The process then identifies (at 2110) an ordering of the nets. As further described below, the 
process 2100 passes through 2110 several times. The process 2100 selects a different net order each 
time it returns to 2110. Different embodiments iterate through the different net orders differently. 
Some embodiments examine all ordering permutations, while other embodiments examine fewer 
permutations. Also, at 2110, some embodiments randomly select new net orderings, while other 
embodiments systematically examine different ordering permutations. 

Some embodiments systematically examine different net orderings by using a heuristic to 
identify a net order the first time the process 2100 passes through 2110, and then specifying a 
different net order during each subsequent iteration of the process 2100 through 2110. For instance, 
in some embodiments, the process 2100 specifies an initial net order based on descending net 
entropies. Entropy provides a measure of difficulty for routing the nets. A high entropy net has many 
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equivalently good options (z'.e., has many shortest-cost routes), while a low entropy net has only a 
few good routes (Le. t has few shortest-cost routes). Some embodiments measure a net's entropy 
according to the approach described above for 815. 

After specifying the initial net order based on a heuristic, the process 2100 in some 
embodiments specifies the subsequent net orders at 2110 by (1) identifying the nets that ended up 
with routes that were longer than anticipated in the previous iteration, and (2) specifying an ordered 
list with the identified nets higher on the list. Other embodiments use a least-discrepancy-search 
("LDS") technique. Under this technique, the process 2100 still specifies the initial net order based 
on a heuristic the first time that it passes through 2110. During each subsequent iteration through 
2110, the process selects a net order that has the fewest differences from the previously selected 
orderings. For instance, after the initial net order, the next series of iterations through 2110 select 
orderings that are identical to the initial ordering except that the positions of two of the nets are 
transposed. Similarly, the subsequent series of iterations through 2110 would explore orderings that 
are identical to the initial one except that the number of transposed net positions would be three, 
four, five, etc. Branch-and-bound techniques can be used to terminate iterating through a series of 
different orderings {e.g., terminate examining different ordering permutations that have the position 
of four nets transposed with respect to the initial ordering). Some embodiments that explore fewer 
than all possible permutations of net orders limit the number of permutations that they explore based 
on the number of nets in the received group of nets. 

After specifying an ordered list of nets at 21 10, the process rips out (at 21 15) any topological 
route previously specified in the sub-region for any nets in the received group of nets. At 2115, the 
process also initializes a variable Solution_Cost to zero. Next, from the ordered list of nets specified 
at 21 10, the process selects (at 2120) a net to route. It then calls (at 2125) a route-generation process 
to identify a route for the selected net. This route-generation process will be described further below 
by reference to Figure 22. 
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Next, the process determines (at 2130) whether the route generation engine was able to 
identify a route for the selected net. If not, the process transitions to 2155, which will be further 
described below. Otherwise, the process increments (at 2135) the Solution_Cost by the cost of the 
route identified at 2130. Next, the process determines whether it should stop (at 2140) identifying 
5 routes for the current net order specified at 2110. In some embodiments, the process terminates (at 
2140) finding routes for the current net order when the cost of the nets routed for the current net 
order (Le., the current Solution_Cost) plus the estimated cost of the unrouted nets in the current net 
order exceeds the Best_Cost encountered thus far. 

In some embodiments, the estimated cost of the unrouted nets equals the sum of the estimated 
10 cost for each unrouted nets. In some embodiments, each unrouted net's estimated cost equals the 
length of the net's minimum octilinear spanning tree or Steiner tree. Several manners for constructing 
octilinear spanning and Steiner trees for a net are described in United States Patent Application 
entitled "Method and Apparatus for Considering Diagonal Wiring in Placement" and having serial 
number 09/731,891. 

15 As mentioned above, the Best_Cost is initially set to a large number at 2105 for the group of 

nets. This cost is set to a lower number when the process 2100 finds a routing solution for the initial 
net order, assuming that at least one topological route can be defined for each net in the order 
irrespective of the route's cost. Each time the process 2100 finds a new solution for all the nets based 
on a different ordering, it sets the Best_Cost to the new solution's cost if the new solution has a lower 

20 cost than the previously defined BesLCost. 

If the process determines (at 2140) that it should stop identifying routes for the current net 
order, it transitions to 2155, which is further described below. Terminating a search for a particular 
net order when it is apparent that the search will not provide a better solution than the best solution 
already encountered speeds up the operation of the process 2100 so that it can examine more net 

25 orderings. 
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If the process determines (at 2140) that it should not terminate its search for the current net 
order, it determines (at 2145) whether it has examined all the nets on the current ordered list 
specified at 21 10. If not, the process returns to 2120 to select another net and to repeat its operations 
for this newly selected net. Otherwise, the process determines (at 2150) whether the Solution_Cost 
5 for the current net order is better than the Best_Cost recorded thus far. If so, the process sets (at 
2150) the Best_Cost to the Solution_Cost and records the solution (/.<?., the set of topological routes) 
for the current order of the nets as the best solution (£.«., as the best set of topological routes). 

The process then determines (at 2155) whether it has sufficiently examined different net 
orders. This determination depends on the technique that the process is using to search through the 
10 space of different net orders. These different techniques were described above. If the process 
determines (at 2155) that it should examine other net orders, it returns to 2110 to select another net 
order and then repeats its subsequent operations to identify a solution based on this order. Otherwise, 
the process determines (at 2160) whether it found a solution that included a route for all the nets in 
the received group of nets. If not, the process returns (at 2165) its failure to find a solution and the 
15 best partial solution that it encountered. 

Otherwise, the process determines (at 2170) whether the solution embedded for the last net 
order was the best solution it examined. If so, the process transitions to 2180, which is further 
described below. Otherwise, it rips up (at 2175) the last embedded solution and embeds the best 
solution that it identified by examining the various net orderings. This embedding operation is 
20 similar to an embedding operation performed by the route generation process, which will be further 
described below. At 2180, the process returns the embedded solution and then ends. 
(2) Route Generation Process 
The solving engine directs (at 2125) a route generation process to construct the lowest-cost 
route for a particular net in a particular region. Figure 22 illustrates a route generation process 2200 
25 of some embodiments. The process 2200 starts (at 2205) by initializing a variable, Route_Cost, to 



zero. The process uses this variable to specify the cost of the route that it tries to construct for the 
net. 

Next, at 2210, the process orders the pins of the net. Some embodiments order the pins based 
on their estimated distance (in the appropriate wiring model metric) to a chosen source pin. The 
process then specifies (at 2215) one pin of the net as the target of the path search and the remaining 
pins as sources of the search. At 2215, the process then specifies the nodes, exits, and/or holes (if 
any) of the target pin as the target particles, and specifies the nodes, exits, and holes (if any) of the 
selected source pins as the target particles. Some embodiments require the routed particles for each 
pin to be from the same port of that pin. 

After specifying the source and target particles for a path search, the process directs (at 2220) 
the Q* path search process to identify and embed the lowest-cost path between the specified source 
and target sets. If the path search process embeds the lowest-cost path, the path search process 
increments the Route.Cost by the cost of the embedded path. The path search process is further 
described below by reference to Figure 35. 

At 2225, the route-generation process determines whether the path search process was able to 
identify and embed a path between the specified source and target sets. If not, the process 2200 has 
failed to find a route for the net. Accordingly, it returns (at 2230) a notification specifying its failure 
and then ends. 

On the other hand, if the process determines (at 2225) that the Q* path search process 
identified and embedded a path, it determines (at 2235) whether it has routed all the pins of the net. 
If so, the process 2200 notifies (at 2240) the solving engine that it has embedded a route for the net 
and provides this route's associated Route_Cost. 

If the process 2200 determines (at 2235) that it has not routed all the pins, the process 
specifies (at 2245) new source and target particles for another path search. In some embodiments, the 
process specifies the nodes, exits, and holes of the net's pins that are currently unrouted as source 
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particles, and specifies the nodes, holes, exits, and/or walls of the net's previously routed pins and 
paths as target particles. One of ordinary skill will realize that other embodiments might specify 
source and target sets differently at 2245. For instance, some embodiments might specify the 
unrouted topological particles as the set of targets while specifying the routed particles and paths as 
the set of sources. 

After specifying the source and target sets at 2245, the process 2200 returns to 2220 to direct 
the Q* path search process to identify and embed the lowest-cost path between the specified source 
and target sets. The operation of the process 2200 from 2220 was described above. 
(3) Q* Path Search Process 

The Q* path search process is a search process that can identify the lowest-cost path between 
source and target states in a graph with non-zero dimensional states. Such a path often (although not 
always) traverses through intermediate states. In a graph with non-zero dimensional states, a zero- 
dimensional state is a point, a one-dimensional state is a line, a two-dimensional state is a surface, a 
three-dimensional state is a volume, etc. 

In the embodiments described below, the tessellated region, which was described above, is a 
multi-layer graph that the Q* path search explores. In the tessellated region, the source, intermediate, 
and target states are zero-, one-, or two- dimensional topological particles, such as nodes, lines, and 
surfaces. In other embodiments, the source, target, and intermediate states might be higher 
dimensional states (e.g., three-dimensional volumes, etc.). 

In the embodiments described below, the Q* search process performs a best-first search 
operation. The following sections provide (i)an introduction to best first searches, (ii)a brief 
introduction to the Q* search, (iii) an introduction to the functions utilized by the Q* search in some 
embodiments, (iv) the overall flow of the Q* search in some embodiments, and finally (v) an 
example of a Q* search. 
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(i) Introduction to Best First Searches: A* 

A best-first search is an iterative search that at each iteration tries to extend the partial 
solution with the best estimated cost. Best-first searches have often been used to identify the shortest 
path between source and target points in a multi-point search space. One such best-first search is the 
A* search. To identify a path between source and target points, the A* search starts one or more 
paths from the source and/or target points. It then iteratively identifies one or more path expansions 
about the lowest cost path, until it identifies a path that connects the source and target points. The 
typical cost of a path expansion in an A* search is an F cost, which is the cost of the path leading up 
to the path expansion plus an estimated cost of reaching a target point from the path expansion. 

Figures 23-28 provide an example of an A* path search that uses such an F cost. In this 
example, the search process has to find the shortest path between a source point 2305 and a target 
point 2310 in a region 2320. The source and target points are part of a multi-point grid 2315 that is 
imposed over the region. 

As shown in Figure 24, the search process initially identifies four path expansions from the 
source point 2305 to four points 2402-2408 that neighbor the source point in the Manhattan 
directions. In Figures 24-28, the search process represents each path expansion by using a path 
identifier, called a "drop." More specifically, the search process represents each path expansion from 
one grid point (a start point) to another point (a destination point) by (1) specifying a drop, 
(2) associating the drop with the expansion's destination point, and (3) defining the specified drop's 
previous drop to be the drop of the expansion's start point. Drops allow the search to keep track of 
the paths that it explores. 

The search process specifies four drops 2410-2416 for the expansions to the four points 
2402-2408, as illustrated in Figure 24. It also specifies a drop 2418 for the source point 2305. The 
source point drop 2418 is the previous drop of drops 2410-2416. The source point drop's previous 
drop is null, as it is the first drop in the path search. 
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For each drop 2410-2416, the search process computes an F cost based on the following 
formula: 

F =G+ H. 

where (1) G specifies the cost of a path from the source point to the drop's grid point through the 
5 sequence of expansions that led to the drop, and (2) H specifies the estimated lower-bound cost 
from the drop's grid point to the target point. When computed in this manner, the F cost of a drop is 
the estimated cost of the cheapest path that starts at the source point, traverses through the sequence 
of expansions that led to the drop, and traverses from the drop to the target point. 

To simplify the description of the example illustrated in Figures 23-28, the distance between 
10 each pair of horizontally or vertically adjacent grid points is 1. Accordingly, in Figure 24, the G cost 
of each drop 2410-2416 is 1, as the grid point of each of these drops is one grid unit away from the 
source point. In Figures 24-28, a drop's H cost is computed as the Manhattan distance between the 
drop's point and the target point. Hence, the H cost of drops 2410, 2412, 2414, and 2416 are 
respectively one, three, three, and three, as these distances are respectively the Manhattan distances 
15 of the points of these drops from the target point. 

After costing these drops, the search process stores the drops 2410-2416 in a priority queue 

that is sorted based on their F costs. It then retrieves the drop with the lowest F cost from the 
priority queue. This drop is drop 2410. Since this drop's corresponding point (Le. 9 point 2402) is not 
the target point, the search process then identifies a path expansion from the retrieved drop's point 
20 2402 to point 2505, as shown in Figure 25. As shown in this figure, this expansion is the only viable 
expansion from the retrieved drop's point 2402 as the search has previously reached all other 
unblocked neighboring grid points (i.e., grid points that are not blocked by an obstacle 2515) through 
less expensive paths. The search process specifies a drop 2510 for this expansion, and computes this 

drop's G, H , and F costs, which are respectively 2, 2, and 4. It then stores this specified drop in the 
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priority queue. 

After storing drop 2510 in the priority queue, the search process might retrieve either drop 
2510, drop 2412, or drop 2414 from the priority queue, as each of these drops has an F cost of 4. 
However, if the search process retrieved drop 2510, it will not expand from this drop to its 
neighboring points that are not blocked by obstacle 2515 since all these neighboring points were 
previously reached less expensively. Also, if the search process retrieves drop 2414, it will identify 
drops that will be more expensive than drop 2412. 

When the search process retrieves drop 2412 from the priority queue, it checks whether this 
drop's point is the target point. When it discovers that it is not, the process (1) identifies expansions 
to three neighboring points 2602-2606 about this drop's point, as shown in Figure 26, (2) specifies 
three drops 2608-2612 for the three identified expansions, as shown in Figure 26, (3) computes each 
specified drop's G, H , and F costs, (4) defines each specified drop's previous drop (which in this 
case is drop 2412), and (5) stores each newly specified drop in the priority queue based on its F cost. 

Next, as illustrated respectively in Figures 27 and 28, the search process performs these six 
operations first for drop 2608 and then for drop 2710, since these two drops are the ones with the 
lowest F costs during the next two iterations of the search process. As illustrated in Figure 27, the 
drop 2710 is specified for an expansion about drop 2608's point 2602. 

As shown in Figure 28, one of the expansions about drop 2710 reaches the point 2310. The 
search process creates, costs, and stores a drop 2815 for this expansion. It then retrieves this drop in 
its next iteration, and then realizes that this drop's point is the target point. Accordingly, at this 
juncture, it terminates its path search operation. It then commences a path-embedding, back-trace 
operation that uses the previous-drop references of the drops 2815, 2710, 2608, 2412, and 2418 to 
identify the sequence of drops that reached the target point 2310 from the source point 2305. This 
operation embeds a path along the grid points associated with the identified sequence of drops. 
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Figures 29 and 30 illustrate the back-trace operation and the resulting embedded path, 
(ii) Brief Introduction to Q* 

The A* search is not suitable for finding the lowest-cost path in a graph with non-zero 
dimensional states. This is because the A* search computes a single cost value for the expansion to 
any state in the graph, while the actual cost can vary across a non-zero dimensional state. 
Accordingly, some embodiments use the Q* search, which can identify the lowest-cost path in a 
graph with non-zero dimensional states. In the embodiments described below, the source, 
intermediate, and target states in the graph are zero-, one-, or two- dimensional topological particles, 
such as nodes, lines, and surfaces. 

In the embodiments described below, the Q* search process has two phases: (1) a path 
exploration phase, during which the process identifies a path between the specified source and target 
particles, and (2) a path-embedding phase, during which the process embeds the identified path. 
During the path exploration phase, the Q* search starts one or more paths from the source and/or 
target particles. It then iteratively identifies one or more expansions about the lowest cost path, until 
it identifies a path that connects the source and target states. Each identified expansion about a path 
is from a "current particle" (also called "start particle") reached by the path to a "destination particle" 
that neighbors the current particle. 

Like the A* search example described above, the embodiments described below cost each 
expansion based on an F cost that can be expressed as: 

F = G + H . (4) 
In the A* search example, the F , G, and H costs always provide a single value for a state in a 
search space. However, in a Q* search, these costs can be functions that are defined over the entire 
destination particle of an expansion or one or more portions of this destination particle. 

Accordingly, in equation (4), the F cost is a function that expresses the estimated cost of a 
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path that traverses from a source through the expansion's destination particle to a target. The G cost 
is a function that expresses the cost of the path that has reached the expansion's destination particle. 
The H cost is a function that expresses an estimated cost of a path from the expansion's destination 
particle to the target set. In the embodiments described below, the H cost function expresses the 
lower-bound estimate of the shortest path from the expansion's destination particle to the target set. 
Accordingly, in these embodiments, the F cost function expresses the estimated cost of a lowest- 
cost path from the source through the expansion's destination particle to a target. Also, in these 
embodiments, the G function and hence the F function account for several different types of path 
costs, such as a distance cost, route-piercing cost, via cost, wire-shoving costs, etc. Hence, the F 
function of equation (4) allows the embodiments described below to identify the lowest-cost path 

between the source and target sets. Other embodiments, however, might utilizes a different F 
function, as further described below. 

(iii) PLF's and Filter Function 

In the embodiments described below, the G, H , and F functions are convex piecewise linear 
functions ("PLF"), although in other embodiments they might be other types of functions. In the 
embodiments described below, each PLFs domain is either a point, a line, or a surface in the 
tessellated region. In the discussion below, PLFs that are defined over a point are called point PLFs, 
PLFs that are defined across lines are called line PLFs, and PLFs that are defined across surfaces 
are called surface PLFs. 

A point PLF maps the point over which it is defined to a single value. A line PLF is a PLF 
with a domain that is a line. Such a PLF can be specified by a sequence of vertices, called knots, that 
represent the endpoints of its line segments. Figure 31 illustrates an example of a line PLF 3110 that 
has four line segments. This PLF specifies a PLF-value for each point P that is offset by an amount Q 
along an edge 3105 (Le. t for each real number Q that defines a point P on the edge at Q • u, where u 
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is a unit vector of L). Each knot of a line PLF can be specified in terms of an offset value Q and a 
PLF-value V. Also, the knots of a line PLF can be sorted in an order based on their offset values. 
Accordingly, the five knots K1-K5 that are at the endpoints of PLF 3110's four line segments can 
represent this PLF of Figure 31. These five knots can be represented by the following ordered list of 
5 offset and PLF-value pairs: (Ql, VI), (Q2, V2), (Q3, V3), (Q4, V4), (Q5, V5). 

A surface PLF is a set of one or more planar surfaces, called facets. Figure 32 illustrates an 
example of a surface PLF. This PLF 3210 has four facets (F1-F4), each of which has a different 
slope. This PLF is defined across a surface 3205. For each x-y value in the surface 3205, the surface 
PLF 3210 provides a PLF-value (V). 

10 Each surface PLF can be represented by a set of knots, a set of edges between the knots, and a 

set of facets defined by the edges. Using this approach, the surface PLF 3210 of Figure 32 can be 
represented by a list of knots K1-K10, a list of edges E1-E13, and a list of facets F1-F4. Some 
embodiments represent (l)a surface-PLF knot with an x,y coordinate, a PLF-value at that 
coordinate, and a list of edges that are incident upon the knot, (2) a surface-PLF edge by a pair of 

15 references to two knots that the edge connects, and (3) a surface-PLF facet by a list of edges, a 
normal vector (e.g., x, y, z coordinate values), and a z-intercept, where z is the axis in which the 
PLF-values are defined. For instance, in the example in Figure 32, the list of edges for the knot Kl 
specifies edges El and E2, the list of knots for edge E4 specifies knots K3 and K4, and the list of 
edges for facet Fl identifies edges E1-E4. For facet Fl, a normal vector Nl and a z-intercept are also 

20 specified. 

In the embodiments described below, the Q* search maintains another function, called the 
filter function, for each topological particle that can serve as a source, intermediate, or target state in 
the tessellated region. In the embodiments described below, the filter function is a PLF. A particle's 
filter function expresses the lowest path cost that the search process has been able to identify from 
25 the source set to each point on the particle during a path search. As further described below, the Q* 
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search in the embodiments described below uses the particle filter functions to determine whether a 
potential expansion to a destination particle is a viable one (Le., whether the expansion provides a 
cheaper path to any portion of the destination particle than previously identified expansions to the 
destination particle). This search makes this determination for an expansion when it initially 
identifies the expansion and also later if it expands about it. 

Filtering is when a first PLF Fl (a filtered PLF) is compared with a second PLF F2 (a filter 
PLF) to determine whether any portion of the Fl needs to be discarded. The filter and filtered PLFs 
have overlapping domains (i.e., domain(Fl)C\domain(F2)*null). In the embodiments described 
below, filtering discards the portion of the filtered PLF Fl that is larger than the corresponding 
portion of the filter PLF (Le., discards every portion of Fl where F1(V)> F2(V)). Figure 33 
presents an example of filtering a first filtered PLF 3305 by a second filter PLF 3310, where both 
PLFs are defined across a line 3350. Each PLF is represented by its sequence of knots, which is 
sorted by the domain-offset values. Knots K1-K5 specify the first PLF while knots K6-K10 specify 
the second PLF 3310. As shown in this figure, this filtering discards portion 3315 of PLF 3305 as the 
PLF-values of this portion are larger than the PLF-values of the corresponding portion 3320 of filter 
PLF 3310. After the filtering operation, two portions 3325 and 3330 of the PLF 3305 remain. Knots 
K1-K3 and Kll specify the first remaining portion 3325, while knots K12 and K5 specify the second 
remaining portion (where knots Kl 1 and K12 are at the intersection of the PLF's 3305 and 3310). 

Filtering PLFs will be further described below. The filtering that is described below not only 
filters the first PLF, but also defines the second filter PLF to be the minimum of the first and second 
PLFs. A minimum of two PLFs Fl and F3 is another PLF F3 that specifies a PLF-value for each 
location V in the intersection of Fl's domain and F2's domain that is equal to the smallest PLF-value 
specified by the first and second PLF's for that location (Le., F3(V) = min(Fl(V),F2(V)). Figure 34 
illustrates the minimum PLF 3405 for the two PLFs 3305 and 3310 of Figure 33. The portion 3325 

and 3330 of this minimum function corresponds to the remaining portion of the filtered PLF 3305, 
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while the portion 3320 of this minimum function is from the original filter PLF 3310. Knots K1-K3, 
Kl 1, K8, K9, K12, and K5 specify the minimum function 3405. 

(iv) The Q* Path Search Process 
The route-generation process 2200 directs (at 2220) the Q* search process to identify and 
embed the lowest-cost path between specified source and target particles in the tessellated region. 
Figure 35 illustrates a Q* search process 3500 of some embodiments of the invention. As shown in 
this figure, the process 3500 initially sets (at 3502) a variable Current_Drop to null. As in the above- 
described A* search example, the process 3500 uses drops to represent path expansions. Specifically, 
this process represents an expansion from one topological particle (called a "start particle" or 
"current particle") to another topological particle (called a "destination particle") by a "drop" that 
associates with the destination particle and refers back to the drop of the start particle. One of 
ordinary skill will realize that other embodiments might not use drops or might implement drops 
differently. 

At 3504, the process 3500 initializes filter and H PLFs for each topological particle in the 
tessellated region that can serve as a start, intermediate, or target particle for a path search. For each 
such topological particle, the process 3500 maintains (l)the filter PLF to express the lowest path 
cost that the process has been able to identify from the source set to the particle during a path search, 
and (2) the H PLF to express the estimated distance between the particle and the target set. The 
process 3500 stores the H PLF for each particle so that it only has to compute the H PLF once for 
each particle. In some embodiments, the process initializes the filter and H PLFs to "infinite". Also, 
in the embodiments described below, nodes, exits, and holes are the topological particles that can 
serve as start particles during a path search, while nodes, holes, exits, and walls are the topological 
particles that can serve as intermediate and target particles during the path search. 

Next, for each source particle that the process 2200 specified for the current path search, the 
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process 3500 (at 3506) identifies and sets the particle's H PLF and sets the particle's filter PLF to 

zero. The generation of a particle's H PLF is described below. At 3506, for each source particle, the 
process 3500 also specifies a drop, defines the source particle as the drop's particle, defines the drop's 

prior drop as null, sets the drop's G PLF to zero, and sets the drop's F PLF equal to the source 

particle's H PLF. The process stores the specified drops in a storage structure. In the embodiments 
described below, the storage structure is a priority queue (e.g., a heap) that is sorted based on the 
minimum of the F PLF of each drop. 

At 3508, the process then retrieves from the priority queue a drop with the smallest minimum 

F value. Next, the process filters (at 3510) the retrieved drop's G PLF with the filter PLF of this 
drop's particle. As further described below, the process performs this filtering operation to ensure 
that the retrieved drop is still valid. After the process 3500 stored the retrieved drop in the priority 
queue, it might have created and stored other drops that represent cheaper expansions to the retrieved 
drop's particle. These cheaper expansions would have modified the filter PLF of the retrieved drop's 
particle. 

Figure 36 illustrates one such situation. This figure illustrates a G PLF 3605 of a first 
expansion to a line 3610. This G PLF is also the filter function of the line 3610 after the first 
expansion. Figure 36 illustrates a G PLF 3615 of a second expansion to the line 3610. The second 
expansion's G PLF 3615 is smaller than the first expansion's G PLF 3605 over the entire line 3610 
except for the portion 3625 of the line. After the second expansion is identified, the filter function of 
the line 3610 is set to the minimum of the second expansion's G PLF and the filter function's original 
PLF (which, as mentioned above, is identical to the first expansion's G PLF). This new filter function 
is the PLF 3620 in Figure 36. 

Accordingly, after retrieving a drop from the priority queue, the process 3500 filters (at 3510) 
the drop's G PLF with its particle's filter function to ensure that the drop still represents the best valid 
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expansion to one or more portions of its particle. In the example illustrated in Figure 36, the filtering 
of the G PLF 3605 of the first expansion to (Le. t the first drop on) the line 3610 with the line's filter 
PLF 3620 after the second expansion results in the PLF 3630. This PLF 3630 corresponds to the 
portion of the original PLF 3605 of the first expansion that is smaller than the second-expansion's 
5 PLF 3615. This PLF 3630 is defined over a much smaller domain (i.e., segment 3625) than the 
original PLF 3605 for the first expansion. Filtering one function against another will be further 
described below. 

At 3512, the process determines whether any portion of the CurrentJDrop's G PLF remains 
after the filtering operation. If the process determines (at 3512) that at least one portion of the 

10 CurrentJDrop's G PLF remains after the filtering, the process determines (at 3514) whether the 
filtering operation at 3510 resulted in two or more convex pieces of the retrieved drop's G PLF. 

If the filtering did not result in two or more convex pieces, the process specifies (at 3516) the 
retrieved drop's F PLF again, as some portions of this drop's G PLF might have been discarded 
because of the filtering operation, and such a modification would discard some portions of this drop's 

15 F PLF. Next, the process determines (at 3518) whether the retrieved drop's F PLF minimum is 
greater than the lowest F PLF minimum of the drops that are currently stored in the priority queue. 
If so, the process stores (at 3522) the retrieved drop again in the priority queue, and then transitions 
back to 3508 to select another drop. Otherwise, the process specifies (at 3520) the retrieved drop as 
the CurrentJDrop and then transitions to 3532, which will be further described below. 

20 If the process determines (at 3514) that the filtering at 3510 resulted in two or more convex 

pieces of the retrieved drop's G PLF, the process specifies (at 3524) a drop for each remaining piece 

and sets their parameters as follows. The process defines each specified drop's particle as the 

retrieved drop's particle. It also sets each specified drop's previous drop identically to the retrieved 

drop's previous drop (which might be null). The process also sets each specified drop's G PLF equal 

25 to the portion of the retrieved drop's G PLF for which the drop was specified. It also sets each 
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specified drop's F PLF equal to the sum of (1) the specified drop's G PLF, and (2) the portion of the 
retrieved drop's H PLF with the same domain as the specified drop's G PLF. 

At 3526, the process then determines whether any of the drops created at 3524 has the lowest 
F PLF minimum of all the drops stored in the priority queue. If not, the process stores (at 3528) the 
drops specified at 3524 in the priority queue based on the minimum of the F PLF of each drop. 
From 3528, the process transitions back to 3508 to select another drop. On the other hand, if the 
process determines (at 3526) that at least one drop specified at 3524 has the lowest F PLF minimum 
of all the drops stored in the priority queue, the process identifies (at 3530) as the Current_Drop a 
drop that was specified at 3524 and that has the lowest F PLF minimum. At 3530, the process also 

stores the remaining specified drops in the priority queue based on the minimum of the F PLF of 
each drop. 

From 3530, the process transitions to 3532. The process also transitions to 3532 from 3520. 
At 3532, the process determines whether the Current JDrop's particle is a target. If not, the process 
tries (at 3534) to expand the path search about the Current_Drop. Specifically, at 3534, the process 
tries to identify one or more potential expansions about the Current_Drop. Some embodiments 
identify potential expansions about the Current_Drop by (1) identifying a set of valid spaces to which 
a path can expand from the Current_Drop's particle, and (2) identifying destination particles in each 
identified space. A valid space is one that contains the Current_Drop's particle but does not contain 
the particle of the prior drop in a path search (i.e., does not contain the particle of the CurrentJDrop's 
previous drop). Figure 37 illustrates an example of a valid expansion space for a planar expansion. 
This figure illustrates a path search 3705 that has a last drop 3710 and a second-to-last drop 3715. 
The particles of both these drops are part of space 3725. The particle of drop 3710 is also part of 
space 3720. Consequently, space 3720 is a valid expansion space for drop 3710, but space 3725 is 
not. There can be multiple viable expansion spaces as a retrieved drop's particle (such as a node) can 
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be in multiple spaces. One of ordinary skill will realize that other embodiments might identify 
potential expansions about the Current_Drop differently. For instance, some embodiments might 
define a valid space for an expansion differently. One such embodiment might not require that the 
Current JDrop's particle to be part of the valid space. 
5 After the process identifies (at 3534) a set of valid spaces to which the path can expand from 

the CurrentJDrop's particle, it identifies (at 3534) potential destination particles in each identified 
space. In some embodiments, a path can expand towards exits, walls and holes, as well as the 
selected net's nodes, in a valid expansion space. One of ordinary skill will realize that other 
embodiments might specify other potential expansion particles in a valid expansion space. 

10 In the example illustrated in Figure 37, the path search can expand from drop 3710 to exit 

3745, hole 3750, and wall 3725. Wall 3725 belongs to a route 3730 of another net. The path 
generation process 3500 might allow a path expansion to the wall of another net's previously defined 
path because it might allow later defined routes to rip up earlier defined routes. In such cases, 
expanding to the walls of another net's path is assessed a penalty, as further described below. 

15 In some situations, the process 3500 cannot identify (at 3534) any potential expansions to 

another particle. However, when the process identifies one or more potential expansions at 3534, the 

process performs the following four operations at 3534. First, it specifies the H PLF of each 
potential expansion's destination particle if it had not previously been set for the current path search 

(i.e., for the current source and target sets). The generation of the H PLF is described below. 

20 Second, the process specifies a G PLF for each potential expansion. The generation of the G PLF for 

an expansion (i.e., the costing of an expansion) is described below. Third, it filters the G PLF of each 

potential expansion with the filter PLF of the expansion's destination particle. This filtering 

operation also sets the destination particle's filter PLF equal to the minimum of the filtered G PLF 

and the destination particle's previous filter PLF. Filtering two PLFs is described below. Fourth, the 

25 process specifies a drop for each convex piece (if any) of a G PLF of a potential expansion that 
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remains after the filtering operation. For each specified drop, the process (1) sets the drop's G PLF 
equal to the remaining piece of the filtered G PLF for which the drop was created, (2) associates the 
drop with its expansion's destination particle, (3) sets the drop's previous drop to the Current-Drop, 
and (4) sets the drop's F PLF to the sum of the drop's G PLF and the portion of its destination 
particle's H PLF that is defined over the domain of the drop's G PLF. The process stores each drop it 
creates at 3534 in the priority queue. 

From 3534, the process transitions to 3536. The process also transitions to 3536 if it 
determines at 3512 that no portion of a retrieved drop's G PLF remains after the filtering at 3510. At 
3536, the process determines whether the priority queue that stores the drops is empty. If so, the 
process has failed to find a path between the specified source and target sets. Accordingly, it returns 
(at 3538) a notification specifying its failure and then ends. On the other hand, when the process 
determines (at 3536) that the priority queue is not empty, the process transitions back to 3508 to 
retrieve the next lowest-cost drop from this priority queue and then to perform the subsequent 
operations for this drop. 

The process has found a path between the source and target sets when it determines (at 3532) 
that the Current_Drop's particle is a target. In this situation, the process transitions from 3532 to 
3540. At 3540, the process 3500 increments the Route_Cost by the cost of the identified path 
between the source set and the Current_Drop's particle. This cost is the minimum cost of the 
Current_Drop's G PLF. 

After 3540, the process topological^ embeds (at 3542) the identified path between the source 

set and the target. Starting at the CurrentJDrop on the target, the topological embedding "back 

traces" the sequence of drops that reached the target and generates an orderd list of topological 

particles that define the topological path. Generation of such an ordered list entails creation of wall 

particles between each pair of non-wall, non-via particles in the topological path, and can also entail 

the creation of joints, Steiners, and via nodes. Also, each time the topological-embedding process 
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defines a joint on an edge, it marks all other joints on the edge as "dirty," which signifies that their 
locations need to be later updated. The topological embedding process is further described below. 

After the topological embedding operation, the process triangulates (at 3544) further any face 
that contains a via node that was defined during the topological embedding operation. Specifically, 
5 the process embeds a topological path by replacing each hole with one via node on each layer 
traversed by the hole. The process then triangulates each face that contains one of the via nodes (Le. 9 
decomposes each face that contains a new via node into three triangles by defining an edge between 
the via node and each of the face's nodes). Via nodes that replace via holes are moveable in some 
embodiments. In these embodiments, they can be shoved during routability checking to alleviate 

10 congestion about them. 

At 3544, the process also might perform a follow-up edge-flipping operation on the new 
triangles that were created at 3544. This edge-flipping operation is similar to the edge-flipping 
operation at 935. In addition, at 3544, the process also creates one or more spaces for the new 
triangles and modifies the topological description of the routes of nets that cross the new tessellated 

15 edges. At 3544, the process marks the joints of the routes that cross the new triangulated edges as 
"dirty" so that their positions can be later updated. At 3544, the process also tries to create holes for 
the new spaces that it defined at 3544, and computes the capacity of each new edge that it created at 
3544. 

After 3544, the process defines (at 3546) a loose geometric path for the topological path that 
20 it specified at 3542. This entails specifying a location of each joint, via node, and Steiner on the path. 
At 3546, the process also re-computes the location of the joints that were marked dirty during the 
topological embedding and edge flipping operations of 3542 and 3544. This computation also affects 
the location of holes and Steiners. The geometric embedding operation is further described below. 

One of ordinary skill will realize that the Q* path-generation process might be implemented 
25 differently in other embodiments. For instance, some embodiments might utilize non-convex PLFs. 
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Also, in some embodiments, the H cost function might not specify a lower bound on the shortest 
path between a drop's particle and a target set. In addition, some embodiments might compute the F 
function slightly differently. For instance, some embodiments might express the F function as: 

F =G + 2*H. 

Such a function would bias the search process to expand about the drops that are closer to the target 
set. Alternative embodiments might express the F function as: 

F = G + H + J, 

where J represents the estimated computational effort needed to complete the path from the current 

drop. The embodiments that use alternative F function do not satisfy the admissibility requirement 
(Le., they do not produce consistently the lowest-cost path between source and target sets). On the 
other hand, the embodiments that use the F function of equation (4) do satisfy this requirement, 
(v) Example of the Q* Search Process 

Figure 38 presents an example of the Q* path search. In this example, the Q* search process 
is used to construct the lowest-cost path between a single source node 3805 and a set of target nodes 
3810 in a triangulated graph 3800. To simplify this example, the only source and target particles are 
nodes, and the only intermediate particles are edges. Also, this example accounts only for the 
distance costs, and ignores spacing constraints on edges due to obstacles. 

Figure 38 illustrates several sets of expansions identified during a path exploration phase. 
The expansions are represented by drops la, lb, 2, 3a, 3b, 4a, 4b, 5-7, 8a, 8b, 9, 10, and 11. In this 
example, the Q* path search starts by identifying and storing drops la and lb about the source node 
3805. This process then successively identifies the sequence of drops 2, 3a and 3b, and 4a and 4b. It 
identifies this sequence by successively (1) identifying each of the drops la, 2, and 3a as the drops in 

the priority queue with the lowest F function minimum, and (2) identifying viable expansions about 
each of the edges of these drops. Once the process creates and stores expansion drops 4a and 4b as 
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viable expansions from drop 3a, the process identifies drop lb as the drop with the lowest F 
function minimum in the priority queue. Accordingly, it then successively identifies the sequence of 
drops 5, 6, 7, 8a, 8b, 9, 10, and 11 that reaches the target set 3810. The process identifies this 
sequence by successively (1) identifying each of the drops 5, 6, 7, 8a, 9 as the drops in the priority 

queue with the lowest F function minimum, and (2) identifying viable expansions about the edge of 
each of these drops. 

In this example, two expansions are identified to each of the following four edges 3825, 
3830, 3835, and 3840. The search process identifies the first expansions (i.e., identifies drops 3a, 3b, 
4a, and 4b) to these edges when it expands about drops la, 2, 3a, 3b, and 4a, while it identifies the 
second expansions {i.e., identifies drops 6, 7, 8a, and 8b) to these edges when it expands about drops 
5, 6, and 7. The second set of drops (Le. 9 drops 6, 7, and 8b) that are identified for the edges 3825, 
3830, and 3840 are only specified over portions 3845, 3850, and 3855 of these edges. This is because 
the filtering operations that were performed to assess the viability of these second expansions to 
edges 3825, 3830, and 3840 resulted in the filtering of their G PLFs outside of the portions 3845, 
3850, and 3855. However, the second drop 8a on edge 3835 is defined over the entire edge, as the 
filtering operation that the process performs before identifying this drop does not filter any portion of 
this expansion's G PLF. This is because the entirety of the edge 3835 can be reached more cheaply 
through the sequence of drops lb, 5, 6, and 7 than through the sequence of drops la, 2, and 3a. 

In the example illustrated in Figure 38, the lowest-cost path is the one represented by drops 

11, 9, 8a, 7, 6, 5, lb, and 12, where drop 12 is the drop created for the source node 3805. After 

identifying this sequence of drops, the Q* search process embeds a path for this sequence. An 

approximation of this embedded path is illustrated in Figure 39. As shown in Figure 39, joints are 

defined at the intersection of the embedded path and the edges that this path intersects. 

(4) Computing G PLF 

When the process 3500 identifies (at 3534) a potential expansion from the Current_Drop's 
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particle to a destination particle, this process specifies (at 3534) a G PLF for such a potential 
expansion. In some embodiments, the process 3500 computes this G PLF by (1) propagating the 
Current_Drop's G PLF to the destination particle, and (2) for certain expansions, adding one or more 
penalty costs to the PLF resulting from the propagation. 
5 (i) Propagation 

Propagating the Current_Drop's G PLF to an expansion's destination particle generates a PLF 
that expresses the distance and penalty costs of reaching the CurrentJDrop's particle plus the distance 
cost of reaching the expansion's particle from the Current.Drop's particle. In other words, the 
propagation accounts for the extra distance cost of going from the Curoent_Drop's particle to the 
10 expansion's particle. 

Some embodiments limit each propagation operation to only the portion of the expansion's 
destination particle to which a path can expand from the CurrentJDrop. The discussion below uses 
the phrase "destination state" to refer to the portion of an expansion's destination particle to which a 
path can expand from the CurrenMDrop. This discussion also uses the phrases "the Current_Drop's 
15 domain" or "the start domain" to refer to the domain of the Current_Drop's G PLF. 

Nine different propagation operations are described below for nine pair-wise combinations of 
expansions between points, line, and surfaces. In these propagation operations, points represent 
nodes, lines represent edge- and wall-segments to which a path can expand, and surfaces represent 
portions of holes to which a path can expand. 
20 As mentioned above, a hole is specified between each pair of overlapping spaces when the 

intersection of the two polygons of the overlapping spaces is larger than a threshold size. Such a hole 
represents a potential via that can be anywhere in the polygonal intersection of the two polygons of 
the two overlapping spaces. In some embodiments, a path for a net can expand to each portion of a 
hole. These embodiments might specify the same width and spacing constraints for all nets being 
25 routed at a given time, or might construct each hole for the worst case of constraints (e.g., construct 
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the polygons of the hole's overlapping spaces for the worst case of constraints). 

Alternatively, some embodiments construct each hole based on the best case of constraints 
(e.g., construct the polygons of the hole's overlapping spaces for the best case of constraints). For an 
expansion of a path of a net to a hole, these embodiments then identify a sub-set of the hole's 
5 polygon (i.e., the polygonal surface represented by the hole) to which the net's path can expand given 
the spacing constraints for the net. These embodiments identify such a sub-set by identifying within 
the hole's polygon a polygon that specifies the legal area for the center of a via for the net. These 
embodiments identify such a polygon within a hole polygon in the same manner that they identified 
the hole polygon in the first place. They identify the closest point within the hole's polygon that a via 

10 can be to each obstacle on the layer, accounting for the via size and/or shape, the required spacing to 
obstacles, the required spacing to other routes crossing the face, etc. If the result is disjoint or non- 
convex, these embodiments approximate the polygon as the convex hull of the legal locations. In 
some cases, such a polygon might be identical to the hole's polygon. 

(1st) Expansion from point to point 

15 To perform a propagation operation for an expansion that goes from an expansion start point 

to an expansion destination point (e.g., goes from one node to another node), the process 3500 
identifies the distance between the two points and adds this distance to the G PLF of the start point 
(i.e., to the Current_Drop's G PLF). In some embodiments, the process measures the distance 
between two points in an IC layout using only a specified set of wiring directions that are specified 

20 by the layout's wiring model. For instance, when two points are on an octilinear layer (like any of the 

layers of Figure 1) that can have horizontal, vertical, and ±45° diagonal interconnect lines, the 

distance between the two points is the shortest distance that can be traversed by using horizontal, 

vertical, and ±45° diagonal interconnect lines. As disclosed in United States Patent Application 

entitled "Hierarchical Routing Method and Apparatus that Utilize Diagonal Routes," and having 

25 serial number 10/013,813, the shortest distance between two points for such an octilinear wiring 
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model: 

Distance = L + S * (sqrt(2) - 1), (5) 
where L and S respectively represent the long and short sides of a rectangular bounding box that is 
aligned with the layout's x-y axes and that encompasses the two points. This manner of computing 
5 the shortest distance does not disfavor or penalize any preferred wiring direction over another 
preferred wiring direction. 

Numerous other operations below require the computation of the distance between two 
points. The embodiments described below compute each such distance according to the above- 
described equation (5). 

1 0 (2nd) Expansion from point to line or point to surface 

Figure 40 illustrates a process 4000 that propagates a PLF that is defined over a point to a 
line or a surface. This process is described below by reference to Figure 41, which illustrates an 
example of propagating a PLF from a point P to a line L, and Figure 42, which illustrates an 
example of propagating a PLF from a point P to a surface S. 

15 As show in Figure 40, the process 4000 initially projects (at 4005) vectors in all available 

interconnect directions from the Current_Drop's point. These vectors are referred to below as 
propagation vectors. The embodiments described below utilize the wiring model illustrated in 
Figure 1. Accordingly, these embodiments project vectors in the 0°, 45°, 90°, 135°, 180°, 225°, 
270°, and 315° directions. Figures 41 and 42 illustrate eight propagation vectors emanating from 

20 their points P in 0°, 45°, 90°, 135°, 180°, 225°, 270°, and 315° directions. Some embodiments do not 
project vectors in all interconnect directions. Instead, they project only the propagation vectors that 
will intersect the destination state. These propagation vectors are the vectors that fall within a 
triangle defined by the start-state point and the leftmost and rightmost vertices of the destination 
state. 

25 At 4010, the process then identifies the intersection of the projected propagation vectors and 
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the destination line or surface. A propagation vector intersects with a line at a point that is the 
location of a knot of the destination's line PLF. Figure 41 illustrates the intersection of the 
propagation vectors 4105, 4110, and 4115 at three points 4120, 4125, and 4130 along line L. A 
propagation vector intersects a surface either (1) at only a vertex of the surface, or (2) at two points 
5 on the surface boundary and an edge that runs through the surface. Figure 42 illustrates the 
intersection of two propagation vectors 4215 and 4220 with the surface S along edges 4270 and 
4275, which respectively terminate on boundary point pair 4225 and 4230 and point pair 4235 and 
4240. 

Next, for the expansion to the destination state, the process 4000 then specifies (at 4015) a G 
10 PLF with knots at the intersections identified at 4010. Specifically, for a destination line, the process 
4000 specifies (at 4015) a knot at each point of the destination line that a propagation vector 
intersects. Figure 41 illustrates a G PLF 4145 with three knots 4150, 4155, and 4160 for the three 
intersection points 4120, 4125, and 4130. For a destination surface, the process 4000 specifies (at 
4015) a knot at each surface vertex intersected by a propagation vector and a knot at each surface 
15 boundary point intersected by a propagation vector. Figure 42 illustrates a G PLF 4260 with four 
knots 4280, 4282, 4284, and 4286 that are defined for four boundary intersection points 4225, 4230, 
4235, and 4240. 

At 4015, the process 4000 sets the PLF-value of each knot that it specifies at 4015. The PLF- 
value of a knot specified at 4015 equals (1) the Current_Drop's PLF-value at the start point, plus 
20 (2) the distance between the knot's x,y coordinates and the start point P, where this distance is 
measured along the direction of the propagation vector that was used to identify the knot. For 
example, in Figure 41, the PLF-value of knot 4155 equals the PLF-value at point P plus the distance 
D between point P and intersection point 4125 along the direction of propagation vector 4110. In 
Figure 42, the PLF-value of knot 4286 equals the PLF-value at point P plus the distance D between 
25 point P and intersection point 4240 along the direction of propagation vector 4220. 
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At 4020, the process specifies knots for the expansion's G PLF at the location of any 
unexamined vertex of the destination state. The vertices of a destination line are its endpoints. A line 
endpoint is an unexamined vertex if none of the projected propagation vectors intersect the 
destination line at the endpoint. The PLF-value of an unexamined endpoint of a destination line 
equals (l)the Cument_Drop's PLF-value at the start point, plus (2) the distance (according to the 
equation (5)) between the start point and the endpoint. For instance, in Figure 41, the process 
specifies two knots 4165 and 4170 at the endpoints 4175 and 4180 of the line L, and specifies the 
PLF-value of each of these knots as the PLF-value at point P plus the distance (according to equation 
(5)) between point P and each knot location 4175 or 4180. 

The vertices of a destination surface are the vertices of the edges that define the surface's 
external boundary. Such a vertex is an unexamined vertex if none of the propagation vectors 
intersect the destination surface at the vertex. In Figure 42, the destination surface S has five vertices 
4265, all of which are "unexamined" as they are not at the intersection of the propagation vectors and 
the surface. Accordingly, five knots 4290 are specified at the unexamined vertices 4265. The PLF- 
value for each such knot equals (1) the distance (according to equation (5)) between the start point 
and the surface's vertex that corresponds to the knot, plus (2) the Current_Drop's PLF-value at the 
start point. For instance, the PLF-value of knot 4290a equals the PLF-value at the start point P plus 
the distance (according to equation (5)) between point P and vertex 4265a. 

For a destination surface, the process 4000 uses (at 4025) the knots specified at 4015 and 
4020 to specify the edges and facets of the surface PLF that is defined over the destination surface. 
For each facet, the process defines a normal and a z-intercept. Standard plane-defining techniques 
can be used to derive the normal and z-intercept values of a facet from three points on the facet. For 
each facet, the process 4000 identified at least three knots at 4015 and 4020. 

(3rd) Line to Line or Surface to Line 

Figure 43 illustrates a process 4300 for propagating a PLF from a line to another line or from 
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a surface to a line. This process is described by reference to Figure 44, which illustrates the 
propagation of a line PLF from a line 4405 to another line 4410, and Figure 45, which illustrates the 
propagation of a surface PLF from a surface 4502 to a line 4504. 

As shown in Figure 43, the process 4300 initially identifies (at 4305) the propagation vectors 
5 that emanate from the locations on the CurrentJDrop's domain that are locations of knots in the 
CurrentJDrop's G PLF. The identification of these propagation vectors are further described below 
by reference to Figures 46-51. In Figure 44, knots are located at points 4415 and 4420 on line 4405. 
In Figure 45, knots are located at vertices 4506-4518 of surface 4502. 

Next, at 4310, the process 4300 projects the propagation vectors identified at 4305. Figure 

10 44 illustrates the projection of five propagation vectors from knot-location 4415 and five propagation 
vectors from knot-location 4420. In Figure 45, three propagation vectors are projected from each of 
the vertices 4506, 4512, and 4514, two propagation vectors are projected from each of the vertices 
4508 and 4518, and one propagation vector is projected from each of the vertices 4510 and 4516. In 
some embodiments, the process 4300 does not project (at 4310) vectors in all interconnect directions. 

15 Instead, it projects only the propagation vectors that will intersect the destination state. In these 
embodiments, this process projects propagation vectors that fall within a triangle defined by the 
knot-location and the leftmost and rightmost vertices of the destination line or surface. 

The process 4300 identifies (at 4315) the intersection of the destination line and the 
propagation vectors that were projected at 4310. Figure 44 illustrates the intersection of the 

20 propagation vectors 4430 and 4435 and the line 4410 at points 4450 and 4455. Figure 45 illustrates 
the intersection of the propagation vectors 4522, 4524, and 4526 and the line 4504 at points 4530, 
4532, and 4534. 

Each intersection identified at 4315 is a knot in the expansion's G PLF. Accordingly, for the 

expansion to the destination line, the process specifies (at 4320) a line G PLF with a knot at each 

25 intersection identified at 4315. At 4320, the process computes the PLF- value of each knot specified 
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at 4320. The PLF-value of each destination-state knot equals (l)the PLF-value of the 
Current_Drop's knot that was used to identify the destination-state knot, plus (2) the distance 
between the x,y coordinates of the Current_Drop and destination-state knots, where this distance is 
measured along the projected propagation vector that identified the destination-state knot. 

For instance, in Figure 44, the PLF-value of knot 4472, which is defined at the intersection 
4450 of the propagation vector 4430 and the line 4410, equals (l)the PLF-value of the 
Current_Drop's G PLF at 4420 on the line 4405, plus (2) the distance Dl between points 4420 and 
4450 along the direction of the propagation vector 4430. Similarly, in Figure 45, the PLF-value of 
knot 4544, which is defined at the intersection 4532 of the propagation vector 4524 and the line 
4504, equals (l)the PLF-value of the Current.Drop's G PLF at 4512 on the surface 4502, plus 
(2) the distance D2 between points 4512 and 4532 along the direction of the propagation vector 
4524. 

At 4325, the process specifies knots for the expansion's G PLF at the location of unexamined 
vertices of the destination line, and then terminates. As mentioned above, an unexamined vertex of a 
destination line is an endpoint of the line that does not intersect any of the projected propagation 
vectors. An unexamined destination-line vertex can be in either (1) a "wedge" that is defined by two 
propagation vectors that emanate from the same location on the Current_Drop's domain, (2) a 
"channel" that is defined by two parallel propagation vectors that emanate from two different 
locations on the Current_Drop's domain. The PLF-value of a knot specified for an unexamined 
vertex that is within a wedge defined by two propagation vectors emanating from the same start-state 
location equals (1) the PLF-value of Current_Drop's G PLF at the start-state location, plus (2) the 
distance (according to the equation (5)) between the start-state location and the unexamined vertex. 
On the other hand, the PLF-value of a knot specified for an unexamined vertex that is within a 
channel defined equals (1) the length of a line segment that is parallel to the two channel-defining 
vectors and that starts at the unexamined vertex and terminates at the start domain, plus (2) the PLF- 
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value of the Current_Drop's G PLF at the point that the line segment terminates on the start domain. 
The line segment terminates on the start domain on a second line segment that is between the two 
knot locations from which the two channel-defining vectors emanate. When the start domain is a 
surface, the second line segment (1) is an edge on the boundary of the surface if the two knot 
5 locations are boundary vertices of the surface, and (2) is a line segment within the surface if the two 
channel-defining knot locations are within the surface. 

For instance, in Figure 44, endpoint 4445 of the line 4405 falls within a channel defined by 
propagation vectors 4425 and 4430. The distance between endpoint 4445 and line 4405 is the length 
D3 of a line segment 4440 that is parallel to vectors 4425 and 4430. Accordingly, the PLF-value of 

10 the knot 4476 specified at endpoint 4445 equals the length D3 plus the PLF-value of the 
Current_Drop's G PLF at point 4480, which is the location that line segment 4440 intersects line 
4405. The endpoint 4460, on the other hand, is within a wedge defined by two propagation vectors 
4435 and 4482 that emanate from the start-state location 4420 on the line 4405. Accordingly, the 
PLF-value for the knot 4478 that is specified at 4460 equals (1) the PLF-value of the Current_Drop's 

15 G PLF at the start-state location, plus (2) the distance between points 4460 and 4420 according to the 
equation (5). 

In Figure 45, endpoint 4554 falls within a channel defined by propagation vectors 4520 and 
4522. The distance between endpoint 4554 and surface 4502 is the length D4 of a line segment 4560 
that is parallel to vectors 4520 and 4522. Accordingly, the PLF-value of the knot 4550 specified at 

20 point 4554 equals the length D4 plus the PLF-value of the Current_Drop's G PLF at point 4562, 
which is the location that the line segment 4560 intersects the surface boundary. The endpoint 4556, 
on the other hand, is within a wedge defined by two propagation vectors 4526 and 4564 that emanate 
from the start-state location 4512. Accordingly, the PLF-value for the knot 4552 that is specified at 
4556 equals (1) the PLF-value of the CurrentJDrop's G PLF at the start-state location 4512, plus 

25 (2) the distance between points 4512 and 4556 according to the above-described equation (5). 
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Figures 46-51 illustrate how some embodiments identify the propagation vectors that 
eman'ate from the knot locations of a line PLF or surface PLF. These embodiments identify the 
propagation vectors based on the following observations. At most eight propagation-vector wedges 
can be defined about a CurrentJDrop's domain when the octilinear wiring model such as the one 
illustrated in Figure 1 is used. Knots can be the only sources for wedges. 

As mentioned above, each wedge is defined by two propagation vectors that emanate from 
the same knot location. Two wedges are abutting wedges when they share a propagation vector, 
while two wedges are adjacent wedges when they have parallel propagation vectors. For instance, 
Figure 46 illustrates eight wedges 4650-4664 that are defined about five knot locations 4625-4645 
from a surface PLFs domain 4600. In this figure, there are three abutting wedge pairs, and five 
adjacent wedge pairs. For instance, wedges 4658 and 4656 are abutting wedges as they share vector 
4666, while wedges 4654 and 4656 are adjacent wedges as their vectors 4668 and 4670 are parallel. 

The parallel propagation vectors of two adjacent wedges define a freeway. For instance, a 
freeway 4672 exits between adjacent wedge pairs 4662 and 4664 in Figure 46. This freeway either 
(1) defines a channel when no other propagation vectors emanates from a knot location that falls 
within the freeway, or (2) contains two or more channels when other propagation vectors that are 
parallel to the freeway emanate from knot location(s) that fall(s) within the freeway. At most there 
are eight adjacent wedge pairs about a Current_Drop's domain. Consequently, there are at most eight 
freeways between the adjacent wedge pairs. 

Some embodiments define the direction of the propagation vectors that emanate from an a 
CurrentJDrop's domain by performing the following three operations. First, they identify the knot 
location for each wedge. Second, these embodiments identify one or more freeways between 
adjacent wedge pairs. Third, for each identified freeway, these embodiments determine whether to 
treat the freeway as one channel, or to define additional channels within the freeway by defining one 
or more propagation vectors that are parallel to the freeway-defining vectors and that emanate from 
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knot location(s) that fall within the freeway. The first operation (/.«., the identification of the knot 
location for each wedge) is described below by reference to Figure 47 and 48. The second and third 
operations (i.e., the identification of freeways and additional channels within the freeways) are then 
described by reference to Figures 49-51. 

5 Figure 48 illustrates a process 4800 that identifies the knot location for each propagation- 

vector wedge that emanates from the Current_Drop's domain, which can be a line or a surface. This 
process will be described by reference to an example illustrated in Figure 47. This example 
illustrates how one of the wedges of the PLF-surface domain 4600 of Figure 46 is identified. 

The process 4800 initially selects (at 4805) a wedge. There are eight wedges when the 

10 octilinear wiring model of Figure 1 is used. The process then selects (at 4810) a knot location as the 
first candidate location for the source of the selected wedge, and records this knot location as the 
current best (Current_Best) location for the selected wedge. Figure 47 illustrates a selected wedge 
4660 and a first candidate knot location 4610. Next, at 4815, the process selects a second candidate 
knot location for the source of the selected wedge and records it as Next_Candidate location. In 

15 Figure 47, the second candidate location is location 4625. 

At 4820, the process then determines whether the Current_Best knot location is a better 
candidate location than the Next_Candidate location. This determination ensures that all locations in 
the wedge are best reached from the source of that wedge, when both distance and PLF-value is 
considered. To make this determination, the process performs the following four operations. First, 

20 the process uses the selected wedge's "vector," which in some embodiments is defined to be a unit 
vector that bisects the wedge (i.e., it is midway between the two vectors that define the wedge). 
Figure 47 illustrates the vector 4705 of the selected wedge 4660. 

Second, the process computes the dot product of the selected wedge's vector and a vector that 
starts from Current_Best knot location and terminates on Next_Candidate knot location. This 

25 computation quantifies whether the Next_Candidate knot location is closer to an arbitrary point in 
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the wedge than the Currentjtest knot location. Figure 47 illustrates a vector 4710 from the 
Current_Best location 4610 to the Next_Candidate location 4625. 

Third, after computing the dot product, the process computes a Cost_Delta, which is the 
difference in the PLF-values of Next_Candidate location and the Current JBest location according to 
5 the Current_Drop's G PLF (CostJDelta = G PLF(Next_Candidate) - G PLF(CurrentJBest)). For the 
example in Figure 47, the CostJDelta is the difference in the PLF-values at locations 4625 and 4610. 

Fourth, the process determines whether the computed dot product is greater than the 
Cost_Delta. If not, the Current_Best location is better than the Next_Candidate location, and the 
process transitions from 4820 to 4830, which is further described below. If so, the Next_Candidate 

10 location is a better location for the wedge than CurrentJBest location, and the process transitions 
from 4820 to 4825. At 4825, the process sets the CurrentJBest location equal to the Next_Candidate 
location, and then transitions to 4830. At 4830, the process determines whether it has examined all 
knot locations for the wedge selected at 4805. If not, the process selects (at 4835) another knot 
location, sets this knot location as the Next_Candidate, and transitions to 4820 to compare this newly 

15 selected Next_Candidate with the Current_Best. 

When the process determines at 4830 that it has examined all the locations of knots in the 
Current_Drop's G PLF, the process defines (at 4840) CurrentJBest as the knot location for the 
selected wedge. The process then determines (at 4845) whether it has examined all the wedges. If 
not, it returns to 4805 to select another wedge and to repeat its operations for this wedge. Otherwise, 

20 it ends. 

After identifying the locations of the wedges, the Q* engine has to specify the channels 
between the wedges. When two wedges abut (i.e., when they share a propagation vector), no channel 
can exist between the wedges. However, when two wedges are adjacent wedges {i.e., when they have 
parallel propagation vectors), one or more channels can be defined in the freeway defined by the 
25 parallel vectors of the adjacent wedge pairs. 
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If the CurrentJDrop's domain is a line, the Q* engine examines each adjacent wedge pair that 
is defined along the line. If no knot location exists on the line segment between the adjacent wedge 
pair, then the Q* engine defines the freeway between the adjacent wedge pair as a channel. If one or 
more knot locations exist on the line segment between an adjacent wedge pair, the Q* engine 
5 examines each knot location to determine whether they are sources of propagation vectors. The 
engine initially sorts in descending order the PLF-values of the knot locations between the adjacent 
wedge pair. If there is only one knot location between the adjacent wedge pair, the engine simply 
adds this PLF-value to its sorted list. The engine then selects the largest PLF-value on the sorted list. 
The Q* engine then identifies the channel that contains the knot location corresponding to the 
10 selected PLF-value. This channel is the freeway formed by the adjacent wedge pair when the knot 
location is the first location between the wedge pair that is examined. When the knot location is not 
the first location between the wedge pair, the channel might be formed by one or two propagation 
vectors that the engine specified for knot locations (between the wedge pair) that it previously 
examined. 

15 The engine next determines whether the selected PLF-value for the knot location is less than 

the PLF-value that can be obtained by linearly interpolating between the values at the knot locations 
of the vectors that define the channel that contains the selected PLF- value's knot location (Le. t 
determines if the selected value is below a line that connects the PLF-values at the knot locations 
from which the identified-channel's vectors emanate). If so, the engine specifies a channel-defining 

20 vector at the knot location associated with the selected PLF-value. This specified vector is parallel to 
the parallel vectors of the adjacent wedge pair. 

Figure 49 illustrates an example of a knot location 4930 that is between two adjacent wedge 
pairs (where one pair is formed by wedges 4910 and 4915 and one pair is formed by wedges 4920 
and 4925). This knot location is examined for each of these adjacent wedge pairs. When the engine 

25 examines location 4930 for adjacent wedges 4910 and 4915, it determines whether the PLF-value of 
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this location is smaller than the PLF-value that can be obtained for this location by linearly 
interpolating between the PLF-values at knot locations 4945 and 4950 that serves as the emanating 
location of the wedges 4910 and 4915. In Figure 49, it is assumed that the PLF-value of location 
4930 is less than the PLF-value that can be obtained through the linear interpolation. Accordingly, 
5 two propagation vectors 4935 and 4940 are defined for the knot location 4930. The propagation 
vector 4935 is parallel to the parallel vectors of adjacent wedges 4910 and 4915, while the 
propagation vector 4940 is parallel to the parallel vectors of adjacent wedges 4920 and 4925. 

Figure 50 illustrates a process 5000 for defining channels between adjacent wedge pairs 
when the CurrentJDrop's domain is a surface. This process will be described by reference to 

10 Figure 51A, which illustrates the domain 5100 of a surface PLF. The process 5000 initially 
identifies (at 5005) all adjacent wedge pairs. It then selects (at 5010) one of the adjacent wedge pairs. 
In Figure 51 A, the selected pair of adjacent wedges are wedges 5105 and 51 10. 

Next, at 5015, the process projects each location of a knot (in the Current JDrop's G PLF), 
which falls within the freeway that the selected wedge pair defines, towards a line that connects the 

15 knot locations from which the wedge pair emanates. The projection is in a direction that is parallel to 
the adjacent-wedge-pair vectors that define the freeway. In Figure 51A, knot locations 5120, 5125, 
5130, 5135, 5140, and 5142 fall within the freeway 5115 between the adjacent wedges 5105 and 
5110. Also, the two wedges emanate from knot locations 5145 and 5150, and an edge 5170 exists 
between these two locations. Knot location 5120 lies on the edge 5170, and hence does not need to 

20 be projected towards this line segment. On the other hand, knot locations 5125, 5130, 5135, 5140, 
and 5142 need to be projected towards this line segment in the direction of the vectors 5155 and 
5160 that define the freeway 5115. Figure 51A illustrates the projection of knot location 5125 
towards this edge. 

At 5015, the process then computes a cost at the intersection of each projection with the line. 
25 The cost at each intersection point equals the sum of the PLF-value of the G PLF at the intersection 
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point plus the distance between the projected knot location and the intersection of the knot-location's 
projection and the line. In Figure 51 A, the intersection of the projection of knot location 5125 is 
point 5165 on the edge 5150. Its computed cost equals the distance D between knot location 5125 
and point 5165 and the value of the G PLF at the point 5165. The distance D between a projected 
5 knot location and the intersection of its projection and the line might be a positive or a negative value 
depending on position of the knot location with respect to the line that connects the knot locations 
from which the adjacent wedge pair emanates. For instance, Figure 51B illustrates an edge 5175 
within the domain 5180 of a surface PLF. This edge connects to adjacent wedge pair 5182 and 5184. 
Figure 51B also illustrates two knot locations 5186 and 5188 that fall within the freeway 5190 

10 (defined by the wedge pair 5182 and 5184) on opposing sides of the edge 5175. The distance Dl 
between knot location 5186 and this knot's projection point 5192 on the edge 5175 is positive, but 
the distance D2 between the knot location 5188 and this knot's projection point 5194 on the edge 
5175 is negative. The distance D2 is negative because the knot location 5188 is on the side of the 
edge 5175 that is in the direction of the freeway-defining wedge vectors. 

15 At 5020, the process sorts the computed PLF- values in descending order. It then selects (at 

5025) the largest PLF-value on the sorted list. The process then identifies (at 5030) the channel that 
contains the intersection point corresponding to the PLF-value selected at 5025. This channel is the 
freeway that the adjacent wedge pair selected at 5010 defines at the first iteration through 5030 for 
the selected wedge pair. In the subsequent iterations through 5030, the propagation vectors that the 

20 process 5000 might identify for the knot locations between the adjacent wedge pair might define this 
channel partially or completely. 

After 5030, the process determines whether the selected PLF-value for the intersection point 
is less than the PLF-value that can be obtained by linearly interpolating between the values at the two 
knot locations of the vectors that define the channel that contains the intersection point (£.*, 

25 determines if the selected value is below a line that connects the PLF-values at the knot locations of 
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the identified-channel ! s vectors). For instance, in Figure 51 A, the process 5100 has the PLF- value at 
the knot locations 5145 and 5150. From these two PLF- values the process can linearly interpolate the 
value of any point on the line 5150 between these two locations. Accordingly, if the first selected 
PLF-value is the PLF-value at the intersection point 5165 for the knot location 5125, the process 
determines (at 5130) whether the cost computed for this knot location's projection onto point 5165 is 
less than the value that can be obtained for this point by linearly interpolating between the values 
specified at knot locations 5145 and 5150. 

If the process determines (at 5135) that the selected PLF value is not less than the value that 
can be obtained through the linear interpolation, the process transitions to 5045, which is further 
described below. Otherwise, the process specifies (at 5040) a channel-defining vector at the knot 
location associated with the selected PLF-value. This specified vector is parallel to the parallel 
vectors of the adjacent wedge pair. From 5040, the process transitions to 5045. 

At 5045, the process determines whether it has examined all the knot locations that fall 
within the freeway defined by the adjacent wedge pair. If not, the process selects (at 5050) the next 
largest PLF-value in the sorted list that it created at 5020. The process then performs 5030 and 5035 
for the intersection point that corresponds to the PLF-value selected at 5050. As mentioned above, 
the first iteration through 5030 identifies the freeway as the channel that contains the intersection 
point corresponding to the PLF-value selected at 5025. However, a subsequent iteration through 
5030 from 5050 might identify a different channel that contains the intersection point corresponding 
to the PLF-value selected at 5050. This is because, if the process iterates through 5040 one or more 
times, it defines more propagation vectors that break the freeway into smaller and smaller channels. 
When the process 5000 identifies (at 5030) a smaller channel (Le. 9 a channel that covers only a 
portion of the freeway) that contains the intersection point corresponding to the selected PLF-value, 
it derives (at 5035) the interpolated value based on PLF-value of the knot locations from which the 
vectors that define the channel emanate. 
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When the process determines at 5045 that it has examined all the knot locations that fall 
within the selected wedge pair f s freeway, the process determines (at 5055) whether it has examined 
all the adjacent wedge pairs. If not, the process transitions back to 5010 to select an unexamined 
wedge pair. Otherwise, the process terminates. 

5 (4th) Line to Surface 

Figure 52 illustrates a process 5200 for propagating a G PLF from a line to a surface. The 
process 5200 is quite similar to the process 4300 of Figure 43. Accordingly, similar reference 
numbers are used for similar operations of the two processes. The process 5200 has only a few minor 
differences with the process 4300. Propagating a PLF to a surface results in a surface PLF. Hence, 

10 after 4325, the process 5200 uses the attributes of the knots that it specifies to define edges and facets 
of a surface PLF, and to define the normal and z-intercept values of the facets. 

Also, in process 5200, a propagation vector intersects the destination surface either at only a 
vertex or along an edge that runs through the surface and connects two points on the boundary of the 
surface. For the case where the propagation vector intersects the surface only at a vertex, the process 

15 5200 would specify (at 4320) a knot at the vertex's x,y coordinates. For the case where the 
propagation vector runs through the surface, the process 5200 would specify (at 4320) two knots at 
the surface boundary points where the propagation vector intersects the surface boundary. 

Figure 53 presents an example that illustrates the propagation of a G PLF from a line 5305 to 
a surface 5310. For such a propagation, the process 5200 would initially identify (at 4305) the 

20 propagation vectors that emanate from the locations on the CurrentJDrop's line that are locations of 
knots in the Current_Drop's G PLF. These propagation vectors are identified based on the process 
described above by reference to Figure 46-51. In Figure 53, the knots are located at points 5315 and 
5320 on line 5305. From each of these points, five propagation vectors are projected. From these 
points, some embodiments would project only the propagation vectors that would intersect the 

25 destination surface, as mentioned above. 



Next, at 4310 and 4315, the process 5200 projects the propagation vectors identified at 4305, 
and identifies their intersections with the destination surface. In Figure 53, the propagation vectors 
5325 and 5330 that emanate from points 5315 and 5320 intersect the destination surface along edges 
5360 and 5362, which respectively terminate on boundary point pair 5335 and 5340 and point pair 
5 5345 and 5350. 

The process 5200 then starts to specify (at 4320) a G PLF that is defined over the destination 
surface. Specifically, at 4320, the process specifies a knot at each intersection of the propagation 
vectors and the boundary of the destination surface. In Figure 53, the process starts to specify a 
surface PLF 5398 by specifying four knots 5352, 5354, 5356, and 5358 at the identified intersections 

10 5335, 5340, 5345, and 5350. The PLF-value of each specified destination knot equals (1) the PLF- 
value of the start knot that was used to identify the destination knot, plus (2) the distance between the 
x,y coordinates of the start and destination knots, where this distance is measured along the projected 
propagation vector that identified the destination knot. For instance, the PLF-value of knot 5354 that 
is specified for the location 5340 equals the distance Dl between 5340 and 5315, plus the PLF-value 

15 of the Current_Drop's G PLF at 53 15. 

At 4325, the process 5200 then specifies knots for the destination's PLF at the location of the 
unexamined vertices of the destination. In Figure 53, the unexamined vertices of the destination 
surface are vertices 5364, 5366, 5368, and 5370. For these vertices, the process 5200 specifies knots 
5372, 5374, 5376, and 5378. The PLF-value of each of these knots is computed based on whether the 

20 knot's corresponding destination-surface vertex falls within a propagation-vector wedge or channel. 

For instance, the unexamined vertex 5370 falls within a wedge defined by two propagation 
vectors 5325 and 5322 that project from 5315. Accordingly, the PLF-value of knot 5378 at vertex 
5370 equals (l)the PLF-value of the CurrentJDrop's G PLF at point 5315, plus (2) the distance 
(according to equation (5)) between points 5315 and 5370. On the other hand, the unexamined vertex 

25 5368 falls within a channel defined by two parallel propagation vectors 5325 and 5330 that project 
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from two different points 5315 and 5320. The distance between the vertex 5368 and the line 5305 is 
the length D2 of a line segment 5380 that is parallel to vectors 5325 and 5330. This line segment 
intersects line 5305 at 5382. Hence, the PLF-value of knot 5376 at vertex 5368 equals (1) the length 
D2, plus (2) the Current_Drop's G PLF-value at 5382. 

5 After 4325, the process 5200 uses (at 5205) the knots specified at 4320 and 4325 to specify 

the edges and facets of the surface PLF that is defined over the destination surface. For each facet, 
the process defines a normal and a z-intercept. Figure 53 illustrates three facets and ten edges that 
define these three facets. 

(5th) Surface to Point or Line to Point 

10 Figure 54 illustrates a process 5400 for propagating a PLF from a line to a point or from a 

surface to a point. This process is described by reference to Figure 55, which illustrates the 
propagation of a line PLF from a line 5505 to a point 5520, and Figure 56, which illustrates the 
propagation of a surface PLF from a surface 5602 to a point 5620. 

As shown in Figure 54, the process 5400 initially identifies (at 5405) the propagation vectors 

15 that emanate from the locations on the Current_Drop's domain that are locations of knots in the 
Current_Drop's G PLF. The identification of these propagation vectors was described above by 
reference to Figures 46-51. In Figure 55, knots are located at points 5510 and 5515 on line 5505. In 
Figure 56, knots are located at vertices 5604-5616 of surface 5602. 

Next, at 5410, the process 5400 projects the propagation vectors identified at 5405. 

20 Figure 55 illustrates the projection of six propagation vectors from knot-location 5510 and four 
propagation vectors from knot-location 5515. In Figure 56, three propagation vectors are projected 
from each of the vertices 5606 and 5614, two propagation vectors are projected from each of the 
vertices 5604, 5608, 5610, and 5612, and one propagation vector is projected from vertex 5616. 

The process 5400 then identifies (at 5415) the propagation-vector wedge or channel that 

25 contains the destination point. As mentioned above, a propagation-vector "wedge" is defined by two 



propagation vectors that emanate from the same location on the start domain, while a "channel" is 
defined by two parallel propagation vectors that emanate from two different locations on the start 
domain. After identifying the wedge or channel in which the destination point falls, the process 
computes (at 5415) the PLF-value at the destination point. The PLF-value at the destination point 

5 that is within a wedge equals (1) the PLF-value of Current.Drop's G PLF at the start domain vertex 
from which the wedge's propagation vectors emanate, plus (2) the distance (according to equation 
(5)) between this vertex and the destination point. On the other hand, the PLF-value at a destination 
point that is within a channel equals (1) the length of a line segment that is parallel to the two 
channel-defining vectors and that starts at the destination point and terminates at the start domain, 

10 plus (2) the PLF-value of the CurrentJDrop's G PLF at the point that the line segment terminates on 
the start domain. The line segment terminates on the start domain on a second line segment that is 
between the two knot locations from which the two channel-defining vectors emanate. When the 
start domain is a surface, the second line segment (1) is an edge on the boundary of the surface if the 
two knot locations are boundary vertices of the surface, and (2) is a line segment within the surface if 

15 the two channel-defining knot locations are within the surface. 

For instance, in Figure 55, the destination point 5520 falls within a channel defined by two 
propagation vectors 5530 and 5535. The distance between the destination point 5520 and line 5505 is 
the length D of a line segment 5550 that is parallel to vectors 5530 and 5535. Accordingly, the PLF- 
value at destination point 5520 equals the length D plus the PLF-value of the Current_Drop's G PLF 

20 at point 5545, which is the location that line segment 5550 intersects line 5505. If the destination 
point was point 5525 that is within the wedge defined by propagation vectors 5530 and 5540, the 
PLF-value at point 5525 would be (1) the PLF-value of the Current_Drop's G PLF at point 5515, 
plus (2) the distance (according to equation (5)) between points 5515 and 5525. 

In Figure 56, the destination point 5620 falls within a channel defined by two propagation 

25 vectors 5624 and 5626. The distance between the destination point 5620 and surface 5602 is the 



length D of a line segment 5630 that is parallel to vectors 5624 and 5626. Accordingly, the PLF- 
value at destination point 5620 equals the length D plus the PLF-value of the Current JDrop's G PLF 
at point 5662, which is the location that line segment 5630 intersects surface 5602. If the destination 
point was point 5618 that is within the wedge defined by propagation vectors 5622 and 5624, the 
5 PLF-value at point 5618 would be (1) the PLF-value of the CurrentJDrop's G PLF at point 5610, 
plus (2) the distance (according to equation (5)) between points 5610 and 5618. 
After 5415, the process 5400 terminates. 

(6th) Expansion from Surface to Surface 
A G PLF is propagated from a surface to another surface when, for example, the expansion is 
10 from one hole to another hole. In the embodiments described below, such an expansion would define 
a topologically stacked via, which is a topologically defined via that would start and end on two non- 
adjacent layers. A topologically defined stacked via does not always result in a geometrically stacked 
via. 

A viable expansion from one hole to another would require that the start and destination 
15 surfaces (i.e., the portions of the hole polygons to which a path can expand) have sufficient overlap. 
Some embodiments define the G PLF of such an expansion only over the region of the destination 
polygonal surface that overlaps the start polygonal surface {i.e., only at the intersection of the 
domains of the start and destination surfaces). In this region, the expansion's G PLF would be 
identical to the Current JDrop's G PLF for the corresponding region of the start surface. Of course, 
20 the expansion's G PLF might have new knots and might have modified edge and facet descriptions to 
account for the boundary of the overlapping regions. 

Figure 57 presents an example that illustrates an expansion from a start surface to a 
destination surface. Two spaces in layers 2 and 3 might contain the start surface, while two spaces in 
layers 3 and 4 might contain the destination surface. Figure 57 illustrates a start polygonal surface 
25 5705 and a destination polygonal surface 5710. This figure also illustrates a region 5715 that the 



polygons 5705 and 5710 overlap. It also illustrates the G PLF 5720 that is defined over the start 
surface 5705. This figure illustrates a projected view of the G PLF 5720 onto the x,y plane, in order 
to simplify the visual presentation of this PLF. 

In Figure 57, the overlap region 5715 is sufficiently large. Hence, the expansion is viable, 
5 and a G PLF needs to be computed for the expansion. Figure 57 illustrates that the expansion's G 
PLF 5770 is defined only over the overlapping portion 5715. This PLF 5770 is identical to the 
portion of the CurrentJDrop's G PLF 5720 that is defined across the overlapping region 5715. 
However, the PLF 5770 has new knots and modified edge and facet descriptions in order to account 
for the boundary of the overlap region. In particular, the expansion's G PLF 5770 will include two 

10 facets 5725 and 5730. Facets 5725 and 5730 have the same normal as the two facets 5740 and 5745 
of the Current_Drop G PLF. However, these facets have slightly different knots and edges. 
Specifically, facet 5725 includes knots 5750, 5754, 5756, 5764, and 5766, while facet 5730 includes 
knots 5756, 5758, 5760, 5762, and 5764. All these knots specify the boundaries of facets 5725 and 
5730, which are smaller than the facets 5740 and 5745 of the CurrentJDrop's PLF. Also, the edge 

15 5772 between facets 5725 and 5730 is shorter than the edge between facets 5740 and 5745. 

(ii) Penalty Cost 

Propagating the CurrentJDrop's G PLF to a potential destination particle specifies an initial G 
PLF for a potential expansion. If the potential expansion needs to be penalized, the process 3500 
adds to the initial G PLF one or more penalty costs associated with the potential expansion. Some 
20 embodiments penalize expansions to holes, expansions to overcongested edges and expansions to 
walls of other nets. Some embodiments also penalize expansions that shove routes of other nets 
along edges or within faces. 

(1st) Penalty for Expansion to overcongested edges 
Process 3500 allows a path to expand to the exits of overcongested edges. The route flow 
25 across an edge equals the width of the net routes crossing the edge plus the spacing between the 
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crossing net routes and between the net routes and the edge nodes. Figures 58A and 58B illustrate 
how to compute the flow across an edge after a potential expansion. Specifically, Figure 58A 
illustrates the center-lines of topological routes for two nets that were previously inserted across an 
edge 5800. The two routes have widths Wl and W2 and spacing SI and S3, respectively, towards 

5 their adjacent edge nodes. Also, the spacing between the two routes is defined as S2. Figure 58B 
illustrates how the edge 5800 would look like after a third route is inserted across it. The flow of this 
edge equals the sum of the following: (1) the minimum spacing SI between net 1 and its adjacent 
node, (2) the width Wl of net 1, (3) the minimum spacing S4 between nets 1 and 3, (4) the width W3 
of net 3, (5) the minimum spacing S5 between nets 3 and 2, (6) the width W2 of net 2, and (7) the 

10 minimum spacing S3 between net 3 and its adjacent node. 

When the route flow is equal to or less than the capacity of the edge, the expansion is allowed 
and legal. On the other hand, when this flow is larger than the edge capacity, the expansion is still 
allowed, but it is recorded as an illegal expansion. Some embodiments compute the capacity of each 
edge according to the processes described above by reference to Figures 19 and 20. 

15 Some embodiments pretabulate in a storage structure the widths for each net on each layer 

and the spacing requirements between each net and every other net, pin, or obstacle on each layer. 
Accordingly, in these embodiments, the path-generation process 3500 retrieves these values from the 
storage structure when it needs them (e.g., when it is computing edge capacities). 

An expansion to an edge that is overcongested before or after the expansion is assessed an 

20 extra penalty cost. This penalty is called a congestion penalty. In some embodiments, this congestion 
penalty is a pre-defined penalty constant. In other embodiments, this penalty is a constant that is 
derived based on the congestion of the edge. This constant is linearly or non-linearly (e.g., 
quadratically) proportional to the amount of over-congestion on the edge. 

(2nd) Expansion to holes 

25 Expansions to holes are also assessed a penalty cost. This penalty can be different for holes 



between different layer pairs if the cost introduced by a via between different layer pairs is different. 
A via introduces a variety of costs in a design, and the penalty can take into account all these costs. 
For instance, from a resistance or delay point of view, a via might cost fifty times more than the 
resistance or delay of a unit length of a wire. Accordingly, the penalty might be 50 (Le. 9 the via can 
be counted as a wire that is 50 units long). Alternatively, each via might reduce the manufacturing 
yield by some amount. This reduction can also be accounted for by equating it to a wirelength and 
adding it to the cost function. 

(3rd) Expansion to Walls of Other Nets 

As mentioned above by reference to Figure 37, the process 3500 allows the path for one net 
to expand to the wall of another net {i.e., to specify a drop for an expansion to the wall of another 
net). The process, however, adds a rip penalty to the cost G PLF) of an expansion to a wall of 
another net. This is because if such a drop is ever retrieved from the storage structure, the expansions 
that are defined about this drop will require piercing of the other net's route. In some embodiments, 
the rip penalty is larger than two via penalties. Such a penalty ensures that the path search only rips 
another net's route when it has reasonably exhausted all other viable expansions on the same layer or 
in the layers above or below. 

(4th) Expansion that Shoves Previously Defined Routes 

Some embodiments also penalize expansions that shove routes of other nets along edges or 
within faces. 

• Shoving Other Routes Along Edges 
A potential expansion to an edge might cause the other net routes on the edge to have to be 
shoved along that edge. Figure 59 illustrates a process 5900 that generates a line PLF that expresses 
the extra shoving cost of an expansion to an exit on an edge with one or more net routes crossing it. 
This process is described by reference to Figure 60, which illustrates an expansion from an exit 6002 

on a first edge 6004 to an exit 6006 on a second edge 6008. In Figure 60, a previously defined route 
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6010 of another net crosses the second edge 6008. This route also crosses edges 6004 and 6028. 
Joints 6018, 6030, and 6028 specify this route's crossing of edges 6004, 6008, and 6028. 

As shown in Figure 59, the process 5900 initially determines (at 5905) whether the route of 
at least one other net intersects the edge that contains the destination exit. If not, the process 
5 terminates. Otherwise, the process 5900 selects (at 5907) one of the routes that crosses the 
destination exit's edge. 

Next, at 5910, the process identifies the vectors to project from the selected route's joints on 
the edges that the shoved route crosses before and after the destination edge (/.<?., before and after the 
edge containing the expansion's destination exit). Some embodiments project vectors in each 

10 available interconnect-line direction that falls within a four-sided polygon that is formed by the 
selected route's two joints and by two points that define the exit's domain (Le. t maximum possible 
boundaries) along its edge. The boundary-defining points are defined based on the spacing and width 
constraints and obstacle constraints on the destination exit's edge. For instance, in Figure 60, points 
6016 and 6024 are the two boundary-defining points of the exit. The destination exit's boundaries 

15 cannot be defined beyond these two points without violating required spacing and width constraints 
between the net routes crossing the edge and between the nets and obstacles near the edge. 

The points 6016 and 6024 define a four-sided polygon along with joints 6018 and 6026. 
Within this polygon, two vectors 6012 and 6014 are projected from joint 6018 on edge 6004, and one 
vector 6020 is projected from joint 6026 on edge 6028. By projecting vectors from the joints 6018 

20 and 6026, the process 5900 assumes that these joints will remain fixed as the route's joint 6030 is 
shoved along the edge 6008. Figure 60 also illustrates a vector 6022 that emanates at a 135° angle 
from joint 6026. This vector is not projected as it does not fall within the polygon defined by points 
6016, 6018, 6024, and 6026 (/.<?., this vector would intersect the edge 6008 outside of the destination 
exit's domain). 

25 After identifying (at 5910) the vectors to project from the previous and next joints of the 



shoved route, the process 5900 identifies (at 5915) the intersection of the identified vectors and the 
destination edge. In Figure 60, projection vectors 6012, 6014, and 20 intersect the destination 
edge at points 6034, 6036, and 6038. 

At 5915, the process specifies a line PLF that is defined over the destination exit's domain 
5 and that expresses the cost of shoving the route along the destination exit's domain. The expansion to 
the destination exit necessarily requires the shoving of one or more routes on the edge when this 
exit's edge does not have any portion to which the path can expand without shoving the selected 
route. In this situation, the process specifies a knot in the line PLF for each intersection point 
identified at 5910. In this situation, the process also specifies two knots for the two points that define 

10 the exit's maximum possible boundaries along its edge. For each knot that the process specifies in 
this situation, the process computes a PLF-value that equals (1) the distance (according to equation 
(5)) between the knot's location and the location of the selected route's joint on the previous edge, 
plus (2) the distance (according to equation (5)) between the knot's location and the location of the 
selected route's joint on the next edge, minus (3) the current length (according to equation (5)) of the 

15 route between its previous and next joints. 

On the other hand, when the destination exit's edge has a segment to which the path can 
expand without shoving any route on the edge, the process specifies a knot at each identified 
intersection point that is outside the segment that can be expanded to without shoving other routes. 
In this situation, the process also specifies a knot at each of the two boundary points of this segment 

20 (i.e., the segment to which the path can expand without shoving any route) and sets the PLF-value of 
these knots to zero. The PLF-value of each of the other knots that the process specifies at this stage 
equals (1) the distance (according to equation (5)) between the knot's location and the location of the 
selected route's joint on the previous edge, plus (2) the distance (according to equation (5)) between 
the knot's location and the location of the selected route's joint on the next edge, minus (3) the 

25 current length (according to equation (5)) of the route between its previous and next joints. 
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In the example illustrated in Figure 60, the PLF-value of a knot specified at point 6038 is 
(1) the distance (according to equation (5)) between point 6038 and joint 6018, plus (2) the distance 
(according to equation (5)) between point 6038 and joint 6026, minus (3) the length (according to 
equation (5)) of the route between its joints 6018 and 6026. 
5 The line PLF that the process defines at 5915 expresses the additional wirelength due to 

shoving of the selected route along the edge. After 5915, the process determines (at 5920) whether it 
has generated such a line PLF for each route that crosses the destination exit's edge. If not, the 
process returns to 5907 to select another route and to repeat operations 5910-5920 for this newly 
selected route. When the process determines (at 5920) that it has generated a line PLF for each route 

10 that crosses the destination exit's edge, it adds (5925) the PLFs that it generated if it generated more 
than one PLF. After 5925, it then terminates. 

• Shoving Other Routes in a Face 
A potential expansion to a surface (e.g., a hole or a portion of a hole) within a face might 
cause the other net routes that intersect the face to be shoved within the face. Figure 61 illustrates a 

15 process 6100 that generates a surface PLF that expresses the extra shoving cost of an expansion to a 
surface that has one or more other net routes crossing its face. This process is described by reference 
to Figure 62, which illustrates an expansion from an exit 6202 on a first edge 6204 to a hole 6206 in 
a face 6208. In Figure 62, a previously defined route 6210 of another net crosses the face 6208. This 
route crosses the first edge 6204 and a second edge 6212 at two joints 6214 and 6216. 

20 As shown in Figure 61, the process 6100 initially determines (at 6105) whether the route of 

at least one other net intersects the face that contains the destination surface. If not, the process 
terminates. Otherwise, the process 6100 selects (at 6110) a route that intersects the destination 
surface's face. It then identifies (at 6115) two locations from which to project vectors on the edges of 
the destination surface's face. In some embodiments, these two locations are a distance S away from 

25 the selected route's joints along the edges of these joints. The distance S accounts for the via width 
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and the minimum spacing between the current net's route (i.e., the route of the net currently being 
routed) and the selected route. Figure 62 illustrates two points 6220 and 6222 that are a distance S 
away from joints 6214 and 6216 of the route 6210. 

Next, the process 6100 projects (at 6120) vectors from the locations identified at 6105. Some 

5 embodiments project vectors in each available interconnect-line direction that falls within the 
destination surface's face. Figure 62 illustrates (1) three vectors 6224, 6226, and 6228 that are 
projected from location 6220 on edge 6204, and (2) three vectors 6230, 6232, and 6234 that are 
projected from location 6222 on edge 6212. By projecting vectors from locations 6220 and 6222 that 
are identified from joints 6214 and 6216, the process 6100 assumes that the joints 6214 and 6216 

10 will remain fixed as the route is shoved within face 6208. 

After projecting vectors at 6120, the process 6100 identifies (at 6125) the intersection of the 
projected vectors (1) with the boundary of the destination polygonal surface and (2) with each other 
within the boundary of the destination polygonal surface. In Figure 62, the projection-vector 
intersections are at points 6240 and 6242. The projection vectors intersect the boundary of the 

15 destination surface at 6246-6256. 

At 6130, the process specifies a surface PLF that is defined over the destination surface and 
that expresses the cost of shoving the route in this surface. The process specifies a knot at each 
identified intersection point, unless the intersection point is within a portion of the destination 
surface to which the path can expand without shoving the selected route. The process sets the PLF- 

20 value of each knot that is specified at such an intersection point equal to (1) the distance (according 
to equation (5)) between the knot's location and the location of the selected route's joint on one edge 
of the destination surface's face, plus (2) the distance (according to equation (5)) between the knot's 
location and the location of the selected route's joint on the other edge of the destination surface's 
face, minus (3) the current length (according to equation (5)) of the route between the two joints. For 

25 instance, the PLF-value of a knot specified at point 6240 is the (1) the distance between point 6240 



and location 6220, plus (2) the distance between point 6240 and location 6222, minus (3) the length 
of the route between joints 6214 and 6216. Similarly, the PLF-value of a knot specified at 
intersection point 6246 equals the sum of the distance between this point and points 6220 and 6222, 
minus the length of the route between joints 6214 and 6216. 
5 The intersection of one or more projected vectors with each other or with the destination- 

surface boundary might fall within a portion of the destination-surface to which a path can expand 
without shoving the route. Of these intersection points, knots are specified only for those points that 
might define the boundary between the destination-surface portion that a path can expand to without 
shoving the selected route and the destination-surface portion that a path can only expand to by 
10 shoving the selected route. In general, the process 6100 specifies three or more knots to define the 
facet for the destination-surface portion that a path can expand to without shoving the selected route. 
This facet is a zero-shove cost facet, and hence the PLF-value of its knots will be zero. In Figure 62, 
the points 6220, 6254, 6260, and 6266 specify the portion of the destination surface to which a path 
can expand without shoving the route 6210. In some cases, a destination surface will not have any 
15 portion to which a path can expand without shoving the selected route. At 6125, the process also 
specifies knots at the remaining unexamined vertices of the destination surface. For instance, in the 
example illustrated in Figure 62, the process specifies knots at the vertices 6262 and 6264. 

After specifying all the knots, the process specifies (at 6130) the edges and facets of the 
surface PLF from the specified knots and the vectors that were used to define these knots. The 
20 process also computes and specifies (at 6130) the normal and z-intercept values for each specified 
facet. After 6130, the process determines (at 6135) whether it has generated a surface PLF for each 
route that crosses the destination-hole's face. If not, the process returns to 6110 to select another 
route and to repeat operations 6115-6135 for this newly selected route. When the process determines 
(at 6135) that it has generated a surface PLF for each route that crosses the destination surface's face, 
25 it adds (6140) the PLFs that it generated if it generated more than one PLF. After 6140, it then 
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terminates. 

(5) Adding PLF's 

The process 3500 adds point, line, and surface PLFs. Adding two PLF's at a point is trivial, 
as their values are simply added. A constant function can be added to a line or surface PLF by 
5 incrementing the PLF-values of all knots of the line or surface PLF by value of the constant function. 

(i) Adding Two Line PLF's 

Adding two convex line PLFs that are defined across the same domain can also be efficiently 
done by taking advantage of the piecewise linear nature of the functions. At any point Q, the sum of 
two convex line PLFs PLF1 and PLF2 that are defined over the same line is the sum of PLF-values 

10 of PLF1 and PLF2 at that point, as illustrated by the equation below. 

V = PLF1(Q) + PLF2(Q). 
Accordingly, in some embodiments, the process 3500 adds two line PLFs by performing the 
following two operations. First, it initializes to null a PLF3 that is to represent the sum of the two 
PLFs. Second, it specifies a knot in PLF3 at each unique knot location in either PLF1 or PL2. The 

15 value of each specified knot in PLF3 is the sum of the PLF-values in PLF1 and PLF2 at the specified 
knot's location. In some embodiments, the process 3500 examines the knots of PLF1 and PLF2 based 
on the order that they appear, i.e. the process starts from one end of the line over which the PLFs are 
defined and traverses to the other end. If a first line PLF (Fl) is defined over a larger portion of a line 
than a second line PLF (Fl), the first and second PLFs can be summed over their overlapping 

20 portion of their domains by (1) specifying a third line PLF (F3) that is equal to the first PLF but is 
only defined over the overlapping portion, and (2) performing the above described approach to 
obtain the sum of the second and third PLFs (F2 and F3). 

(ii) Adding Two Surface PLF's 

Figure 63 illustrates a specific process 6300 for adding two surface PLFs that have 
25 overlapping but potentially different domains. To add two surface PLF's, this process performs a line 
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sweep operation that identifies the intersection of the two polygonizations. This process is described 
by reference to Figure 64, which illustrates an example where two surface PLFs 6405 and 6410 are 
added to produce a third surface PLF 6415. In Figure 64, the PLFs are projected onto the x,y plane 
in order to simplify their visual presentation. 
5 As shown in Figure 63, the process 6300 initializes (at 6305) a new PLF ("result PLF") that 

will eventually represent the sum of the two PLFs. Next, at 6310, the process identifies the leftmost 
knot location in either PLF, and sets a variable CurX to the x-coordinate of the identified knot 
location. The process then selects (at 6315) a knot in either PLF that has an x-coordinate equal to 
CurX. There might be several knots with a particular x-coordinate in each PLF. To address this 
10 situation, some embodiments use the convention that when two or more knots in a PLF have the 
same x-coordinate, the process should select (at 6315) the knot that has the smallest y-coordinate and 
that has not yet been examined. At 6315, the process specifies a knot in the result PLF at the x,y 
location of the selected knot, if one such knot has not already been specified. At 6315, the process 
sets the PLF-value of the specified knot as the sum of the PLF-values of the two PLFs at the 
15 specified knot's x,y coordinates. 

As further described below, each time the process 6300 encounters a knot that "opens" an 
edge (i.e., encounters a new edge) in its left-to-right sweep, it adds the edge to a list of open facets 
for the selected knot's PLF. Some embodiments use the convention that a knot opens an edge if 
(1) its x-coordinate is smaller than the x-coordinate of the other knot of its edge, or (2) its y- 
20 coordinate is smaller than the y-coordinate of the other knot of its edge when both edge knots have 
the same x-coordinate. 

Before adding a new edge to the open list, the process 6300 determines (at 6320) whether the 
selected knot closes any edge in the selected knot's PLF. In some embodiments, a knot closes an edge 
if (1) its x-coordinate is larger than the x-coordinate of the other knot of its edge, or (2) its y- 
25 coordinate is larger than the y-coordinate of the other knot of its edge when both edge knots have the 
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same x-coordinate. If the process determines (at 6320) that the selected knot does not close any 
edges, the process transitions to 6345, which is further described below. If the selected knot does 
close an edge, the process selects (at 6325) an edge that the selected knot closes. It then performs (at 
6330) an intersection operation between the closed edge and each edge on the open list of the other 
5 PLF. 

Figure 65 illustrates the intersection process 6500. This process initially compares (at 6505) 
the closed edge with each edge on the open list of the other PLF in order to identify each intersection 
of the closed edge with another edge. It then determines (at 6510) whether it identified at 6505 the 
intersection at a point of the closed edge with any of the open edges of the other PLF. 

10 If the process determines (at 6510) that the closed edge does not intersect any open edge, the 

process adds (at 6515) the closed edge to the result PLF and then terminates. Otherwise, the process 
specifies (at 6520) a knot at each point that the closed edge intersects one of the open edges of the 
other PLF. The value of each specified knot is the sum of the two PLFs at the location of the 
specified knot. For instance, in Figure 64, when the edge between knots 6458 and 6460 in the first 

15 PLF 6405 closes, the process 6500 identifies this edge's intersection with the edge between knots 
6462 and 6464 in the second PLF 6410. At the intersection of these two edges, the process 6500 
identifies a knot 6455 in the third PLF 6415. The value of the knot 6455 equals the sum of the values 
of PLFs 6405 and 6410 at this knot's x, y coordinates. 

After specifying one or more knots at 6520, the process 6500 defines (at 6525) two or more 

20 col linear edges to replace the closed edge, and adds the collinear edges to the result PLF. In Figure 
64, the process identifies (at 6525) two new edges, one between knots 6450 and 6455, and one 
between knots 6455 and 6454. It then adds these two edges to the definition of the third PLF 6415, in 
place of adding the closed edge between knots 6450 and 6454. 

For each particular edge on the open list of the other PLF that intersects the closed edge, the 

25 intersection process 6500 (at 6530) (1) defines two new edges, (2) adds one of the two edges to the 



result PLF, and (3) replaces the particular edge with the other of the two edges. The process uses the 
knot specified at 6520 at the intersection of the closed edge and an open edge to define the two new 
edges for the open edge. For instance, in Figure 64, the process (at 6530) would define two edges for 
the open edge between knots 6462 and 6464. One of the two edges is between knots 6462 and 6455, 
5 and one edge is between knots 6455 and 6464. The process would add the edge between knots 6455 
and 6464 to the result PLF. It would also replace the edge between knots 6462 and 6464 on the open 
list of PLF 6410 with the edge between knots 6455 and 6464. After 6530, the process terminates. 

After the process 6300 performs the intersection process at 6330, the process 6300 removes 
(at 6335) the selected closed edge from its open list. It then determines (at 6340) whether it has 

10 examined all the edges that the selected knot (i.e., the knot selected at 6315) closes. If not, the 
process transitions back to 6330 to select another closed edge. 

Otherwise, the process adds (at 6345) any edge that the selected knot opens to the list of 
"open" edges for the selected knot's PLF. After 6345, the process determines (at 6350) whether it has 
examined each knot in each PLF with an x-coordinate CurX. If not, the process transitions back to 

15 6315 to select another knot with the x-coordinate CurX. If so, the process transitions to 6355, where 
it determines whether it has examined all the unique x -coordinates of knots in either PLF. If not, the 
process sets (at 6360) CurX to the next leftmost x-coordinate of any knot in either PLF. 

When the process determines (at 6355) that it has examined all the unique x-coordinates of 
knot locations in both PLFs, it examines (at 6365) the edges added to the result PLF and identifies 

20 facets formed by related sets of edges. At 6365, the process also identifies the normal and z-intercept 
of each facet based on the knot values of the facet. When one (Fl) of the PLFs being summed is 
defined over a larger portion of a surface than the other PLF (F2), the process 6300 can limit (at 
6365) the surface PLFs definition that represents the sum of these two PLFs to cover only the 
overlapping domain of the two PLFs. Alternatively, before 6305, the process 6300 could generate a 

25 third PLF (F3) that is equal to the first PLF but is only defined over the domain of the second PLF 
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(i.e., only defined over the portion of the surface that both PLFs Fl and F2 are defined). 

Figure 64 illustrates the result of the addition of the two PLFs 6405 and 6410. In this figure, 
the first PLF 6405 has five boundary knots 6430-6438 that are at the same location as five boundary 
knots 6440-6448 of PLF 6410. For these five pairs of knots, the process identifies five knots 6420- 
5 6428 in the third PLF 6415. In the third PLF 6415, the process also identifies four knots for the other 
four unique knots 6458-6462 in the two PLFs 6405 and 6410. As described above, it also specifies a 
knot 6455 at the intersection of two non-boundary edges of the two PLFs. In turn, these knots 
specify 13 edges that define four facets 6466-6472. 
(6) Filtering PLF's 

10 After identifying the G PLF of a potential expansion, the path-generation process 3500 filters 

(at 3534) the expansion's G PLF with the filter PLF of the destination particle. This filtering 
operation also sets the destination particle's filter PLF equal to the minimum of its original filter PLF 
and the expansion's G PLF. Filtering a point PLF for an expansion to a point is trivial. The expansion 
is discarded when its PLF-value is greater than the filter function value at that point. Otherwise, it is 
15 kept, and the point's filter function is set to the expansion's PLF-value. 

(1st) Filtering two line PLF's and Identifying their minimum 
Figure 66 illustrates a process 6600 that performs the filtering and minimum operations for 
an expansion to a line (i.e., an expansion to an exit or a wall). Specifically, this function filters a 
filtered PLF (PLFl) by a filter PLF (PLF2), which may be non-convex. In the embodiments 
20 described below, the filter PLF2 is defined across a topological particle, while the filtered PLFl 
might be defined across only a portion of the topological particle. Accordingly, the process 6600 
performs its filtering and minimum operations for only the portion of the particle over which the 
filtered PLFl is defined. Other embodiments might perform the filtering over the entire particle. In 
the embodiments described below, the filtering operation discards the portions of the filtered PLF 
25 that are larger than the corresponding portions of the filter PLF. Other embodiments, however, might 
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not discard these portions, but rather might set the PLF-values in these portions to infinite. These 
embodiments then would not need to define several drops when several pieces of the filtered PLF 
remain after the filter. 

As shown in Figure 66, the process 6600 initially sets (at 6605) a state variable to 
"Unknown" and sets two PLFs PLF3 and PLF4 to null. The process uses PLF3 to represent the 
minimum of the filtered and filter PLF, while it uses PLF4 to keep track of the remaining portions of 
the filtered PLF. As the minimum, PLF3 will represent at each offset Q along the line the minimum 
of PLF1(Q) and PLF2(Q). PLF3 might end up being a non-convex PLF. As PLF1 and PLF2 are, 
PLF3 is defined over the line and will be represented by a sequence of knots that is sorted based on 
the offset values of the knots. 

Next, at 6610, the process identifies the unique knot locations in both PLF1 and PLF2 and 
generates a set of unique knot locations that are sorted based on the offset values of the knots in this 
set. The process then selects (at 6615) the first location in the sorted set. The process then determines 
(at 6620) whether the filtered PLF1 is greater than the filtered PLF2 at the location selected at 6615. 

If not, the process determines (at 6625) whether the state is PLFl_Lower or unknown. When 
the process determines at 6625 that the state is PLFl_Lower or is unknown, the process (at 6635) 
specifies a knot at selected location, sets its value to the value of PLF1, and adds the specified knot 
to PLF3 and PLF4. At 6635, the process also sets the state to PLFlJLower if the state was previously 
unknown. In this situation, the process performs these operations because (1) if the state is unknown, 
this knot is at the first knot location and it is from the filtered PLF1, and (2) if the state is 
PLFl_Lower, this knot is from the same PLF as the last knot added to PLF3 (i.e., PLF1 and PLF2 
have not crossed between the last knot added to PLF3 and this knot, and during this internal PLF1 
has been the smaller PLF). From 6635, the process transitions to 6655, which will be further 
described below. 

On the other hand, when the process determines at 6625 that the state is PLF2_Lower, the 
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knot at selected location is from a different PLF than the last knot added to PLF3. Hence, PLF1 and 
PLF2 have crossed between the last knot added to PLF3 and this knot. During this internal, PLF2 
started smaller and ended larger. Accordingly, the process transitions to 6630 to perform several 
operations to account for this crossing. At 6630, the process initially sets the state to be 
PLFl_Lower. It then identifies the offset value Q along the edge where PLF1 and PLF2 intersected 
between the last knot added to PLF3 and the knot at the location selected at 6615. This intersection 
can be easily computed as it is an intersection between two line segments representing PLF1 and 
PLF2 between the last knot added to PLF3 and the knot at selected location. The process then 
specifies a knot at intersection point Q and adds this knot to PLF3 and PLF4. The PLF-value of this 
added knot is the PLF-value of either PLF1 or PLF2 at the offset-value of the intersection point Q. If 
the PLF-value of a knot at the selected knot location is smaller than the other PLFs value at this 
location, the process specifies a knot at the selected knot location, sets its value to the value of PLF1, 
and adds the specified knot to the description of PLF3 and PLF4. From 6630, the process then 
transitions to 6655. 

If the process determines (at 6620) that the filtered PLF1 is greater than the filter PLF2 at the 
location selected at 6615, the process determines (at 6640) whether the state is PLF2_Lower or 
unknown. When the process determines at 6640 that the state is PLF2_Lower or unknown, the 
process (at 6645) specifies a knot at selected location, sets its value to the value of PLF2 at this 
location, and adds the specified knot to PLF3. At 6645, the process also sets the state to PLF2_Lower 
if the state was previously unknown. The process performs these operations at 6645 because (1) if 
the state is unknown, this knot is the first knot selected, and it is from the filter PLF2, and (2) if the 
state is PLF2_Lower, this knot is from the same PLF as the last knot added to PLF3 (i.e., the two 
PLFs have not crossed between the last knot added to PLF3 and this knot, and during this interval 
PLF2 has been the smaller PLF). From 6645, the process transitions to 6655. 

On the other hand, when the process determines at 6640 that the state is PLFlJLower, the 
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knot at selected location is from a different PLF than the last knot added to PLF3. Hence, the two 
PLFs have crossed between the last knot added to PLF3 and the knot selected at 6615. During this 
interval, PLF1 started smaller and ended larger. Accordingly, the process transitions to 6650 to 
perform several operations to account for this crossing. At 6650, the process initially sets the state to 
be PLF2_Lower. It then identifies the offset value Q along the edge where PLF1 and PLF2 
intersected between the last knot added to PLF3 and the knot selected at 6615. This intersection can 
be easily computed, as it is an intersection of two line segments representing PLF1 and PLF2 
between the last knot added to PLF3 and the knot at selected location. The process then specifies a 
knot at intersection point Q and adds this knot to PLF3 and PLF4. The PLF-value of this added knot 
is the PLF-value of either PLF1 or PLF2 at the offset-value of the intersection point Q. If the PLF- 
value of a knot at the selected knot location is smaller than the other PLFs value at this location, the 
process specifies a knot at the selected location, sets its value to the value of PLF2 at this location, 
and adds the specified knot to the description of PLF3 and PLF4. The process then specifies a new 
PLF equal to PLF4 and then reinitializes PLF4. The new PLF represents a remaining portion of the 
filtered PLF. From 6650, the process transitions to 6655. 

At 6655, the process determines whether it has examined all the unique knot locations in the 
sorted set produced at 6610. If not, the process returns to 6615 to select the next unique knot location 
in the sorted set and to perform the subsequent operations for this newly selected knot location. On 
the other hand, when the process determines (at 6655) that it has examined all the unique knot 
locations in the combined set, the process uses (at 6660) PLF3 to define the new filter function. As 
mentioned above, the filtered function PLF1 might not be defined over the entire domain over which 
the filter function PLF2 is defined. Because of this, the process 6600 performs its minimum 
operations for only the portion of the filter function's particle over which the filtered PLF1 is defined. 

When the filtered function PLF1 and the original filter function PLF2 have the same domain, 
the process sets (at 6660) the particle's filter function to PLF3. On the other hand, when the filtered 
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PLF1 is defined over a smaller portion of the filter function's particle, the process replaces the 
portion of the filter function (PLF2) that is defined over the filtered function's (PLFl's) domain with 
PLF3. The process leaves the portion of PLF2 outside of PLFl's domain intact. A boundary knot of 
PLF3 will be removed if it is an original knot of the filter PLF2 or if it is on a PLF segment in PLF3 

5 that has the same slope as the original filter PLF2 at the boundary knot. At 6660, the remaining 
portion or portions of the filtered PLF1 are the PLF's that the process 6600 defines at 6650. 

As mentioned above, the path-generation process 3500 filters (at 3512) the G PLF of a drop 
that it retrieves from the priority queue with the filter function of the drop's particle. For this filtering 
operation, the path-generation process uses a process that is identical to the process 6600 with the 

10 exception that unlike process 6600, the filtering process at 3512 does not maintain PLF3 as it does 
not need to identify the minimum of the filtered and filter PLF's at 3512. 

(2nd) Filtering two Surface PLF's 
Figure 67 illustrates a process 6700 for filtering two surface PLF's that have overlapping but 
potentially different domains. To filter two such surface PLF's, the process 6700 performs a plane 

15 sweep operation that identifies the intersection of the two polygonizations. This process starts by 
initializing (at 6705) a Min PLF that it will use to record the minimum of the two surface PLF's that 
it is filtering. 

It then identifies (at 6710) the leftmost unique knot location {i.e., the unique knot location 
with the smallest x-coordinate) in either PLF, and sets a variable CurX to the x-coordinate of the 

20 identified knot location. The process then selects (at 6715) a knot in either PLF that has an x- 
coordinate equal to CurX. There might be several knots with a particular x-coordinate in each PLF. 
To address this situation, some embodiments use the convention that when two or more knots in a 
PLF have the same x-coordinate, the process should select (at 6715) the knot that has the smallest y- 
coordinate and that has not yet been examined. 

25 As further described below, each time the process 6700 encounters a knot that "opens" a facet 



(i.e., encounters a new facet) in its left-to-right sweep, it adds the facet to a list of open facets for the 
selected knot's PLF. To address the case where the two knots of a facet have the same x-axis 
coordinate, some embodiments use the convention that a knot opens a facet if (1) it is the one with 
the smallest x-coordinate when all the facet's knots have different x-coordinates, or (2) it is the one 
5 with the smallest x-coordinate and the smallest y-coordinate when two facet's knots have the smallest 
x-coordinate of all knots of the facet. 

Before adding a new facet to the open list, the process 6700 determines (at 6720) whether the 
selected knot closes a facet in the selected knot's PLF. In some embodiments, a knot closes a facet if 
(l)the knot is the one with the largest x-coordinate when all the facet's knots have different x- 

10 coordinates, or (2) it is the one with the largest x-coordinate and the largest y-coordinate when two 
facet's knots have the largest x-coordinate of all knots of the facet. If the process determines that the 
selected knot does not close any facet, the process transitions to 6745, which is further described 
below. However, if the selected knot closes at least one facet, the process selects (at 6725) a facet 
that the selected knot closes. It then performs (at 6730) a min operation between the closed facet and 

15 each facet on the open list of the other PLF. 

Figure 68A illustrates the process 6800 for identifying the minimum of two facets Fl and F2 
from two different PLFs. This process initially identifies (at 6805) a polygon that represents the 
intersection of the domains of the two facets. Next, at each vertex of the polygon, the process (at 
6810) identifies and records the PLF that provides the minimum PLF- value. In case both PLF's 

20 provide the same value at a vertex, the process records the identity of the filtered PLF for the vertex 
at 6810. At 6810, the process also specifies a knot at each vertex of the intersection polygon. The 
PLF-value of a specified knot for a vertex is the value of the recorded PLF at the vertex. 

Next, the process determines (at 6815) whether one PLF provided all the lowest values at all 
the vertices of the polygon identified at 6805. If not, the two facets intersect and the process performs 

25 6820-6830 to account for this intersection. These operations will be described by reference to 
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Figure 68B, which illustrates two facets 6852 and 6854 that intersect. At 6820, the process 6800 
identifies an edge along which the two PLFs intersect, and identifies or specifies two knots that 
represent this edge. Figure 68B illustrates the identified edge 6856 between the two facets 6852 and 
6854, and the two knots 6858 and 6860 that define this edge. 
5 At 6825, the process then (1) uses the knots identified at 6820 to specify two new facets F3 

and F4 that represent the minimum of the intersected facets Fl and F2, and (2) adds these new facets 
F3 and F4 to the Min PLF. The two new facets will include the two knots identified at 6820 and one 
or more knots that were specified (at 6810) at the intersection polygon vertices. Each new facet will 
have edges between its knots. 

10 Figure 68B illustrates the two new facets 6862 and 6864. Facet 6862 is the portion of facet 

6854 that is below than the corresponding portion of the facet 6852. Similarly, facet 6864 is the 
portion of facet 6852 that is below the corresponding portion of the facet 6854. The facet 6862 is 
represented by edges between knots 6858, 6860, 6866, and 6868, while the facet 6864 is represented 
by edges between knots 6858, 6860, 6870, and 6872. The process can obtain the normal and z- 

15 intercept attributes of each specified facet 6862 and 6864 from the open or closed facet 6854 and 
6852 that corresponds to the specified facet. Alternatively, the process can compute the normal and 
z-intercept of each facet that it specifies at 6825 from its knots. 

Next, at 6830, the process adds the facet that it specified at 6825 for the smaller portion of the 
filtered PLF to a list that stores all remaining facets of the filtered PLF. If the facet 6854 in 

20 Figure 68B is from the filtered PLF, the process would add (at 6830) the facet 6862 to the list of the 
remaining facets of the filtered PLF. 

If the process 6800 determines (at 6815) that one PLF provided the lowest values at all 
vertices of the polygon identified at 6805, the process specifies (at 6835) a new facet that is defined 
over the polygon and that is identical to the PLF over this polygon. It sets the values at these knots 

25 and edges from the smaller PLF, and it sets the specified facet's normal and z-intercept identical to 
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the facet of the smaller PDF. The process then writes (at 6840) the facet specified at 6835 and its 
edges and knots into the Min PLF. If the smaller PLF is the filtered PLF, the process also adds (at 
6845) this information (i.e., the facet specified at 6835 and its edges and knots) to the facet list that 
stores the remaining facets of the filtered PLF. After 6830 and 6845, the process 6800 terminates. 

After the process 6700 performs (at 6735) the min process 6730 for each combination of the 
closed facet and the open facets on the other PLFs list, the process 6700 removes (at 6735) the 
selected closed facet from its open list. It then determines (at 6740) whether it has examined all the 
facets that the selected knot (Le. 9 the knot selected at 6715) closes. If not, the process transitions back 
to 6725 to select another closed facet. Otherwise, the process adds (at 6745) any facet that the 
selected knot opens to the list of "open" facets for the selected knot's PLF. The process then 
determines (at 6750) whether it has examined each knot in each PLF with an x-coordinate CurX. If 
not, the process transitions back to 6715 to select another knot with the x-coordinate CurX. If so, the 
process transitions to 6755, where it determines whether it has examined all the unique x-coordinates 
of knots in either PLF. If not, the process sets (at 6760) CurX to the next leftmost x-coordinate of 
any knot in either PLF. 

When the process determines (at 6755) that it has examined all the unique x-coordinates of 
knot locations in either PLF, it examines (6765) the facets added to the filtered list, and identifies 
connected sets of facet. Each connected sets of facets represents a remaining portion of the filtered 
PLF. It associated each connected sets of facet so that the set represents one remaining portion of the 
filtered PLF. It also removes any unnecessary edges and knots (if any) from each associated set and 
modifies the definition of one or more facets of the set accordingly. 

At 6765, the process also defines the filter PLF as the Min PLF. This Min PLF might include 
unnecessary edges and knots (e.g., at the boundaries of the facets that the process 6800 intersected). 
If so, before setting the filter PLF to the Min PLF, the process removes these unnecessary edges and 
knots at 6765 and modifies the definition of one or more facets of the Min PLF accordingly. 
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As mentioned above, the filter PLF might be defined over a larger domain than the filtered 
PLF. As an alternative to the above-described approach, some embodiment might perform the 
process 6700 for only the portion of the surface over which the filtered PLF1 is defined. In these 
embodiments, the filter PLF would be defined over a larger domain than the Min PLF. Accordingly, 
for such cases, the process 6700 only sets the portion of the filter PLF that is defined over the same 
domain as the filtered PLF to the Min PLF. The process retains the portions of the filter PLF that 
falls outside of the domain of the filtered PLF. To account for the boundary between retained and 
replaced portions of the filter PLF, the process might need to add knots and edges to the filter PLFs 
description and might need to modify the filter PLFs description of facets at the boundaries. 

As mentioned above, the path-generation process 3500 filters (at 3512) the G PLF of a drop 
that it retrieves from the priority queue with the filter function of the drop's particle. For this filtering 
operation, the path-generation process uses a process that is identical to the process 6700 with the 
exception that unlike process 6700, the filtering process at 3512 does not maintain Min as it does not 
need to identify the minimum of the filtered and filter PLFs at 3512. 
(7) Computing the H PLF 

In some embodiments, the process 3500 specifies the H PLF of a particle as the lower-bound 
distance between the particle and a bounding octagon that encompasses the set of target particles for 
a path search. Some embodiments identify such a lower-bound distance by (1) identifying an axis- 
aligned rectangular box (i.e. f a box that is aligned with the x-y coordinate axes of the layout) and a 
45° rotated rectangular box (i.e., a box that is aligned with an s-t coordinate axes that is rotated by 
45° with respect to the layout's x-y axes) that enclose the set of target particles, (2) identifying two 
PLFs that express the distance between the particle and each of these boxes, and then (3) defining 
the H PLF as the maximum of the two identified PLFs. 
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(i) Identifying the Bounding Boxes 

In some embodiments, the process identifies an axis-aligned bounding box by (1) identifying 
the minimum and maximum x- and y-coordinates (X M in, X M ax, Y M in, and Y M ax) of the set of points 
in the design layouts coordinate space, and (2) specifying the four vertices of the box as (X M in, 
Ymin), (Xmin, Ymax), (Xmax, Y M in), and (X M ax, Ymax). In some embodiments, the process 3500 
identifies the minimum and maximum x- and y-coordinates of the set of points by examining the x- 
and y-coordinates of vertex points in the set of target particles. Each standalone point in this set is a 
vertex point. Also, each point in the target set that is a vertex of a polygonal shape in the set is a 
vertex point of the set. 

The rotated bounding box is defined with respect to a coordinate system that is rotated by 45° 
with respect to the layout's coordinate system. The rotated coordinate system includes an s-axis that 
is at 45° with respect to the x-axis of the layout's coordinate space, and a t-axis that is at 135° with 
respect to the x-axis of the layout's coordinate space. To identify the 45° rotated bounding box, the 
process 3500 in some embodiments first maps the coordinates of the vertex points of the set of target 
particles from the layout's coordinate system to the 45° rotated coordinate system. The process 
performs this mapping by using the equation below. 

s i = *^, and 

T . X,+Y, 
V2 ' 

In this equation, (1) Xj and Yi are the coordinates in the layout's x- and y-axes for the i-th vertex 
point in the target set, and (2) Sj and Tj are the coordinates in the rotated s- and t-axes for the i-th 
vertex point in the target set. If the rotated bounding box is rotated counterclockwise at an arbitrary 
angle 6 with respect to the layout's coordinate system, the following equation could be used to 
perform the mapping. 
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After performing the mapping, the process 3500 then identifies a 45° bounding box that is 
defined with respect to the 45° rotated coordinate system and that surrounds the set of routed points 
in the s- and t-axes. In some embodiments, the process (1) identifies the minimum and maximum s- 
5 and t-coordinates (S M in, S M ax, T M in, and T M ax) of each vertex point in the set of routed points, and 
(2) specifies the four vertices of the rotated box as (Smin, Tmin), (Smin, Tmax), (Smax» Tmin)> and 
(Smax, Tmax)- 

(ii) Computing the PLF's 

In some embodiments, the process 3500 computes a PLF that represents the distance between 
10 a particle and a bounding box in an analogous manner to the propagating of a zero-cost function that 
is defined over the bounding box to the particle. Figure 69 illustrates this computation for some 
embodiments that use the wiring model of Figure 1. Specifically, this figure illustrates an axis- 
aligned bounding box 6905 and a 45° rotated bounding box 6910. Three vectors are projected from 
each vertex of each bounding box. Figure 69 illustrates the three vectors that are projected from each 
15 vertex of the boxes 6905 and 6910. Some embodiments do not project vectors in all interconnect 
directions. Instead, they project only the propagation vectors that will intersect the destination 
particle. These propagation vectors are the vectors that fall within a triangle defined by the bounding 
box's vertex point and the leftmost and rightmost vertices of the destination particle. 

If the topological particle is a point, the H PLF is a single value that is defined for the point. 
20 This value is the distance between the point and the bounding box. This distance is computed by 
(1) identifying the wedge or channel that contains the point, and (2) computing the distance 
according to the wedge and channel distance rules described above. For instance, in Figure 69, the 
H PLF-value of the point 6915, which falls within a wedge defined by vectors 6920 and 6922 that 
emanate from vertex 6924, is the distance between the point 6915 and the vertex 6924. On the other 



hand, if the particle is point 6936 that falls within the channel defined by vectors 6926 and 6928 that 
emanate from the rotated box 6910, H PLF-value is the length D of a line 6938 that is parallel to 
vectors 6926 and 6928 and that goes from the point 6936 to the edge 6932. 

If the topological particle is a line, the H PLF is a line PLF that is defined over the line. This 
line PLF will have knots at the intersections of the projected vectors and the line. The H PLF-value 
of each knot is simply the straightline distance between the intersection point and the vertex from 
which the intersecting vector projected. This straightline distance is along the direction of the 
projected intersecting vector. Also, a knot is specified at each unexamined endpoint of the line, and 
the value of such a knot depends on whether the endpoint falls within a channel or a wedge. 

For instance, when the particle is line 6940, two knots are specified at the intersections 6942 
and 6944 of the projection vectors 6946 and 6948 with the line. The PLF-value of the knot at 
location 6942 is the distance Dl to the vertex 6950 along the direction of vector 6946, while the 
PLF-value of the knot at location 6944 is the distance D2 to the vertex 6950 along the direction of 
vector 6948. Two knots are also specified at the endpoints 6952 and 6954 of the line 6950. The 
endpoint at 6952 falls within a wedge. Hence, the PLF-value of the knot specified at 6952 is based 
on the distance between point 6952 and vertex 6950. The endpoint at 6954, however, falls in a 
channel, and therefore the PLF-value of the knot defined at this point is the length D3 of the 135° 
line segment 6956 between the point 6954 and edge 6958. 

If the topological particle is a surface, the H PLF is a surface PLF that is defined over the 
surface. This surface PLF will have knots at the intersections of the projected vectors and the 
boundary of the surface. It will also have knots at each unexamined vertex of the surface. It will 
further have edges between the knots, and one or more facets defined by the edges. A normal and z- 
intercept will also be specified for each facet. The values of each facet's normal and z-intercept will 
be determined by the PLF-values of the knots that define each facet. 



101 



The PLF-values of each knot that is defined at an intersection of a projected vector is simply 
the straightline distance between the intersection point and the vertex from which the intersecting 
vector projected. This straightline distance is along the direction of the projected intersecting vector. 
Also, the PLF-value of each knot that is specified at an unexamined vertex of the surface depends on 
5 whether the vertex falls within a channel or a wedge. 

For instance, when the particle is a surface 6960 and the bounding box is box 6905, four 
knots are specified at the intersection 6962, 6964, 6966, and 6968 of the projection vectors 6970 and 
6970 with the surface. The PLF-value of each of these knots is the distance between the knot's 
location and the vertex that was used to define the knot. For instance, the PLF-value of the knot at 

10 location 6962 is the distance D5 between the vertex 6974 and point 6964 along the direction of 
vector 6970. In this example, five other knots are specified at the unexamined vertices 6976, 6978, 
6980, 6982, and 6984 of the surface. Of these vertices, only vertex 6980 falls within a channel. 
Accordingly, the PLF-value of the knot specified at the vertex 6980 is the length of the line segment 
6986 that is parallel to vectors 6970 and 6972. The PLF-value of the knot defined at each of the other 

15 four unexamined vertices is the distance between the knot's location and the bounding box vertex 
that projects the vector that defines the wedge containing the knot location. For instance, the PLF- 
value of the knot at 6982 is the distance between vertex 6982 and the vertex 6924 of the box 6905. 

The H PLF that is defined over the surface 6960 will have eleven edges. Nine of these edges 
will be defined about the boundary of the surface. Two of the edges 6988 and 6990 will be defined 
20 across the surface. These two edges in conjunction with the other nine boundary edges define three 
facets. The normal and z-intercept of each facet can be derived from the PLF-values of the knots that 
define the facet, by using standard plane-defining techniques. 

(Hi) Identifying the Maximum 

After defining the H PLFs for the two bounding boxes, some embodiments identify the 

25 maximum of the two PLFs. Obtaining the maximum of two PLFs is analogous to obtaining the 
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minimum of two PLFs, except that the portions that are discarded are the portions with the smaller 
PLF-values. Accordingly, the processes for computing the maximum of two line PLFs or two 
surface PLFs are analogous to the above-described processes for computing the minimum of such 
PLFs, except that the processes for computing the maximum retain the line segments or surface 
pieces that have the larger PLF-values. 

(iv) Different Computation of the H PLF 
Some embodiments specify the H PLF of a particle slightly differently. These embodiments 
first identify a bounding polygon by identifying and intersecting two bounding rectangular boxes that 
enclose the set of target particles. In some of these embodiments, one box is aligned with the x-y 
coordinate axes of the layout while the other box is aligned with an s-t coordinate axes that is 45° 
rotated with respect to the layouts x-y axes. Figure 70 illustrates an example of an octilinear 
bounding polygon 7025 that is identified by intersecting an axis-aligned box 7030 and a 45° rotated 
box 7035. 

After identifying such an octilinear bounding polygon, these embodiments identify a H PLF 
that expresses the distance between the particle and the octilinear bounding polygon. The 
identification of this PLF is analogous to propagating a zero-cost function that is defined over the 
octilinear bounding polygon to the particle. Specifically, one or more vectors are projected from each 
vertex of the octilinear bounding polygon. In Figure 70, four vectors are projected from vertex 7020, 
three vectors are projected from each of the vertices 7010 and 7015, and two vectors are projected 
from vertex 7005. These vectors are identified based on the following approach. For each vertex, two 
directions are identified that are perpendicular to the edges incident upon the vertex and that point 
away from the polygon, that is, they points left from an incoming clockwise edge, or right from an 
incoming counter-clockwise edge. Vectors are then projected from each vertex (1) in the two 
directions identified for the vertex, and (2) in any other direction that falls between the two identified 
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directions. Some embodiments do not project vectors in all interconnect directions. Instead, they 
project only the propagation vectors that will intersect the destination particle. These propagation 
vectors are the vectors that fall within a triangle defined by the bounding polygon's vertex point and 
the leftmost and rightmost vertices of the destination particle. 

After projecting the vectors from the vertices of the octilinear bounding polygon, the process 
then identifies a point PLF, a line PLF, or a surface PLF in the exact same manner as described 
above for the case of the axis-aligned and 45°-rotated bounding boxes. In other words, the same 
wedge and channel approach that was described above for the bounding boxes, can be used to 
identify the PLFs for the octilinear bounding polygons. One of ordinary skill will realize that the 
above-described H computation techniques (e.g., the wedge and channel distance computation 
technique) can be used with any type of bounding polygon, such as a convex hull or approximate 
convex hull. 

(8) Specifying Topological Route 
Figure 71 illustrates a process 7100 that the process 3500 performs (at 3542) to specify a 
topological path after it identifies a series of drops that connect the source and target sets of a path 
search. By the time the process 3500 reaches 3542, the Current.Drop is a drop that is on a target 
particle. Starting at the Current.Drop on the target, the process 7100 back traces the sequence of 
drops that reached the target and generates an ordered list of topological particles that define the 
topological path. 

The process 7100 initially defines (at 7105) a Previous_Particle to be null. At 7110, the 
process (1) examines the Current.Drop's particle, (2) if necessary, specifies one or more topological 
particles, and then (3) defines a new Previous.Particle. The CurrentJDrop's particle might be a node, 
an exit, a hole, or a wall. 

When the CurrentJDrop's particle is a node that is a target particle, the process does not 
specify any additional topological particles at 7110. When the Cun-entJDrop's particle is a node that 
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is not a target particle, the process (at 7110) specifies a wall that connects the node and the 
Previous_ParticIe, and then defines the node as the Previous_Particle. 

When the CurrentJDrop is an exit, the process creates two new exits and a joint between the 
two new exits at 71 10. In the edge of the CurrentJDrop's exit, the process replaces this exit with the 
two exits and the joint that it specifies at 71 10. If the CurrentJDrop's exit was not a target particle, 
the process also (at 7110) specifies a wall that connects the specified joint and the Previous_Particle 
and then defines the joint as the Previous_Particle. 

When the CurrentJDrop's particle is a hole, the process 7100 creates (at 7110) one via node 
on each layer traversed by the hole. It also specifies a wall that connects the Previous_Particle and 
the via node that is on the same layer as the Previous_Particle. The process then specifies the other 
via node (Le., the via node that is on a different layer from the Previous_Particle) as the 
Previous_Particle. When the process specifies (at 7110) the via nodes, it also specifies a location for 
the via nodes as the process 3500 performs triangulation based on these nodes at 3544. In some 
embodiments, the process defines the x- and y- coordinates of the via nodes that it creates for a 
particular hole as a location in the hole's polygon with a minimum F PLF value. A hole's F PLF 
represents the actual cost at the hole of the cheapest path that traverses from a source to the target 
through the series of particles (including the hole) that the path search process identified between the 

source and target. A particle's F PLF is different than its F function. The computation of a particle's 

F PLF is described below in next sub-section. 

When the CurrentJDrop's particle is a wall, the wall is necessarily a target particle. In this 

situation, the process 7100 might create (at 7110) a Steiner to embed topologically an expansion to 

the wall, which represents a previously specified portion of the route of the net that is currently being 

routed. Before specifying a Steiner for such an expansion, the process 7100 determines whether the 

wall is adjacent to a previously defined Steiner. If so, the process does not create a Steiner for the 

expansion to the wall, but rather specifies the previously defined Steiner as the Previous_Particle. On 
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the other hand, if no Steiner neighbors an expansion to a wall, the process 7100 (1) creates two walls 
and a Steiner in between the two walls, (2) replaces the previous wall (i.e., the wall to which the path 
expanded) with the created two walls and Steiner, and (3) specifies a wall that connects the created 
Steiner to the Previous_ParticIe. 

Next, if the process 7100 defined one or more new particles at 7110, the process associates 
(at 7115) these particles with their respective spaces. If the process specified a wall at 7110 that 
divided a previously defined space into smaller spaces, the process 7100 (at 7115) defines the 
smaller spaces, sets the holes of the newly defined spaces, and associates the particles of the 
previously defined space with the newly defined smaller spaces. The process 7100 uses the process 
1700 illustrated in Figure 17 to set the holes of the newly defined spaces. 

When the process 7100 inserts a joint in an edge at 7110, it specifies (at 7120) as dirty all 
previously defined joints on that edge. Next, at 7125, the process determines whether it has 
examined all the drops in the series of drops between the source and target sets. If not, the process 
specifies (at 7130) the Current.Drop's previous drop as the Current_Drop. It then transitions back to 
71 10 to repeat operations 71 10-7125 for this drop. When the process determines (at 7125) that it has 
examined all the drops between the source and target sets, it terminates. 
(9) Loosely Geometrizing Topological Route 
At 3546, the process 3500 defines a loose geometric path for the topological path that it 
specified at 3542 and possibly modified at 3544. This entails specifying a location of each joint and 
Steiner on the path. For each joint, the process (at 3546) (1) identifies a G PLF and an H PLF across 
the joint's domain on the joint's edge (i.e., across the portion of the joint's edge that can contain the 
joint), (2) adds these two PLFs, and (3) identifies the location of the joint at a minimum of the added 
function. A joint's G PLF specifies the cost of the topological path that reached the joint from a 
source particle through the set of path particles between the joint and the source. This cost is a PLF 
that is defined across the joint's domain. This G PLF can be computed in the same manner as 
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discussed above for G PLFs of drops for a path search. A joint's H PLF specifies the actual (not 
estimated) distance from the joint's domain to the target through the topological path's sequence of 
path particles between the joint and the target. This H PLF can be computed like the G PLF except 
that the only cost that the H PLF accounts for is the length cost. 
5 The topological path that the process 3500 specifies at 3542 and possibly modifies at 3544 

might include one or more Steiners. Each Steiner connects to three or more other particles (e.g., 
joints, nodes, etc.) through three or more walls. In some embodiments, the process 3500 specifies the 
location of each Steiner at the centroid of the three or more other particles to which it connects (e.g., 
at the centroid of the joint and nodes that it connects to). 

10 At 3546, the process also re-computes the location of the joints that were marked dirty during 

the topological embedding and edge flipping operations of 3542 and 3544. This computation might 
also affect the location of Steiners on the paths with the dirty joints. To recompute the location of the 
dirty joints, the process 3500 uses one operation in addition to the three operations discussed above 
for the joints of a newly defined topological path. Specifically, for each topological route that has 

15 one or more dirty joints, the process identifies a path segment that contains all the dirty joints of the 
route and that terminates with two immovable particles (e.g., two nodes). For each joint that is 
between the two immovable ends of the path segment, the process (at 3546) (1) identifies a G PLF 
and an H PLF across the joint's domain on the joint's edge, (2) adds these two PLFs, and 
(3) identifies the location of the joint at a minimum of the added function. If the added function has a 

20 minimum segment, then the process selects a point on the segment such that the route being 
embedded does not cross any other route that was geometrically embedded before the current route 
during the current iteration through JG46. Also, the process identifies the location of each Steiner 
that falls within the two immovable ends of the path segment based on the above-described centroid 
approach. 

25 The locations that the geometric embedding operation specifies are not necessarily fixed 
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locations. The locations are specified in order to enable the Q* search to perform more accurate 
wirelength calculations while generating the lowest-cost paths. One of ordinary skill will realize that 
other embodiments do not specify such locations for the topological routes identified by the Q* 
topological routing engine. 

One of ordinary skill will realize that other embodiments might utilize the Q* engine 
differently. For instance, the embodiments described above propagate cost functions to one- 
dimensional lines and two-dimensional surfaces. Other embodiments might propagate cost functions 
to and from (1) a set of associated collinear points (i.e., a set of points that can be connect by a single 
line passing through the points) that form a one-dimensional state and/or (2) a set of associated co- 
planar points (i.e., a set of points that fall within a two dimensional surface) that form a two- 
dimensional state. 

Some embodiments propagate functions to a set of associated points by examining each point 
independently. Specifically, for each point in the associated set, some embodiments compute a PLF- 
value by identifying the wedge and/or channel containing the point, and then computing the PLF- 
value according to the wedge and/or channel rules described above for 5415 of Figure 54. 

Some embodiments propagate functions from a set of associated points by (1) using the 
above-described approach of Figures 46-51 to identify wedge and channels vectors that emanate 
from the set of points, and (2) identifying the destination state's cost function based on the identified 
channels and wedges. Specifically, the location of each wedge is identified by identifying the point 
that is the best point in the set for the wedge. Also, a channel vector might emanate from a point in 
the set that is on a line segment connecting two points from which two adjacent wedges emanate 
(i.e., from which two wedges with parallel vectors emanate). Such channel vectors can be identified 
based on the approach described above by reference to Figures 49 and 50. After identifying the 
wedge and channel vector emanations from the set of points, the destination state's cost function can 
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be specified according to the approaches described above for destination points, lines, surfaces, and 
set of points. 

2. IDA* Topological Routing Engine 
a. Overview 

Like the Q* topological routing engine, the EDA* topological engine (1) initially decomposes 
each layer of the received sub-region to obtain a decomposed graph that has several nodes, edges, 
and faces, and then (2) uses the nodes, edges, and faces of the generated graphs to define topological 
routes in the received sub-region. 

(1) Overview of the Decomposition Operation. 

Like some embodiments described above, the IDA* engine in the embodiments described 
below uses a triangulation operation to decompose each layer of the received sub-region. In other 
words, the faces resulting from the decomposition operation are triangles and therefore have three 
edges. Each edge is defined to be between two nodes. In the embodiments described below, the 
nodes in each layer's of the decomposed sub-region are defined at the obstacle-geometry vertices, 
pin-geometry vertices, and the four corners of the sub-region on each layer. At the end of the 
triangulation operation, the IDA* topological engine adds, to the sub-region definition, a graph data 
structure for each layer of the sub-region. Each graph data structure includes several node, edge, and 
face data structures. 

Figure 72 illustrates a data structure 7200 that defines a face. As shown in this figure, a face 
is defined to refer {e.g., to point) to its 3 edges and its 3 nodes. The face also includes up to two 
references (e.g., two pointers) to up to two topological items (called face items), which as described 
below are used to define topological vias and Steiner points. Some embodiments limit the number of 
face items in a face to two in order to improve run-time speed. Other embodiments, on the other 
hand, might allow a face to have more than two face items, or might limit a face to have only one 
face item. 
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Figure 73 illustrates a data structure 7300 that defines an edge. An edge can be shared 
between two faces, or it can simply be part of one face. Accordingly, as shown in Figure 73, the 
edge data structure 7300 has two references, each of which can be assigned to refer to a face of the 
edge. 

The edge data structure also specifies the capacity of the edge. It further specifies the wire 
flow across each edge. This flow equals the width of the nets crossing the edge plus the spacing 
between the crossing nets and between the nets and the edge nodes. This data structure also has a 
Boolean flag to indicate whether the edge is a constrained edge. This flag is used during the 
triangulation when the edge is specified as a constrained edge, as described above. 

In addition, the edge data structure 7300 includes a linked list of topological items on the 
edge. During triangulation, this linked list is initialized to start with one of the edge's nodes and end 
with the other. When topological routes are inserted in the triangulated graphs, a topological item 
(called an edge item) is inserted between the end nodes of the edge's linked list for each topological 
route that crosses the edge. The edge items of several topological routes that intersect an edge are 
inserted between the edge nodes in the order that their corresponding topological routes cross the 
edge. This ordering will be further explained below by reference to Figures 77 and 78. The data 
structure for a node will also be described below. 

(2) Embedded Topological Routes 

After triangulation, the IDA* topological engine embeds multi-layer topological routes in the 
triangulated graphs. It uses the nodes, edges, and faces of these graphs to define the topological 
direction of these routes. Specifically, the IDA* topological engine defines a topological route as a 
linked list of topological items that represent various points along the route. These topological items 
include nodes, edge items, and face items. Nodes and face items can serve as via locations, Steiner 
points, or both. Nodes can also serve as end points of a topological route. 
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Figures 74-76 illustrate the data structure of nodes, edge items, and face items. As shown in 
Figure 74, a node data structure 7400 includes a net identifier, which, when the node is on a net's 
topological route, specifies this net. This structure also has one or more planar-path references. When 
the node is part of a topological route, a planar-path reference refers to a planar topological item (i.e., 
an item on the same layer as the node) that is adjacent to the node in its topological route. The node 
data structure can have more than one planar-path reference (e.g., has more than one such reference 
when the node serves as a Steiner point). 

In the embodiments described below, a node can serve as the location for a via. Accordingly, 
the node data structure 7400 includes a pair of via-path references. When a topological via is placed 
at the location of a node, one or both via references refer to face items or nodes in the layer above 
and/or below. 

The node data structure further includes a list of edges connected to the node. For each edge, 
it includes an edge reference to the next or previous topological item (i.e., node or edge item) on the 
edge. It also has a reference to the geometry of the node, and a vertex number that identifies the point 
in the geometry corresponding to the node. The node structure further specifies the location of the 
node. 

Figure 75 illustrates the data structure for an edge item. This data structure has a reference to 
its edge. In addition, like a node data structure, an edge item's data structure includes a net identifier 
that specifies the net of the edge-item's topological route. This structure also includes a pair of edge 
references that refer to the next and previous topological item (i.e., node or edge item) on its edge. 
The edge-item data structure also has a pair of planar-path references that refer to the planar 
topological items (i.e., items on the same layer as the edge item) that are adjacent to the edge item in 
its topological route. 

Figure 76 illustrates the data structure for a face item. The face-item data structure 7600 has 
a reference to its face. In addition, like the data structure of a node or an edge item, a face item's data 
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structure includes a net identifier that specifies the net of the face item's topological route. This 
structure can also includes up to three planar-path references. Each planar-path reference refers to a 
planar topological item (i.e., an item on the same layer as the face item) that is adjacent to the face 
item in its topological route. The face-item data structure can have up to three planar-path references 
since the face-item can serves as a Steiner point for up to three topological paths. 

In the embodiments described below, a face item can serve as the topological position of a 
via. Accordingly, the face item data structure 7600 includes a pair of via-path references. When a 
topological via is placed at a face item, one via-path reference can refer to a node or face item in the 
layer above the face item's layer, while the other via-path reference can refer to a node or face item in 
the layer below the face item's layer. The face-item data structure also specifies a bounding convex 
polygon that approximates the region where the center of the face item can be legally placed in the 
face containing the face item. In some embodiments, this data structure also specifies a set of 
constraining points and distances for the face item in its layer. The use of the bounding polygon and 
the constraining points and distances will be further described below. 

The layer of the topological items can be identified as follows. A node's layer can be 
identified by referring to its geometry. A face's or edge's layer corresponds to the layer of its nodes. 
In addition, a face item's layer corresponds to the layer of its face, and an edge item's layer 
corresponds to the layer of its edge. 

Figures 77 and 78 illustrate one example of topological routes. Figure 77 presents two 
topological routes 7770 and 7772. The route 7770 is a multi-layer topological route that (1) starts at 
node 7702 of layer 2, (2) intersects edges 7720, 7728, and 7736 on layer 2, (3) vias up to layer 3 
through face items 7738 on layer 2 and 7750 on layer 3, (4) intersect edge 7754 on layer 3, and 
(5) terminates on node 7758 on layer 3. This route's intersections with edges 7720, 7728, 7736, and 
7754 are specified by edge items 7716, 7726, 7732, and 7756. 

The route 7772 is a topological route that only traverses layer 2. This route (1) starts at node 
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7704, (2) intersects edges 7712, 7720, 7728, and 7736, and (3) terminates on node 7740. This route's 
intersections with edges 7712, 7720, 7728, and 7736 are specified by edge items 7710, 7718, 7724, 
and 7734. 

Figure 78 illustrates the data structure representation of these topological routes 7770 and 
7772. Each of these routes is specified by two sets of linked lists. One set includes a linked list of 
path references, and the other set includes linked lists of edge references. The path linked list of route 
7770 starts at node 7702, goes through edge items 7716, 7726, 7732, face items 7738 and 7750, edge 
item 7756, and terminates on node 7758. The path linked list of route 7772 starts at node 7704, goes 
through edge items 7710, 7718, 7724, and 7734, and terminates on node 7740. 

As shown in Figure 78, the nodes and edge items are inserted on their edge's linked list in the 
order that they are placed on their respective edge. For instance, the edge list 7764 of edge 7728 
starts with node 7714, traverses through edge items 7726 and 7724, and then terminates on node 
7722. This order of these nodes and edge item on edge list 7764 matches the order (illustrated in 
Figure 77) of these nodes and edge items on edge 7728. 

b. Overall Flow of IDA* Topological Router 

Figure 79 illustrates a process 7900 that provides the overall flow of the IDA* topological 
engine in some embodiments of the invention. As shown in this figure, the process initially 
triangulates (at 7905) the sub-region defined at 205. This triangulation operation is similar to the one 
described above by reference to Figure 9. There are few minor differences, however. For instance, 
the IDA* triangulation in some embodiments decomposes the port geometries. Also, the IDA* 
triangulation does not define holes or spaces. United States Patent Application entitled "Method and 
Apparatus for Generating Multi-Layer Routes," and having the serial number 10/076,121, discloses a 
triangulation technique that is used for the IDA* engine in some embodiments. This application is 
incorporated herein by reference. 
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After the triangulation, the process groups (at 7910) the nets in the sub-region. Different 
embodiments group the nets in the sub-region differently. Some embodiments group the nets based 
on a clustering approach. Some of these embodiments cluster the nets by (1) defining, for each net, a 
three-dimensional bounding box that contains the net's pins, (2) pair-wise interacting the bounding 
box of each net with each of the other nets, (3) computing the volume of intersection between the 
intersections, and (4) clustering each net with the other nets that have the greatest overlapping 
volume (i.e., that are likely to have the most amount of overlap). Some of these embodiments order 
the nets in each group. For instance, some embodiments sort each group based on entropy (e.g., 
based on descending order of entropies, where a high entropy net has lots of equivalently good 
options, while a low entropy net has only a few good routes). 

Next, the process selects (at 7915) one of the groups identified at 7910 for routing. Different 
embodiments process the selected group of nets differently. The embodiments described below solve 
all the nets in the selected group before selecting another group of nets. 

Other embodiments, however, might process the nets differently. For instance, some 
embodiments might order the nets based on an entropy value (e.g., might order the nets in a 
descending order of entropies), and select the group of nets based on this order. Some of these 
embodiments then use a two-tiered approach in selecting the group of nets at 7915. First, from the 
top of the set of unsolved nets in the ordered netlist, these embodiments select a fixed number (e.g., 
20) nets. When this list does not have the fixed number of nets remaining, the process selects all the 
remaining nets. Second, after selecting the group of nets, these embodiments return the nets in the 
selected group that they cannot solve within an acceptable cost in a pre-defined number of tries. 
Accordingly, these embodiments statically try to select a fix-number of nets each time at 7915, but 
then dynamically return some of the selected nets as unsolved nets when the selected group proves 
too hard to route entirely at the same time. The returned unsolved nets are then treated like the 
remaining unsolved nets. In other words, these embodiments perform the two operations listed above 
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on all the remaining unsolved nets (i.e., from the list of unsolved nets that includes the returned 
unsolved nets, these embodiments (1) select and try to process a number of nets, and (2) return any 
nets that they cannot solve within an acceptable cost). 

In the embodiments described below, the IDA* topological engine uses an IDA* solving 
engine that deterministically searches for the best routing solution for the selected group of nets. An 
IDA* searching technique involves searching the solution space in an iteratively deepening manner, 
up to a particular depth limit. 

Accordingly, at 7920, the process specifies a particular depth limit for the IDA* solving 
engine. Some embodiments compute the depth limit as the sum of the cheapest topological routes for 
all the unsolved nets in the sub-region. One manner of computing the cheapest topological routes 
will be explained below by reference to Figure 106. 

Some embodiments generate the cheapest route for each unsolved net by ignoring the routes 
of the other unsolved nets. Accordingly, these embodiments generate the cheapest routes for each net 
in the first group of nets being solved in an empty sub-region (i.e., a sub-region that does not contain 
the topological routes of any other nets in the group). However, the cheapest routes for the 
subsequent groups of nets are computed for a sub-region that contains the topological routes of the 
previously routed nets. 

At 7920, the process also computes a maximum depth-limit threshold. In some embodiments, 
the maximum threshold value is 50% greater than the initial depth limit specified at 7920. The 
process also specifies (at 7920) a maximum number of iterative attempts, referred to below as 
pushes, for the IDA* solving engine. 

At 7925, the process then calls the IDA* solving engine to find a topological route for all or 
some of the nets in the selected group within the specified depth limit and max pushes. The 
operation of the IDA* solving engine will be further described below. Next, the process determines 
(at 7930) whether the IDA* solving engine returned a solution for all the nets in the selected group. 
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If so, the process transitions to 7965, which will be described below. If not, the process increments 
(at 7935) the depth limit. Some embodiments increment the depth limit by 10%. Accordingly, the 
solution found by incrementing the depth limit in this manner is always with 10% of the optimal 
solution. 

The process then determines (at 7940) whether the depth limit exceeds a maximum value. If 
not, the process transitions back to 7925 to direct the IDA* solving engine to find a solution within 
the incremented depth limit. On the other hand, if the process determines (at 7940) that the current 
depth limit is greater than the maximum threshold value, the process determines (at 7945) whether 
the IDA* solving engine returned an empty solution set in response to the call at 7925. 

As further described below, the IDA* solving engine returns an incomplete non-empty 
solution when it can only find topological routes for some of the nets in the selected group. If the 
process determines (at 7945) that the IDA* solving engine has returned an incomplete non-empty 
solution, it accepts (at 7955) the topological routes for some of the nets that the IDA* solving engine 
returned. The process then returns to 7920 to try to find topological routes for the remaining 
unsolved nets of the group selected at 7915. In other words, the process (1) specifies (at 7920) a 
depth limit, a maximum depth-limit threshold, and a max number of pushes for the remaining 
unsolved nets, (2) calls (at 7925) the IDA* solving engine for these nets, and then (2) based on the 
solutions returned by the IDA* solving engine, performs some or all the operations 7930-7945 as 
described above. 

On the other hand, if the process determines (at 7945) the IDA* solving engine did not find a 
solution for any net, the process selects (at 7950) one of the nets in the selected group and creates a 
solution for this net based just on this net's cheapest topological route. The process then determines 
(at 7952) whether it has found a topological route for all the nets in the group identified at 7915. If 
not, the process returns to 7920 to try to find topological routes for the remaining unsolved nets of 
the group selected at 7915. Like the transition back to 7920 from 7955, this transition results in the 
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process (1) specifying (at 7920) a depth limit, a maximum depth-limit threshold, and a max number 
of pushes for the remaining unsolved nets, (2) calling (at 7925) the IDA* solving engine for these 
nets, and then (2) based on the solutions returned by the IDA* solving engine, performing some or 
all the operations 7930-7945 as described above. 

On the other hand, if the process determines (at 7952) that it has found a solution for all the 
nets in the group identified at 7915, the process (at 7965) specifies a shape for each face item that the 
IDA* solving engine defined while routing the nets. The face items can form topological vias or 
Steiner points. Different embodiments use different shapes for face items. For instance, the face 
items can be hexagonal, octagonal, circles, squares, diamonds, etc. Some embodiments select the 
shape of the face items based on the wiring model that the geometric engine uses. Several of these 
shapes, and the uses of these shapes for different wiring models, are described in the above- 
incorporated United States Patent Application 10/076,121. 

After assigning face-item shapes, the topological-routing process 7900 then triangulates (at 
7970) the sub-region layers based on the face items. This triangulation is performed by inserting a 
node in the triangulated graph in the location of each face item {e.g., in the center of each face item 
shape defined at 7965), and then triangulating each face that contains a face item based on the 
location of the inserted node. In some embodiments, the process 7900 specifies the location of each 
face item at 7965. For instance, at 7965, the process 7900 can identify the location of a non-via face 
item at the center of its bounding polygon, and identify the location of a via-forming face item at the 
center of the bounding polygon for its via (i.e., at the center of the intersection of the bounding 
polygons of the items that form the via). In other embodiments, the location of each face item is 
specified while routing the group of nets selected at 7915. For instance, as further described by 
reference to Figures 91, a face item's location can be identified and modified by an optimization 
process. Some embodiments might also insert nodes about the periphery of a face item shape in order 
to triangulate the face-item shape into smaller triangles. 
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At 7970, the process also modifies the topological description of the routes that cross the new 
triangulated edges. After this triangulation, the process (at 7975) performs a follow-up edge-flipping 
operation. If the process ends up flipping any edges at 7975, it also modified at 7975 the topological 
description of the routes that crossed the old or now cross the new edge. At 7980, the process 
determines whether there are any unsolved nets in the sub-region. If so, the process returns to 7915 to 
select the next group of nets for solving. If not, the process ends, 
c. IDA* solving Engine 

(1) Jointly Routing Groups of Nets with an IDA* Search Engine 
As mentioned above, the IDA* topological routing engine calls (at 7925) its IDA* solving 
engine to identify topological routes for a group of nets. In the embodiments described below, the 
IDA* solving engine uses an IDA* searching technique to try to identify the optimal set of routes for 
the group of nets. In other words, the IDA* solving engine traverses through the solution space to 
deterministically identify the best possible combination of topological routes for the group of nets. 

The IDA* solving engine traverses the solution space in an iterative-deepening manner up to 
a particular depth limit. Figure 80 pictorially illustrates an example of this engine's BDA*-searching 
operation for a set of three nets. As shown in this figure, the IDA* solving engine initially identifies 
several topological routes 8005a-05c for a first net. One advantage of using topological routes is that 
the IDA* solving engine can define fewer possible routing solutions for each net at this stage. This is 
because there are fewer topological routes than geometric routes for a net, as one topological route 
can represent a plethora of diffeomorphic geometric routes. 

After identifying the solutions for the first net, the IDA* solving engine selects a topological- 
routing solution for the first net. At each level of iteration, the search selects the solutions according 
to a quality metric. In the embodiments described below, this quality metric is the length of the 
topological routes. Other embodiments might use different quality metrics. In the example illustrated 
in Figure 80, the topological routes for each net get longer from left to right (i.e., the solutions get 
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worse going from left to right). 

Accordingly, in the example illustrated in Figure 80, the first selected solution for the first 
net is the best solution 8005a. Next, the IDA* solving engine determines whether embedding the 
selected solution 8005a in the received sub-region (i.e., the sub-region supplied by the EDA* 
topological engine at 7940) would cause the routing to exceed the depth limit specified by the IDA* 
topological engine. In the example illustrated in Figure 80, this embedding does not make the 
routing exceed its depth limit. Hence, the engine embeds the selected solution 8005a in the sub- 
region, and then generates several solutions 8010a-c for the second net in the sub-region that 
contains the embedded solution 8005a of the first net. 

After identifying the solutions for the second net, the IDA* solving engine selects the best 
solution 8010a for the second net, and determines whether embedding this selected solution in the 
current sub-region (i.e., the sub-region containing the embedded topological route 8005a) would 
make it exceed the depth limit specified by the IDA* topological engine. In the example illustrated in 
Figure 80, this embedding does not make the routing exceed its depth limit. Hence, the engine 
embeds the selected solution 8010a in the sub-region, and then generates several solutions 8015a-c 
for the third net in the sub-region that contains the embedded routes 8005a and 8010a of the first and 
second nets. 

After examining all the solutions 8015a-8015b for the third net in the sub-region that 
contains the embedded routes 8005a and 8010a of the first and second nets, the IDA* solving engine 
discover that all these solutions 8015a-8015b would make the route exceed the depth limit specified 
by the IDA* topological engine. Consequently, the combination of the best solutions 8005a and 
8010a for first and second nets did not lead to an acceptable routing solution. 

The IDA* solving engine then determines whether it has examined all the previously defined 
solutions for the net at the previous level of the search. At this stage, this previous example is net 2, 
and the IDA* solving engine has not yet examined other routing solutions 8010b-8010c of net 2. 
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Accordingly, the IDA* solving engine selects the next best solution 8010b for net 2 in the sub-region 
that contains solution 8005a for net 1. It then embeds this route 8010b as this embedding does not 
violate the depth limit, and then generates several routes 8020a-8020b for the third net in the sub- 
region that the embedded routes 8005a and 8010b of the first and second nets. 

However, the IDA* solving engine cannot embed any of these routes 8020a-8020b. In this 
example, the IDA* solving engine then determines that it can embed all other routing solutions of net 
2, but each time it determines that it cannot embed any of the resulting routing solutions for net 3. 
This means that the best solutions 8005a for first net does not lead to an acceptable routing solution. 
Hence, once the IDA* solving engine determines that it has examined all the previously defined 
solutions for the net 2 at the current level of the search, it determines whether it has examined all the 
previously defined solutions for the net at the previous level of the search. At this stage, this previous 
example is net 1, and the IDA* solving engine has not yet examined other routing solutions 8005b- 
8005c of net 1. 

Accordingly, the IDA* solving engine examines the solution space that results from the 
selection of the next best solution 8005b for net 1. Eventually, in this example, the IDA* solving 
engine finds a solution for all three nets, when it determines that embedding the topological route 
8040a for the third net does not cause the routing to exceed the depth limit. This solution includes 
route 8005b for net 1, route 8030b for net 2 in the sub-region that contains route 8005b, and route 
8040a for net 3 in the sub-region that contains routes 8005b and 8030b. 

This solution will be the optimal solution or very close to it so long as the depth limit is 
increased by small increments from its lowest possible value. For instance, when the depth limit is 
increment by 10%, the solution will be within 10% of the optimal solution. 

Some embodiments described below introduce two additional constraints in the operation of 
the IDA* solving engine. First, these embodiments limit the IDA* solving-engine's IDA* search to a 
maximum number of iterative attempts, referred to below as pushes. This is to decrease the run-time 
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of the IDA* solving engine. Limiting the number of pushes might cause the IDA* solving engine to 
return a solution that is not the optimal, but this might be acceptable in some situations as a 
reasonable tradeoff for run-time speed. 

Second, the IDA* solving engine in these embodiments considers the impact of the routes 
selected for the nets in the selected group on the nets not yet selected. In other words, while IDA* 
solving the selected group of nets, the IDA* solving engine considers the routing cost of nets that 
have not yet been selected. In the embodiments described below, the IDA* solving engine generates 
a lower bound estimate on the cost of the unselected nets, and ensures that the depth limit is not 
exceeded by the sum of this estimate and the actual or predicted costs of the nets in the selected 
group. 

(2) Solving 

Figure 81 illustrates a more detailed process 8100 used by the IDA* solving engine in some 
embodiments of the invention. As mentioned above, the process 8100 is performed each time that 
the process 1500 calls it at 1525 and specifies a selected group of nets to route within a given sub- 
region, depth limit, and maximum number of pushes. As shown in Figure 81, the process 8100 
identifies (at 8102) the Current_Net as the first net in the specified group of nets. At 8102, the 
process also sets a Push.Count to 0. For the Current.Net, the process then generates (at 8108) all the 
legal routes within the sub-region that have an acceptable cost. The acceptable cost is the depth limit 
minus the total cost of the routed nets minus the estimated cost of all the nets below the Current_Net, 
plus the estimated cost of the Current_Net. 

When the specified group of nets is not the first group of nets that the IDA* solving engine 
has solved (i.e., if the topological route has called the IDA* solving engine previously to solve a 
different group of nets), the sub-region includes the topological routes of the previously solved nets. 

In addition, when the process 8100 generates the topological routes for the first net of the 
specified group, the sub-region does not include the topological route of any net within this group. 
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However, each time the process 8100 embeds a topological route of one of the earlier specified nets, 
it generates the routes of the subsequent specified nets for a sub-region that contains the routes of the 
earlier embedded nets. In other words, when the process 8100 generates the topological routes for 
specified nets other than the first specified net, the sub-region includes previously embedded 
topological routes of previous specified nets. At 8104, the process 8100 also stores the cost 
associated with each topological route that it generates. The route-generating process will be 
described below by reference to Figure 82. 

After 8104, the process determines (at 8106) whether it was able to generate (at 8104) any 
topological route for the Current_Net. If not, the process 8100 transitions to 8142, which will be 
described below. On the other hand, if the process generated (at 8104) one or more topological routes 
for the Current_Net, the process 8100 stores (at 8108) the generated routes. 

Some embodiments store the routes in an array of N data objects, where N corresponds to the 
number of nets that the IDA* solving engine is trying to solve. In other words, there is one data 
object for each of the N nets. Each data object can store the route solution pool of its net, and 
includes a pointer into this pool that specifies the solution currently being explored. This pointer 
traverses through the solution pool as the IDA* solving engine examines the solutions of its net. 

At 8110, the process selects one of the topological routes. In the embodiments described 
below, the process selects the shortest routes first. The process then increments (at 8112) the 
Push.Count by 1. Next, at 8114, the process embeds the selected route in the current sub-region, and 
subtracts the cost of this embedded route from the depth limit to obtain a new depth limit. At 81 16, 
the process computes a revised total cost to account for the increase in cost due to addition of the 
selected route. This revised cost equals the old total cost minus the estimated cost of the Current.Net 
plus the actual cost of the inserted route. 

Next, the process determines (at 8118) whether the total cost will exceed the depth limit after 
the cost increase in the downstream nets is added to the total cost. In some embodiments, the process 
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adds the increase in the cost of each downstream net in a loop that checks after the addition of each 
downstream net's delta whether the addition has caused the total cost to exceed the depth limit. Such 
an approach allows the process 8100 to quickly identify that the route inserted at 81 12 has caused the 
depth limit to be exceed. In some embodiments, the process 8100 uses the cheapest-route calculation 
process of Figure 106 to compute the new cost of a downstream net. It also stores the downstream 
costs of each net (e.g., stores these costs in a stack), and computes the net's delta downstream by 
subtracting the appropriate stored cost from the computed cost. 

When the process determines (at 8118) that the total cost exceeds the depth limit, it 
transitions to 8132 to decide whether it has also exceeded its maximum number of pushes. This 
determination will be further described below. On the other hand, if the process determines (at 8118) 
that the inserted route did not cause the depth limit to be exceeded, the process determines (at 8120) 
whether the insertion of the selected route at 8116 resulted in the best solution thus far. When the 
process reaches 8120, it has identified routes for a current set of one or more nets in the specified 
group. Hence, at 8120, the solution for the current set of nets is the best solution (1) if the process 
has not previously identified a solution for more than the current set of nets, or (2) if the process has 
previously identified a best solution for as many nets of the specified group of nets, but the sum of 
the cost of the routes in the previous solution is more the sum of the costs of the routes in the current 
solution. 

If the process determines (at 8120) that the current set of embedded routes for the current set 
of nets do not provide the best solution obtained thus far, it transitions to 8124, which is described 
below. Otherwise, the process stores (at 8122) the current set of embedded routes as the best 
solution, and then transitions to 8124. 

At 8124, the process determines whether it's Push_Count equals the maximum number of 
pushes that it received from the IDA* topological-routing process 1500. If the Push.Count equals 
the maximum number of pushes, the process returns (at 8130) the best solution that it recorded 
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during its search, and then ends. 

On the other hand, if the Push_Count is not equal to the maximum number of pushes, the 
process determines (at 8126) whether the CurrentJMet is the last net of the specified group of nets. If 
so, the process returns (at 8130) the best solution that it recorded during its search, and then ends. 

If not, the process (at 8132) sets the Current_Net to be the next net of the specified group of 
nets. The process then transitions back to 8104 to generate all legal routes for this new Current_Net, 
and then perform the subsequent operations to try to embed one of these generated routes in the sub- 
region. 

As mentioned above, the process 8100 transitions to 8132 when it determines (at 8118) that 
the total cost for inserting the selected route exceeds the depth limit. At 8132, the process determines 
whether its Push_Count equals the maximum number of pushes received from the IDA* topological- 
routing process 1500. If so, the process returns (at 8130) the best solution that it recorded during its 
search, and then ends. 

On the other hand, if the process determines (at 8132) that the Push_Count does not equal the 
received maximum number of pushes, the process removes (at 8134) the route inserted at 81 14. The 
process then determines (at 8136) whether it has examined all the topological routes generated for 
the Current_Net at 8104. If not, the process transitions back to 8110 to select the next best route 
(e.g., next shortest route) for the Current_Net from the array described above, and then perform some 
of subsequent operations described above for the newly selected route. 

If the process determines (at 8136) that has examined all the topological routes generated for 
the Current_Net, the process determines (at 8138) whether the Current_Net is the first net in the 
specified group of nets. If so, the process has failed to embed any of the routes of the first net in the 
specified group, and hence returns (at 8144) its failure to route to the topological-routing process 
1500. 

If not, the process clears (at 8140) the Current.Net's solution pool in the above-described 
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array. The process then defines (at 8142) the net before the Current_Net on the specified group as the 
new Current_Net. It next determines (at 8136) whether it has already explored all the routes for the 
Current.Net defined at 8142. Depending on this determination, the process transitions either to 8110 
or 8138, as described above. 

(3) Route Generation 
Figure 82 illustrates a process 8200 that the solving engine uses at 8104 to generate 
topological routes for a net. Like the solving engine, the route generating process uses an IDA* 
search technique to identify the several routes for each net. 

As shown in Figure 82, the process 8200 initially identifies (at 8202) one or more pairs of 
pins for the Current_Net in the sub-region. It then selects (at 8204) a pair of pins for the net. The 
process next identifies (at 8206) a list of source and target port-geometries for the selected pair of 
pins of the Current_Net. 

The process then identifies (at 8208) a DepthJJmit for generating routing path between the 
source and target geometries. In the embodiments described below, the Depth_Limit is set to the 
shortest distance between the closest source and target geometries. After 8208 , the process defines 
(at 8210) Solution_Count to be equal to 0. 

Next, the process generates (at 8212) all legal paths between the source and target geometries 
identified at 8206. The path-generating process will be described below by reference to Figure 84. 
After 8212, the process determines (at 8214) whether it was able to generate (at 8212) any legal path 
between the source and target geometries. If not, the process transitions to 8216, which will be 
described below. 

If so, at 8226, the process (1) records the generated path or paths and the cost for each path, 
and (2) increments Solution_Count by the number of paths generated at 8212. The process next 
determines (at 8228 ) whether SoIution.Count is less than the desired number of solutions. If so, the 
process transitions to 8216, which will be described below. Otherwise, when the number of solutions 
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exceeds the desired number of solutions, the process transitions to 8224, which will be described 
below. 

As mentioned above, the process can transition to 8216 from 8214 or 8228. At 8216, the 
process increments the Depth.Limit. The process then determines (at 8218) whether the 
Depth.Limit exceeds a maximum depth limit. In some embodiments, the initial depth limit is the 
Euclidean distance between closest source and target, and the depth limit is incremented at each 
iteration by 1.1 times the cost of cheapest branch pruned for cost during previous iteration. 

If the Depth_Limit does not exceed the maximum depth limit, the process transitions to 8212 
to generate legal paths between the source and target geometries for the Depth.Limit incremented at 
8216. Otherwise, the process determines (at 8220) whether it was able to record any solutions for the 
pair of pins selected at 8204. 

If the process determines (at 8220) that it was not able to record solutions for the selected pair 
of pins, it returns (at 8222) its failure to find routes for the net, and then ends. Otherwise, the process 
determines (at 8224) whether it has examined the last pin-pair identified at 8202. If not, the process 
transitions back to 8204 to select the next pin-pair. 

On the other hand, if the process determines that it has examined all the pin-pairs, the process 
(at 8232) identifies up to K cheapest routes for the net, and records the costs of these routes. When 
the net has more than two pins in the sub-region, the process 8200 generates more than one set of 
paths for more than one set of pin-pairs. In such a situation, the process has to generate (at 8232) the 
K-cheapest combination of paths for the different pin-pairs to generate the K-cheapest routes. When 
the net has three pins in the sub-region and the process has identified two sets of paths for two pin- 
pairs, the process selects (at 8232) the cheapest combination of paths from the two sets. On the other 
hand, when the net has more than three pins in the sub-region and the process has identified more 
than two sets of paths for more than two pin-pairs, the process (1) initially computes the minimum 
spanning tree ("MST") for the net pins, and (2) selects the cheapest combination of paths for the pin- 
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pairs that are the endpoints of the MST's edges. 

After 8232, the process inserts (at 8234) any applicable Steiner nodes in the generated routes. 
In some embodiments, the process does this by passing through the edges intersected by each 
generated route to determine whether any two consecutive edge items are for the same route. Two 
such edge items are illustrated in Figure 83A. When the process identifies two such edge items, it 
removes one of the edge items and in its place inserts a face item as a Steiner node for merging the 
two paths as shown in Figure 83B. 
After 8234, the process ends. 

(4) Path Generation 
Figure 84 illustrates a process 8400 for generating (at 8212) paths between one or more 
sources and one or more targets for a selected pin-pair. As shown in this figure, this process initially 
defines (at 8402) one of the source nodes as the starting point of a path. At 8402, the process also 
defines the starting point as the Current_Point, and sets the Path_Cost to 0. 

The process then identifies (at 8404) all possible path expansions from the Current.Point. 
Figures 85-89 illustrate the possible expansions from edge items, nodes, and face items. Figure 85 
illustrates that an edge item can expand (1) towards the opposing edges of its face, (2) towards its 
face's node that is not on its edge, and (3) towards its face's face item. 

Figure 86 illustrates that a node can expand (1) towards the opposing edges of two faces that 
abut the node, (2) towards the other nodes of these two faces, (3) towards face items of these two 
faces, and (4) towards one or more nodes and face items above or below it. A node can via up or 
down to more than one node or face item, when more than one triangulated graph is above or below 
the node. Some embodiments explore each potential via expansion possibility of a node. 

Figure 87 illustrates that a face item can expand towards the three nodes and edges of its 
face. Some embodiments do not allow a planar expansion from a face item (/.<?., an expansion to an 
item in the same face as the face item) when the face item was reached through a planar expansion. 
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In other words, these embodiments only allow a via expansion from a face item, when the face item 
was reached through a planar expansion. Like via-path expansion from a node, a path can via up or 
down from a face item. Some embodiments do not allow a via expansion in a particular direction 
(e.g., down) from a face item, when the face item was reached in a direction opposite to the 
particular direction (e.g., up). 

Also, like a node, a face item can via up or down to more than one node or face item above or 
below it, since more than one face of the above/below triangulated graph can be above or below the 
face item. As with a node, some embodiments explore each potential via expansion possibility of a 
face item. Also, as illustrated in Figures 88 and 89, when the face item serves as the destination of a 
via from another layer, the face item has more expansion possibilities if topological routes intersect 
its face. For instance, in Figure 88, one net runs through the destination face item's face, and the face 
item has eight expansion possibilities, with five of them being on one side of the router and the other 
three being on the other side of the route. In Figure 89, two nets run through the destination face 
item's face, and the face item has ten expansion possibilities. 

When it is time to expand from a destination face item, some of the expansion possibilities 
might be quickly eliminated if there is no space for the via to be located in a region that gives rise to 
the expansion possibilities. For instance, in Figure 89, the only viable via location might be in the 
region between the two crossing routes. In such a circumstance, the three and five expansion 
possibilities on the other sides of these route will not be explored. 

After identifying the expansions from the Current.Point at 8404, the process selects (at 8406) 
one of the expansions identified at 8404. The process then determines (at 8408) whether this 
expansion is legal. The legality check will be further described below. 

If the process determines (at 8408) that the expansion is not legal, the process transitions to 
8428, which will be described below. On the other hand, if the expansion is legal, the process (at 
8410) defines the Current.Point to be the selected expansion point, and adds this new CuirentJPoint 
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to the path. Next, the process calculates (at 8412) the cost due to the expansion, and adds this 
calculates cost to the Path_Cost. In some embodiments, the process computes the cost of each legal 
expansion (which when amalgamated for complete routes define the cost of the resulting routes) by 
using a shortest-path approach. Such an approach is disclosed in "Geodesic Path Inside Polygons," 
by Simon Mak and Sung Soo Kang, which can be found at: 

www.cs.mcgi 1 1 .ca/~skang5/cs507/welcom.htm. 
Essentially, this approach can be used to compute the shortest path within a polygon defined by the 
sequence of edges intersected by the route. As this approach simply computes the shortest distance, it 
does not disfavor or penalize any potential wiring directions for one preferred direction. It should be 
noted that some embodiments compute this shortest distance based on the routing directions of their 
wiring model. For instance, some embodiments that expand a route on an octilinear layer calculate 
the cost of an expansion as the shortest distance that can be traversed by only using horizontal, 
vertical, and ±45° interconnect lines. 

At 8412, the process also adds a via cost for a planar expansion to a face item. In the 
embodiments that allow stack vias, the process also adds a via cost for a via expansion (i.e., non- 
planar expansion) to a face item when the previous expansion was also a via expansion to a face 
item. As mentioned above, in some embodiments, this via cost is proportional to the metric cost that 
a via introduces into the design. Specifically, a via introduces a variety of costs in a design, and the 
via penalty can take into account all these costs. For instance, from a resistance or delay point of 
view, a via might cost fifty times more than the resistance or delay of a unit length of a wire. In this 
situation, the via can be costed as 50 units of wirelength (i.e., the via can be counted as a wire that is 
50 units long). Alternatively, each via might reduce the manufacturing yield by some amount. This 
reduction can also be accounted for by equating it to a metric cost and adding this cost to the total 
cost. The via cost can be different for vias that are between different layers if the cost introduced by 
vias between different layer pairs is different. 
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After 8412, the process determines (at 8414) whether the Path_Cost exceeds the depth limit. 
If so, the process transitions to 8422, which will be described below. If not, the process determines 
(at 8416) whether the Current_Point is a target node. If it is not, the process transitions back to 8404 
to generate all expansions for the Current_Point. However, if the Current_Point is a target node, the 
process (at 8420) records the path with its Path_Cost, and increments the Solution_Count by one. 
The process then transitions to 8422. 

At 8422, the process removes that Current_Point from the path. It then defines (at 8424) the 
Current_Point to be the last topological item in the path. Next, the process subtracts (at 8426) from 
the Path_Cost the cost of the expansion that was removed at 8422. At 8428, the process determines 
whether it has examined all the expansions of the Current_Point. If not, the process returns to 8406 
to select another expansion. If so, the process determines (at 8430) whether the Current_Point is a 
source node. 

If the Current_Point is not a source node, the process transitions to 8422 to remove it from 
the path. However, if the Current_Point is a source node, the process determines (at 8432) whether 
there are any other source nodes that it has not yet examined. The process ends when it has examined 
all the source nodes. On the other hand, if there are one or more such nodes, the process (at 8434) 
defines one of such nodes as the starting point of a path, defines the starting point as the 
Current_Point, and sets the Path.Cost to 0. The process then transitions to 8404 to generate all 
expansions from this Current_Point. 

(5) Legality Check 

For some embodiments, Figure 90 provides a table that illustrates the types of legality checks 
performed for different combination of expansions. The vertical axis lists the starting points of the 
expansion, and the horizontal axis lists the destination points of the expansion. As shown in Figure 
90, the are three types of legality checking. These are: planarity, via checking, and edge capacity. 
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(i) Edge Capacity Check 

The edge capacity check is performed each time a path tries to intersect an edge. This legality 
check requires that the flow across the destination edge after the potential move not exceed this 
edge's capacity. The flow computation is like the one described above for the Q* engine. 

(ii) Via Legality Check 

As indicated in Figure 89, the via check is performed for each possible expansion. Different 
embodiments of the invention perform the via-checking operation differently. Figures 91 and 92 
illustrate two different processes for performing this checking operation. 
(1st) Optimization Technique 
Figure 91 illustrates a process 9100 that uses an optimization technique to perform the via- 
checking operation. The process 9100 initially determines (at 9105) whether the expansion face has a 
face item within it. A path expansion can be between source and destination items on the same layer 
(i.e., can be a planar expansion) or it can be between source and destination items on different layers 
(i.e., can be a via expansion). In the embodiments described below, an expansion face (1) is the face 
that contains both the source and destination items in a planar expansion, and (2) is the face that 
contains a face-item destination in a via expansion. 

In the embodiments described below, there is no expansion face when a node is the 
destination of a via expansion. The process 9100 does not specify an expansion face for a via 
expansion to a node in order to improve run-time efficiency. Otherwise, the process 9100 would 
have to examine every face that is connected to the via-destination node. Similarly, an expansion 
face for a planar expansion is defined as the face that contains both the source and destination items, 
in order to reduce the number of faces that are examined and thereby improve run-time efficiency. 

However, some embodiments could examine each face that contains a path-expansion's 
destination item, as this approach has several advantages. For instance, such a check would allow the 
router to detect an illegal via-path expansion when the router imposes a larger minimum-spacing 
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requirement for a destination node that serves as a via than for a destination node that does not. 

If the expansion is a via expansion to a node, there is no face to check at 9105, and the 
process 9100 thereby terminates. Similarly, if the expansion is a planar expansion to an item in an 
expansion face without a face item, then there is no via check necessary and the process ends. 

On the other hand, if the identified expansion is to a face item or, in a planar expansion, to a 
node or edge item in an expansion face with a face item, the process 9100 identifies (at 9110) 
constraints for each face item in each face involved with the path expansion. For a via expansion, a 
face is involved with a path expansion if it includes an item that forms the via. For a planar 
expansion, the expansion face is the face involved with the path expansion. 

For a face item in a face involved with a path expansion, the constraints include constraining 
points and minimum required distance to each constraining point. If the face item is moveable, the 
face item ! s set of constraints also specifies a legal x,y region that can contain the face item. Some 
embodiments allow each face to have up to two face items. When a face item is within a face that 
contains another face item, the process 9100 identifies a repulsion constraint that specifies the 
minimum required distance between the two face items. 

Figures 93-105 present several examples that illustrate the type of constraints identified by 
the process 9100. 

Constraints for A Single Face Item. 

Figure 93 illustrates a potential face item Fl in a face 9300 on a particular layer. This face 
9300 has nodes N1-N3 and edges E1-E3, and it is intersected by four nets a, b, c, and d. The 
potential face item Fl is a potential expansion of net d in the face 9300. For the face item Fl, the 
process 9100 would identify edge and node constraints for the layer on which this face item exists. 

In the embodiments described below, the edge and node constraints of a face item are the 
minimum required distances from the face item to the edges and nodes of the face that contains the 
face item. Specifically, for the face item, the constraint with respect to each node represents the 
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required width and spacing of the minimum set of topological items (if any) that separate the face 
item from the node. For the face item, the constraint with respect to each edge is the required width 
and spacing of the minimum set of topological items (if any) that topologically separate the face item 
from the edge. One manner of computing edge and node constraints is further described below. 

Figure 94 illustrates the edge and node constraints for the potential face item Fl. In this 
example, as well as the other examples illustrated in Figures 95-104, the acronyms can be 
understood as follows: (1) Wx signifies width of a route of a net x, (2) SNY-x signifies minimum 
required spacing between node Y and an adjacent net route x, (3) Sx-z signifies minimum required 
spacing between adjacent routes for nets x and z, (4) SFQ-x signifies minimum required spacing 
between face item FQ and adjacent net route x, and (5) SFQ-NY signifies minimum required spacing 
between node Y and adjacent face item FQ. 

As shown in Figure 94, the face item Fl's edge constraints for edges El and E3 are null, 
while its edge constraints for edge E2 is Wb + SFl-b. Also, this face item's (1) node constraint for 
node Nl is SF1-N1, (2) node constraint for node N2 is Wb + SFl-b + Wc + SN2-C + Sc-b, and 
(3) node constraint for node N3 is Wa + SN3-a + SFl-a. 

The process 9100 would identify the nodes N1-N3 of the face 9300 as the set of constraining 
points for the face item FL This process would also identify each node's constraint as the minimum 
required distance to the constraining point corresponding to the node. Figure 95 illustrates nodes 
N1-N3 as constraining points of the face item FL This figure also illustrates the minimum required 
distance to each of the points N1-N3 as a circular arc that has at its center one of the nodes N1-N3 
and has as its radius the node's corresponding node constraint. 

The process 9100 would also identify a legal region that can contain the face item Fl from 
this face item's edge constraints. This legal region can be obtained by bringing each edge of the face 
9300 towards its center in the direction of the edge's normal by the amount of the edge's constraint. 
In the example illustrated in Figures 93 and 94, the constraints for two edges El and E3 are 0, and 
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the constraint for the edge E2 is Wb + SFl-b. Accordingly, as illustrated in Figure 95, the face item 
Fl's legal region 9500 can be obtained by bringing edge E2 towards the center of the face 9300 in the 
direction of its normal by the amount Wb + SFl-b, and leaving the other two edges El and E3 
unadjusted. This legal region is a triangle defined by nodes Nl\ N2', and N3. As further described 
below, an optimizer can search through the legal region 9500 to identify the optimal position for the 
face itemFl. 

Each time the process 9100 first identifies the constraining points, distances, and/or legal 
region for a face item, it stores these values in the face items data structure. Also, each time the 
process makes a routing operation that changes the constraining distances and/or legal region of the 
face item, it modifies the stored values in the face item's data structure. Storing these parameters in 
the face item data structures speeds up the operation of the router in the instances where it can simply 
retrieve these values from the data structures. 

Constraints for Two Face It ems in the Same Face and Two Face Items Forming a Via. 
Figure 96 illustrates an example identical to the example illustrated in Figure 93, except that 
the face 9600 in Figure 96 includes two face items. One face item is face item Fl, which still 
represents the potential expansion of net d into face 9600. The other face item is a face item F2 that 
connects to (1) a face item F3 in a face 9700 (illustrated in Figure 97) in a layer above the face 
9600's layer, and (2) a face item F4 in a face 91030 (illustrated in Figure 98) in a layer below the 
face 9600's layer. 

The process 9100 would compute node and edge constraints for face items Fl and F2. Also, 
based on these constraints, the process would identify (1) sets of constraining points and distances, 
and (2) a legal region for each face item. The process would also identify a repulsion constraint that 
specifies the minimum spacing between the two face items. 

For the example illustrated in Figure 96, the node and edge constraints for the face item Fl 
are identical to the constraints shown for this face item in the example illustrated in Figure 93. In 
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other words, the face item F2 does not affect the node and edge constraints of the face item Fl. If the 
face item F2 connected to more than one edge of the face 9600, then it would affect the face item 
Fl's constraint for node N2. Figure 99 illustrates one such example. As shown in this figure, the 
connection of the face item F2 to edges E2 and E3 changes face item Fl's constraint for node N2. 
After identifying the node and edge constraints for the face item Fl, the process would identify the 
legal region that can contain this face item. As before, this legal region is the region 9500 illustrated 
in Figure 95. 

Figure 96 illustrates node and edge constraints for the face item F2 with respect to the nodes 
and edges of the face 9600. Like the constraints for the face item Fl, the route connected to the face 
item Fl does not affect the edge and node constraints for the face item F2. These constraints would 
have been affected if the face item Fl connected to more than one edge of the face 9600. 

The face item F2 forms a via with the face items F3 and F4 in the layers above and below. 
The constraints on all three face items F2, F3 and F4 constrain the location of this via. Accordingly, 
the process 9100 would identify node and edge constraints for face items F3 and F4 with respect to 
the nodes and edges of the faces that contain them. These constraints are illustrated in Figures 97 
and 98. 

For the via formed by these three face items, the process 9100 then would identify a legal 
region in each face 9600, 9700, and 91030. The legal region in each of these faces is obtained by 
moving each edge of the face towards the face's center by the amount of the edge's constraint. As the 
constraints of all the edges in faces 9700 and 91030 are null, the legal region for the face item F3 is 
the entire face 9700 and the legal region for the face item F4 is the entire face 91030. The legal 
region for the face item F2, however, is different than the face containing this face item. Figure 100 
illustrates this legal region 100, which is obtained by moving edge El towards the center of face 
9600 by the amount of its constraint Wb + SF2-b. 

Based on the identified legal region and node constraints for each face traversed by the via 
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that F2, F3, and F4 form, the process would specify a combined legal region and set of constraining 
points and distances for this via. The process would generate the combined region by intersecting the 
legal regions identified for face items F2, F3, and F4. Figure 101 illustrates this combined legal 
region 10100. The combined set of points and distances for this via are simply the union of the set of 
points and set of distances identified for the face items F2, F3, and F4. Figure 101 illustrates the 
combined sets of points, which are nodes N1-N9, and the combined sets of distances for these nodes. 

As mentioned above, the process 9100 would also identify a repulsion constraint that 
specifies the minimum required spacing between the face items Fl and F2. The required distance 
between two such face items is the sum of the width of each route between these face items plus the 
required minimum spacing between the route(s) and face items. In this example, the required 
distance between face items Fl and F2 is equal to SFl-b + Wb + SF2-b (i.e., is equal to the sum of 
the required spacing between Fl and net b, the width of net b, and the required spacing between F2 
and net b). 

As mentioned above, some embodiments allow a face to have at most two face items. These 
embodiments limit the number of face items to improve run-time efficiency. Other embodiments, 
however, might allow a face to have more than two face items. These embodiments would identify 
node, edge, and region constraints similarly to the examples described above. They would also 
analogously identify repulsion constraints between face items in a face except that, in a face with 
three or more face items, the required constraint between two face items would be the sum of the 
width of each route and/or face item between these two face items plus the required minimum 
spacing between the route(s) and face items. 

A face could also include a face item that is a Steiner point. Figure 102 illustrates a face 
10200 that is intersected by the routes of nets a, b, and c. This face includes two face items Fl and 
F2. Face item Fl is a Steiner point for a route of a net a, while face item F2 is a potential expansion 
of the route c into the face 10200. Figure 102 illustrates the node and edge constraints for the face 
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item Fl, while Figure 103 illustrates the node and edge constraints for the face item F2. 

Constraints for Face Item that Connects to a Node in a Different Laver. 

The process 9100 treats all vias that start or end with a node as fixed vias. This is the case 
even for the via nodes that are moveable by the routability engine, which is further described below. 
In other words, as far as the process 9100 is concerned, vias that start or end with nodes are fixed 
entities. 

Accordingly, when a via starts or ends with a node, the process 9100 does not specify node 
and edge constraints for the node in the face that contains the node. However, if the node forms a via 
by connecting to a face item above or below it, the process 9100 identifies constraining points and 
distances for the connected face item. The constraining points for such a face item are the nodes of 
the face containing the face item. The constraining distance to each node of the face is its node 
constraint. 

Figure 104 illustrates an example of a via formed between a node N4 and a face item F2 that 
are on two different layers. The constraining points for the face item F2 are the nodes Nl, N2, and 
N3. The constraining distance to each node Nl, N2, or N3 is a node constraint for the face item F2 
and the node. As illustrated in Figure 104, the face item F2's (1) distance constraint for node Nl is 
SF2-N1, (2) distance constraint for node N2 is Wc + SF2-C + SN2-C, and (3) distance constraint for 
node N3 is Wb + SF2-b + SN3-a + Sb-a + Wa. 

Constraints for Face Items in Faces Connected bv Vias to Faces Involved with an Expansion. 

A face item might exist on a face that connects through a via to a face involved with an 
expansion. Figure 105 illustrates an example of such a face item. This figure illustrates a potential 
expansion between face items Fl and F2. The faces involved with this expansion are faces 10505 
and 10510. Accordingly, the process 10500 would identify constraining points and distances for the 
face-item pair Fl and F2 and the face-item pair F3 and F4 on the faces 10505 and 10510. The 
process would also identify repulsion constraint specifying the minimum required distance between 
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face items F2 and F3. 

The face 10515 is a face that is not involved with the expansion but connects to face 10510 
(which is involved with the expansion) through the via formed by face items F3 and F4. The process 
9100 would treat face item F5 in face 10515 as a fixed face item. For the face item F4, it would 
compute a repulsion constraint that specifies the minimum required distance between the face items 
F4andF5. 

Identifying N ode and Hripe Constraints. 

In some embodiments, the process 9100 modifies (at 9110) the node and edge constraints for 
a face item in the following manner. The process traverses in a particular direction (e.g., 
counterclockwise) the edges of the face-item's face starting at the edge item on the edge connected to 
the face item. During its traversal of the edges, the process identifies the minimum required distance 
to the starting edge item for different points along its traversal. 

In some embodiments, the process uses a stack to keep track of the minimum required 
distances as it traverses the edges. In traversing the edges, the process might encounter routes for 
three types of other nets. The first type does not connect to a face item within the face. The second 
type connects to a face item with a degree greater than one (Le., connects to a face item that connects 
to more than one topological item on the face boundary). The third type connects to a face item with 
a degree one (Le., connects to a face item that only connects to one topological item on the face 
boundary) 

The process adds to the stack the required width plus spacing distances of the encountered 
routes for the first- and second-type nets when it encounters these routes for the first time while 
traversing the identified face's edges. The process removes from the stack the distances of each such 
route when it encounters the route for the last time while traversing the identified face's edges. When 
the process encounters a route for a third-type net while traversing the identified face's edges, the 
process does not add any space and width distance to the stack. 
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Additionally, when the process encounters either a second- or third-type route while 
traversing the edges, the process identifies and stores the required repulsion distance between the 
face item that it is currently examining and the face item that is connected to the encountered second- 
or third-type route. 

The process 9100 sets the constraint for a node when it reaches the node while traversing the 
boundary of the identified face. The constraint for a node that does not connect to any route is the 
current value at the top of the stack plus the minimum required distance between the last edge item 
encountered and the node. 

Identifying the constraint for a node that connects to a route is slightly more involved. When 
the process encounters such a node, it determines whether it is the last time that it will encounter the 
route connected to the node before completing its loop traversal about the identified faces edges. 

If so, this node constraint is (1) the value on the stack after the removal from the stack of the 
spacing and width distances of net route that connects to the node, plus (2) the minimum required 
distance between the last edge item encountered and the node. If not, this node constraint is (1) the 
value on the stack before the addition to the stack of the spacing and width distances of net route that 
connects to the node plus (2) the minimum required distance between the last edge item encountered 
and the node. 

Each time the process 9100 reaches a node, it also sets the constraint for the edge that it was 
just traversing to reach that node. This edge constraint is the smallest value that appeared at the top 
of the stack since the last node. 

After identifying constraints at 9110, the process 9100 determines (at 9115) whether it needs 
to run an optimization process to identify the location of any face item within a face that is involved 
with the path expansion. The process needs to run an optimization process when at least one non-via 
face item, or at least one set of via face items, for which it generated a constraint is moveable within 
a legal region that is larger than a threshold amount. 
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If the process determines (at 9115) that it does not need to use the optimization process, it 
determines (at 9120) whether it identified at 91 10 a legal region that was smaller than the threshold 
amount. If so, the process transitions to 9140, which is further described below. If not, the process 
transitions to 9135. For each face item that the process identified a set of constraints at 9110, the 
5 process determines (at 9135) whether the face item satisfies its set of constraints. If the process 
determines (at 9135) that all constraints are not met, it specifies (at 9140) the expansion as illegal, 
and then ends. On the other hand, if the process determines (at 9135) that all constraints are met, it 
specifies (at 9145) the expansion as legal, and then ends. 

If the process determines (at 91 15) that for at least one face item it has defined a region in 
3 which the face item can move, the process formulates (at 9125) a function to optimize. This function 
can include one or more sub-functions in some embodiments. Several different types of sub- 
functions are described below. 

Sub-function for Non-Via Face Ttem jn „ F ace Affect hv an Pv p on»„ 

For each non-via face item in a face involved with the expansion, the process 9100 specifies 

one sub-function based on the point and distance constraints that it identified (at 9110) for the face 

item. The sub-function takes the form of: 

p r 

P=l a p 

where p is a variable that represents a point in the set of constraining points for the face item, P is the 
final point in the set, d„ represents the distance from a variable x,y location of the face item to the 
current point p, and r p represents the identified minimum constraining distance to the point p. This 

sub-function varies with the x,y location of the face item in the region identified for it at 9110. 

Region 9500 of Figure 95 is an example of such a region for a non-via face item in a face affected 

by the expansion. As further described below, the optimizer searches through this region and 

specifies an x,y location of the face item. 
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Sub-function for Via- Forming Face Item in a Face Affected bv an Expansion. 
For a face item that forms a via with a node (/.<?., a face item that has a fixed position for the 
process 9100), the process 9100 identifies a sub-function that has the same form: 

p r 

p=i u p 

5 where p is a variable that represents a point in the set of constraining points for the face item, P is the 
final point in the set, d p represents the distance from the fixed x,y location of the face item to the 
current point p, and r p represents the identified minimum constraining distance to the point p. This 
sub-function is a constant because its x,y location is fixed, and so is its distance d p from each of its 
constraining points. Face item F2 of Figure 104 is an example of such a fixed face item. When a via 
10 is formed by first and second face items in first and second layers and a node in a third layer, the sub- 
function is similar to the one described above, except that the set of constraining points and distances 
include points and distances that constrain both face items. 

For a set of face items that form a moveable via, the sub-function again has the form: 

f't^- (6) 

15 In this equation, p is a variable that represents a point in the set of constraining points for all the face 
items, P is the final point in the set, and r„ represents the identified minimum constraining distance to 
the point p. In this sub-function, d p represents the distance from the x,y location of the current point p 
to a variable x,y location of the via (i.e., of the set of face items) within its region, which is an 
intersection of the legal region of each of the via's face items. Region 10100 of Figure 101 is an 

20 example of such a region for a via that is formed by the face items F2, F3, and F4 of faces 9600, 
9700, and 91030 of Figures 96, 97, and 98. As further described below, the optimizer searches 
through this region 10100 and specifies an x,y location for the via formed by the face items F2, F3, 
and F4. 
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Sub-function for Two face items in the same face. 

When a face examined by the process 9100 includes two face items, the process specifies a 
sub-function for expressing the repulsion constraint between the two face items. This sub-function 
//,2 is as follows, 

/1.2- — , (7) 

"l/o2 

where d„ o2 represents the distance from the x,y location of face item 1 to the x,y location of face item 
2, and ri t0 2 represents the identified minimum required distance between the two face items. 

When both face items are moveable, the x,y location for both face items can vary within 
regions defined for them. For the two face items F2 and F3 illustrated in Figure 105, Equation (7) 

10 above can be used to specify a sub-function / 2 ,. ? . In this example, the face item F2 can be at any x.y 
location within the region that is obtained by intersecting the legal region for face items Fl and F2, 
while the face item F3 can be at any x,y location within the region that is obtained by intersecting the 
legal region for face items F3 and F4. Also, in this example, r 2to3 represents the identified minimum 
required distance between face items F2 and F3. 

15 In Equation (7), a face item's x,y location is fixed when it forms a via with a node. For 

example, for the two face items Fl and F2 illustrated in Figure 104, Equation (7) above can be used 
to specify a sub-function f u . In this example, the face item Fl can be at any x,y location within its 
legal region, while the face item F2 has a fixed x,y location defined as the x,y location of node N4. In 
this example, ri, o2 represents the identified minimum required distance between face items Fl and 
20 F2. 

In Equation (7), a face item's x,y location is also fixed when the face item is on a face that 
(1) is not involved with a path expansion but (2) connects through a via formed by another face item 
to a face that is involved with an expansion. For instance, face item F5 in Figure 105 illustrates one 
such face item. This face item F5 is on a face 10515 that is not involved with the path expansion 
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from Fl to F2 but connects to face , which is involved with this path expansion. The process 9100 
would use Equation (7) to specify a sub-function f 4 . 5 , which expresses the repulsion constraint 
between face items F4 and F5. In this example, the face item F4 can be at any x,y location within the 
region that is the intersection of the legal region for face item F3 and F4, while the face item F5 has a 
fixed x,y location that is its previously defined position. In this example, r 4 ,o5 represents the 
identified minimum required distance between face items F4 and F5. 

At 9125, the process 9100 formulates the cost function as a sum of all sub-functions that it 
specifies to express the location and repulsion constraints for the face items in the faces involved 
with the path expansion. For instance, the process specifies the following function for the example 
illustrated in Figure 105. 

f B f,+fs+fu + fv 

In this function, (1)// is a sub-function according to Equation (6) for specifying a cost based on the 
location of the via formed by face items Fl and F2, (2)/ ? is a sub-function according to Equation (6) 
for specifying a cost based on the location of the via formed by face items F3 and F4, (3) f 2<3 is a sub- 
function according to Equation (7) for specifying a cost based on the repulsion between face items 
F2 and F3, and (4)f 4i5 is a sub-function according to Equation (7) for specifying a cost based on the 
repulsion between face items F4 and F5. 

The function formulated at 9125 can be just one sub-function. For instance, the function is 
just the sub-function express by Equation (7) when the path expands to a non-via face item in a face 
with only one face item. This is the case for the example illustrated in Figure 93. Also, the cost- 
function formulation approach described above works for the embodiments that allow more than two 
face items in a face. The cost function in these embodiments would typically have more sub- 
functions as there are more location and repulsion constraints to consider. 

After formulating the cost function at 9125, the process 9100 (at 9130) has a minimizer 

search through the legal region of each moveable non-via face item or set of via-forming face items 
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to identify an x,y location for each face item that minimizes the function. In some embodiments, the 
process 9100 uses the method of steepest descent for the minimization operation. This method is 
discussed in Numerical Recipes in C: The Art of Scientific Computing, by William Press, et al., 
Second Edition (1992). 

After performing the minimization operation, the process 9100 determines (at 9135) whether 
the location identified for each face item satisfies the point, distance, and/or repulsion constraints 
identified for it. Specifically, the process determines whether each face item is away from each of its 
constraining points by at least the required minimum distance for that point. Also, if the face item is 
in a face with another face item, the process determines whether the face items are away from each 
other by at least their minimum required distance. If a face item does not satisfy its distance or 
repulsion constraint, then the process 9100 transitions to 9140 to specify as illegal the expansion that 
triggered the via check. Otherwise, the process specifies (at 9145) the expansion as legal. 
(2nd) Non-Optimization Technique 
Figure 92 illustrates a process 9250 that performs a simpler via- checking operation than the 
process 9100 of Figure 91. Some embodiments that use this process limit each face to have at most 
one face item. This process initially determines (at 9255) whether the face that contains the 
destination node, face item, or edge item has a face item within it. When the destination item is part 
of several faces, some embodiments execute the process 9250 for each face, and specify an 
expansion as illegal when the expansion fails the via-legality check for just one of these faces. On the 
other hand, like the process 9100 of Figure 91A, other embodiments might not examine each face 
that contains the destination item, but rather might only examine the expansion face as defined above 
for 9105. 

If the process 9250 determines (at 9255) that a face that contains the destination item does 
not contain a face item, the process ends as a via-check is unnecessary. Otherwise, the process 
computes (at 9260) the bounding polygon of the face item. Some embodiments compute this polygon 
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by computing the closest point to each node on each edge of the face that the face item can be (i.e., 
the node, offset along the edge by the flow on that edge between the node and the face item). 

Next, the process computes (at 9265) the intersection of the face item's bounding polygon 
with via destination (if any) in the layer above or below. The process then examines (at 9270) the 
size of the intersection to determine whether there is sufficient space for the face item. If the 
intersection is too small, then the process determines (at 9275) that the expansion is not legal, resets 
(if necessary) the face-item's bounding polygon to its original shape before 9260, and then ends. If 
the intersection is big enough for the face-item, the process indicates (at 9280) via-check legality, 
and then ends. 

(iii) Planarity Check 

As indicated in Figure 90, the planarity check is performed for four of the expansions. For 
the node to node expansion, the planarity check simply ensures that there is no route crossing the 
shared edge between two nodes. All other planarity checks are performed by traversing the linked list 
of edges from the source towards the destination node/edge, checking that no other net's route blocks 
a direct connection between the source and destination items of the path expansion. 
(6) Cheapest Path Calculation 
Figure 106 illustrates a process for computing the cheapest-route cost for a net. As shown in 
Figure 106, the process 10600 initially identifies (at 10602) one or more pairs of pins for the net. It 
then selects (at 10604) a pair of pins for the net. The process next identifies (at 10606) a list of 
source and target port-geometries for the selected pair of pins of the net. 

The process then identifies (at 10608) a Depth_Limit for generating routing path between the 
source and target geometries. In the embodiments described below, the Depth_Limit is set to the 
shortest distance between the closest source and target geometries. After 10608, the process 
generates (at 10610) all legal paths between the source and target geometries identified at 10606. At 
10610, the process 10600 can use a path-generating process similar to the one described above by 
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reference to Figure 84. 

After 10610, the process determines (at 10612) whether it was able to generate (at 10610) any 
legal path between the source and target geometries. If not, the process increments (at 10614) the 
Depth_Limit. The process then determines whether the Depth JLimit exceeds a maximum depth 
limit. If so, the process indicates (at 10618) failure to route the net. Otherwise, the process transitions 
to 10610 to generate legal paths between the source and target geometries for the Depth_Limit 
incremented at 10614. 

If the process determines (at 10612) that it was able to generate (at 10610) a legal path 
between the source and target geometries, the process records (at 10620) the generated path or paths 
and the cost for each path. The process then determines (at 10622) whether it has examined the last 
pin-pair identified at 10602. If not, the process transitions back to 10604 to select the next pin-pair. 

On the other hand, if the process determines that it has examined all the pin-pairs, the process 
(at 10624) identifies the cheapest routes for the net, and records the costs of this route. When the net 
has more than two pins in the sub-region, the process 10600 generates more than one set of paths for 
more than one set of pin-pairs. In such a situation, the process has to identify (at 10624) the cheapest 
combination of paths for the different pin-pairs to generate the cheapest route. The process 10600 
identifies the cheapest combination of paths in the same manner as described above for process 8200 
at 8232. After 10632, the process ends. 
IV. THE COMPUTER SYSTEM 

Figure 107 presents a computer system with which one embodiment of the present invention 
is implemented. Computer system 10700 includes a bus 10705, a processor 10710, a system memory 
10715, a read-only memory 10720, a permanent storage device 10725, input devices 10730, and 
output devices 10735. 

The bus 10705 collectively represents all system, peripheral, and chipset buses that 
communicatively connect the numerous internal devices of the computer system 10700. For instance, 
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the bus 10705 communicatively connects the processor 10710 with the read-only memory 10720, the 
system memory 10715, and the permanent storage device 10725. 

From these various memory units, the processor 10710 retrieves instructions to execute and 
data to process in order to execute the processes of the invention. The read-only-memory (ROM) 
10720 stores static data and instructions that are needed by the processor 10710 and other modules of 
the computer system. The permanent storage device 10725, on the other hand, is read-and-write 
memory device. This device is a non- volatile memory unit that stores instruction and data even when 
the computer system 10700 is off. Some embodiments of the invention use a mass-storage device 
(such as a magnetic or optical disk and its corresponding disk drive) as the permanent storage device 
10725. Other embodiments use a removable storage device (such as a floppy disk or zip® disk, and 
its corresponding disk drive) as the permanent storage device. 

Like the permanent storage device 10725, the system memory 10715 is a read-and-write 
memory device. However, unlike storage device 10725, the system memory is a volatile read-and- 
write memory, such as a random access memory. The system memory stores some of the instructions 
and data that the processor needs at runtime. In some embodiments, the invention's processes are 
stored in the system memory 10715, the permanent storage device 10725, and/or the read-only 
memory 10720. 

The bus 10705 also connects to the input and output devices 10730 and 10735. The input 
devices enable the user to communicate information and select commands to the computer system. 
The input devices 10730 include alphanumeric keyboards and cursor-controllers. The output devices 
10735 display images generated by the computer system. For instance, these devices display IC 
design layouts. The output devices include printers and display devices, such as cathode ray tubes 
(CRT) or liquid crystal displays (LCD). 

Finally, as shown in Figure 107, bus 10705 also couples computer 10700 to a network 10765 
through a network adapter (not shown). In this manner, the computer can be a part of a network of 
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computers (such as a local area network ("LAN"), a wide area network ("WAN"), or an Intranet) or a 
network of networks (such as the Internet). Any or all of the components of computer system 10700 
may be used in conjunction with the invention. However, one of ordinary skill in the art would 
appreciate that any other system configuration may also be used in conjunction with the present 
invention. 

While the invention has been described with reference to numerous specific details, one of 
ordinary skill in the art will recognize that the invention can be embodied in other specific forms 
without departing from the spirit of the invention. Thus, one of ordinary skill in the art would 
understand that the invention is not to be limited by the foregoing illustrative details, but rather is to 
be defined by the appended claims. 
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CLAIMS 



We claim: 



1. 



A method of routing a set of nets, the method comprising: 



a) 



specifying a first order for the set of nets; 



b) 



routing the nets according to the specified first order; 



c) 



specifying a second order for the set of nets, wherein the second order has the 



fewest possible number of differences with the first order; 

d) routing the nets according to the specified second order. 

2. The method of claim 1 further comprising: 

a) specifying a third order for the set of nets, wherein the third order has the next 
fewest possible number of differences with the first order; 

b) routing the nets according to the specified third order. 

3. The method of claim 2 further comprising, after routing the nets according to the 
second order but before routing the nets according to the third specified order, identifying a set of 
orders for the set of nets, wherein each order in the identified set of orders has the same number of 
differences with the first order as the second order. 

4. The method of claim 3 further comprising: 

routing the nets according to each of the orders in the identified set of orders. 

5. The method of claim 4, wherein identifying the set of orders comprises iteratively 
identifying each order in the set and routing the nets at each particular iteration according to the order 
for the particular iteration. 
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6. The method of claim 5, wherein the routing the nets according to at least some of the 
orders results in a complete routing solution that includes a route for each net, the method further 
comprising: 

after at least each particular routing that result in a complete solution, 
5 determining whether the particular routing solution is better than any of the prior complete routing 
solutions examined, and if so, storing the particular routing solution. 

7. The method of claim 6, wherein the set of orders is a first set of orders, the method 
further comprising: 

identifying a second set of orders for the set of nets, wherein each order in the 
10 identified second set of orders has the same number of differences with the first order as the third 
order; 

routing the nets according to each of the orders in the identified second set of 

orders. 

8. The method of claim 7 further comprising: 

15 identifying a third set of orders for the set of nets, wherein each order in the 

identified second set of orders has one more difference with the first order than the number of 
differences between the first and third orders; 

routing the nets according to each of the orders in the identified third set of 

orders. 

20 9. The method of claim 8 further comprising: 

iteratively identifying additional sets of orders for the set of nets, wherein at 
each iteration the number of differences between the orders in the set and the first order is more than 
the number of differences between the orders in the set from the previous iteration and the first order; 
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at each iteration, identifying a routing solution for at least one order in the 

iteration; 

terminating the iterative identification of additional order sets once the number 
of differences between the orders in the next iteration and the first order reaches a certain value. 

10. The method of claim 8 further comprising: 

iteratively identifying additional sets of orders for the set of nets, wherein at 
each iteration the number of differences between the orders in the set and the first order is more than 
the number of differences between the orders in the set from the previous iteration and the first order; 

at each iteration, identifying a routing solution for at least one order in the 

iteration; 

terminating the iterative identification of additional order sets once the number 
of differences between the orders in the last iteration and the first order reaches a certain value. 

11. A computer readable medium that stores a computer program for routing a set of nets, 
the computer program comprising sets of instructions for: 

a) specifying a first order for the set of nets; 

b) routing .the nets according to the specified first order; 

c) specifying a second order for the set of nets, wherein the second order has the 
fewest possible number of differences with the first order; 

d) routing the nets according to the specified second order. 

12. The computer readable medium of claim 11, wherein the computer program further 
comprises sets of instructions for: 
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a) specifying a third order for the set of nets, wherein the third order has the next 
fewest possible number of differences with the first order; 

b) routing the nets according to the specified third order. 

13. The computer readable medium of claim 12, wherein the computer program further 
comprises instructions for identifying, after routing the nets according to the second order but before 
routing the nets according to the third specified order, a set of orders for the set of nets, wherein each 
order in the identified set of orders has the same number of differences with the first order as the 
second order. 

14. The computer readable medium of claim 13, wherein the computer program further 
comprises a set of instructions for routing the nets according to each of the orders in the identified set 
of orders. 

15. The computer readable medium of claim 14, wherein the set of instructions for 
identifying the set of orders comprises a set of instructions for iteratively identifying each order in 
the set and routing the nets at each particular iteration according to the order for the particular 
iteration. 

16. The computer readable medium of claim 15, wherein the routing of the nets according 
to at least some of the orders results in a complete routing solution that includes a route for each net, 
the computer program further comprises instructions for determining, after at least each particular 
routing that result in a complete solution, whether the particular routing solution is better than any of 
the prior complete routing solutions examined, and if so, storing the particular routing solution. 

17. The computer readable medium of claim 16, wherein the set of orders is a first set of 
orders, the computer program further comprising sets of instructions for: 
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' I 1 t 

identifying a second set of orders for the set of nets, wherein each order in the 
identified second set of orders has the same number of differences with the first order as the third 
order; 

routing the nets according to each of the orders in the identified second set of 

5 orders. 

18. The computer readable medium of claim 17, wherein the computer program further 
comprises sets of instructions for: 

identifying a third set of orders for the set of nets, wherein each order in the 
identified second set of orders has one more difference with the first order than the number of 
10 differences between the first and third orders; 

routing the nets according to each of the orders in the identified third set of 

orders. 

19. The computer readable medium of claim 18, wherein the computer program further 
comprises sets of instructions for: 

15 iteratively identifying additional sets of orders for the set of nets, wherein at 

each iteration the number of differences between the orders in the set and the first order is more than 
the number of differences between the orders in the set from the previous iteration and the first order; 

at each iteration, identifying a routing solution for at least one order in the 

iteration; 

20 terminating the iterative identification of additional order sets once the number 

of differences between the orders in the next iteration and the first order reaches a certain value. 

20. The computer readable medium of claim 18, wherein the computer program further 
comprises sets of instructions for: 
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iteratively identifying additional sets of orders for the set of nets, wherein at 
each iteration the number of differences between the orders in the set and the first order is more than 
the number of differences between the orders in the set from the previous iteration and the first order; 

at each iteration, identifying a routing solution for at least one order in the 

5 iteration; 

terminating the iterative identification of additional order sets once the number 
of differences between the orders in the last iteration and the first order reaches a certain value. 
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ABSTRACT OF THE INVENTION 

Some embodiments of the invention provide a method of routing a set of nets. The method 
specifies a first order for the set of nets. It then routes the nets according to the specified first order. 
The method then specifies a second order for the set of nets, where the second order has the fewest 
possible number of differences with the first order. The method then routes the nets according to the 
specified second order. 
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Identify drop refered to by the 
Current_Drop f s back reference. 
Specify the identified drop as 
the CurrentJDrop 



Figure 71 



Face 

-Reference to 3 edges 
-Reference to 3 nodes 
-Up to two references for up to two 
face item 



7200 



Figure 72 



Node 



-Net Identifier 

-One or more planar-path references to 
adjacent topological items in the same 
planar path 

-A pair of via-path references to up and 
down topological via items 
-A references to list of edges connected 
to the node 

-For each edge, an edge reference to 
the next or previous topological item 
on the edge 

-A reference to the geometry of the 
node 

-Vertex number identifying the vertex 
of the geometry 
-Location of the node 



Figure 74 



7301 



7400 



Edge 

-Two references for up to two faces of 

the edge 

-Capacity 

-Flow 

-Constrained 

-Linked list of items on the edge 
starting with one of the edge's nodes 
and ending with its other node 



Figure 73 



7500 



Edge Item 



-Reference to its edge 
-Net Identifier 

-A pair of planar-path references to 
adjacent topological items in the same 
planar path 

-A pair of edge references to the next 
and previous topological item on the 
edge 



Figure 75 



7600 



Face Item 



-Reference to its face 
-Net Identifier 

-Up to 3 planar-path references for 
adjacent topological items in the same 
planar path 

-A pair of via-path references for up 
and down topological via items 
-Bounding polygon that defines legal 
face item locations 
-Constraining Points and Distances 



Figure 76 



Q Start ^ 



Triangulate sub-region 



7905 



7900 



Group nets in the sub-region 



7910 

V 



7915 



Select a group of nets 



HFrom 7980 



Specify a depth limit, a maximum depth-limit threshold, 
and max pushes for the group of nets 



7920 
< 



7925 
/ 



From 7955 
or 7960 



Call solving engine to find topological solution for all or some of the 
nets in the group within the specified depth limit and max pushes. 



No 



7930 

Return 
'solution for all nets, 
.m group?. 

No 



Yes 



Increment depth limit 



V 



7935 




7940 



Fivuve 79- Fisme 79A 

figure 79. Figure 7 9B 



To 7965 



Figure 79 A 



From 7940 



i ' 7945 




To 7920 

t 





Accept returned best 


No 


solution 



Pick a net and create solution ^ 
with just its cheapeast path 



7950 



7955 



7 50 

"All netsV No 
in group 
vSolved?> 



Yes 



Assign face item 
shapes (an d locations) 

i 



-4 



7965 



From 7930 



7970 



Triangulate based on face items, 
and modify route descriptions 



V 



7975 



Perform edge-flipping and 
modify route description 



V 



Yes 



To 7915 



Any net left? 



y 



7980 



No 



Figure 79B 



8100 



Q Start ^ 



Current_Net is first net in the received net list. 
Push Count = 0 



8102 
V 



8104 



For Current_Net, produce all legal routes that 
have an acceptable cost 




-From 8181 



To 8138 



8108 



Store returned solutions 




8110 



Select route for Current Net 



-From 8136 



Increment Push by 1 

l 

To 8114 

Figure 81A 



8112 

V 



Figure 81' Fisure 81A 
riguredi. Figure glR 

Figure 81C 



1 



Insert selected route in region. 



8114 

V 



Compute increase in total cost 
due to addition on of selected route 



8116 



V 




To 8132 



8122 



Store solution 



Yes 



To 8104 



1 



8128 



Set Current_Net to be 
the next net in the net list 



Return best solution seen 
and identify nets solved 



-From 8132 



Q END ^ 



Figure 81B 
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8132 




Return failure to route 



^ END ^ 



( Stol ) 



8200 
/ 





r 


8202 


Break net into one or more pairs of pins 


/ 




r 


8204 


; Select a pair of pins 




■< 


\ 


r 


8206 


Determine list of source and targets for the 
selected pair 


/ 


\ 


r 


8208 


Set Depth JLimit to straight line distance 
between closest source and target 






r 


8210 


Set Solution_Count to 0 






8212 


Generate al 


legal paths 




for the selected pin-pair 







r 





From 8224 



From 8218 



To 8214 



Figure 82 A 



Figure 82- Figure 82 A 
figure dJ. Figure82B 






No 8216 






Increment Depth_Limit 



To 8212-*- 



8222 



Return failure 
to find paths 



To No 
8204 M 



Identify K cheapest routes 

i 



Pass through generated routes 
and create Steiner points whenever possible 


\ 





8234 

V 



•( END ) 



Figure 82B 



• 1 



From 8416 
and 8434 " 



From 8428 



8400 



Set one of the source nodes as starting point of path, define 
Current Point to be this source node, and set Path_Cost to 0. 



8402 
V 



For Current_Point, generate all expansions 



8404 
V 



Select one of the expansions 



8406 
V 



Expansion\ >J k_^ To842g 
legal? 



Yes 



8410 





Define CurrentJPoint to be the selected expansion 
point, and add Current Point to path 








f 




Calculate cost due to the expansion and revise Path_Cost of the path. 



8412 
V 



T 



To 8414 



Figure 84A 



Figure 84A 
Fl S ure 84: Fi2ure 84B 




Yes 



8422 



Remove Current_Point from path 



8424 



Set Current_Point to last topological item in path 



To 3404 
No 



Is 

Current_Point a 
Target? 



8416 



Yes 



8420 
/ 



Record path with cost, 

and increment 
Solution_Count by 1 



8406 



No 



No 



Adjust 
Path Cost 



8426 

V 



8428 




From 8408 



8430 



8404 



Set one of the source nodes as 
starting point of path, define 
CurrentJPoint to be this source 
node, and set Path Cost to 0. 



8434 



Figure 84B 




9250 
/ 



Compute bounding polygon of face item 




r 



9260 



9265 



with bounding polygons above and below 



V 




9275 



\ 


Yes 

r 


Legal 




f 



^ END ^ 



\ No „ 


Not legal; reset 
bounding polygon 




9280 






/ 






U 







Figure 92 



5 



Start ^ 



Break net into one or more pairs of pins 



10600 
/ 

10602 
V 



10604 



Select a pair of pins 



Y 



Determine list of source and targets for the 
selected pair 



10606 
V 



10608 



Set Depth_Limit to straight line distance 
between closest source and target 



10610 



Generate all legal paths 
for the selected pin-pair 



Yes 



Record solutions 



10612 

Is there 
any path for current 
epthJLimit? 

106^ 

V 



No 



1 



10614 



Increment Depth_Limit 



No 




10616 

y 

'Max depth 

limit? y jsj 0 



If multiple pirs, generate cheapest path 



10624 



Yes 



10618 



^ END ^ 





Return failure to find 




paths 


U 





Figure 106 
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